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SECTION  I 
INTRODUCTION 


The  objective  of  this  research  was  to  develop  experimental  techniques  for 
the  quantitative  characterization  of  adhesively  bonded  multilayer  struc¬ 
tures  nondestructively .  These  multilayered  structures  are  very  similar 
to  the  ones  used  in  mechanically  fastened  aircraft  wings. 

Concerning  the  structural  integrity  of  the  wings,  one  is  interested  in 
characterizing  the  inner  layer  by  measuring  the  acoustic  parameters  and 
studying  these  parameters  as  a  function  of  aging.  The  acoustic  reflec¬ 
tion  and  transmission  coefficients  for  a  layered  structure  can  be  computed 
using  the  theoretical  model  of  Brekhovskikh.  Details  of  this  model  and 
the  computer  program  for  calculation  of  reflection  and  transmission  coeffi¬ 
cients  for  layered  structures  developed  during  the  initial  phase  of  this 
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research  program  were  presented  in  an  intermediate  report.  Based  upon 
the  theoretical  model,  acoustic  parameters  such  as  attentuation  coeffi¬ 
cients  are  determined  from  the  experimental  data  using  a  fit  program 
which  has  been  developed  further  since  the  intermediate  report  was  pub¬ 
lished.  The  determination  of  these  acoustic  parameters  using  cw  spec¬ 
troscopy  is  described  in  Section  II. 

Microcracks  in  wing  structures  grow  and  develop  because  of  tremendous  stres 
which  these  structures  experience  in  flight.  The  types  of  flaws  most  fre¬ 
quently  produced  in  aircraft  wing  structures  are  radial  cracks  emanating  f 
the  bolt  hole  either  at  the  faying  surface  of  the  first  layer  or  at  th 
inner  layer.  These  cracks  are  caused  by  fatigue  and  develop  from  the  sur¬ 
face  down.  Structural  components  can  continue  in  service  even  with  a 
crack,  provided  its  dimensions  remain  below  certain  critical  values.  In 
view  of  this  fact,  it  becomes  very  important  to  be  able  to  measure  periodi¬ 
cally  the  size  of  these  radial  cracks  nondestructively. 
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Several  techniques  such  as  amplitude,  time,  and  spectral  methods  can  be  used 
to  measure  the  radial  dimensions  of  these  cracks.  Since  each  crack  has  a 
different  reflectivity,  quite  Inconclusive  Information  regarding  size  is 
extracted  by  comparing  the  amplitude  of  the  backseat tered  echo  with  the 
signal  from  a  reference  reflector  and  thus  determining  the  equivalent  dimen¬ 
sion.  The  amplitude  methods  disregard  several  factors  such  as  orientation 
of  the  crack,  shape,  width  of  the  crack  opening,  and  quality  of  the  coupling — 
all  of  which  affect  measurement  accuracy  drastically.  Timing  methods  measure 
the  crack  dimension  from  the  time  delay  of  a  signal  returning  from  a  crack 
relative  to  the  instant  it  left  the  transducer  or  some  known  scattering 
point  on  the  sample.  Although  these  methods  yield  better  results  with  a 
focused  beam,  they  are  quite  cunfrersome  since  the  best  results  are  obtained 
only  if  the  beam  enters  the  specimen  in  the  plane  of  the  crack — which  is  not 
possible  for  the  case  of  aircraft  wing  geometry. 

Due  to  the  particular  geometry  of  radial  bolt-hole  cracks,  a  novel  method 
based  upon  the  Doppler  shift  was  developed  to  measure  the  radial  crack 
size.  This  method  does  not  depend  upon  the  amplitude  or  time  delay  of 
the  reflected  signal  from  the  crack.  The  radial  size  of  the  crack  is 
proportional  to  the  maximum  Doppler  shift  of  the  signal  backscattered 
from  the  tip  of  the  crack.  The  method  is  far  more  accurate  than  any 
previously  used  method  and  has  potential  for  practical  use  in  the  field. 
Section  III  describes  the  measurement  of  radial  cracks  using  the  Doppler 
effect. 

For  adhesively  bonded  layered  structures,  the  attenuation  a  is  very  large. 
Consequently,  the  sensitivity  of  the  sampled  cw  technique  is  reduced 
drastically.  The  conventional  pulse-echo  technique  cannot  be  used  for 
such  specimens  since  the  received  signal  becomes  comparable  to  the  thermal 
noise  of  the  receiver.  An  alternative  approach  to  solving  the  above  S/N 
problem  was  undertaken. 

Random-signal-correlation  techniques  can  be  used  in  ultrasonic  flaw- 
detection  systems  to  overcome  several  disadvantages  of  conventional 
pulse-echo  systems.  In  a  conventional  pulse-echo  system,  a  high  ratio 


2 


of  peak-to-average  transmitted  power  is  usually  required  to  meet  the  range 
and  resolution  requirements  of  the  system.  Furthermore,  when  strongly 
sound-absorbent  material  is  present,  a  pulse-echo  system  requires  large 
average  transmitted  power  to  maintain  an  acceptable  S/N  ratio.  However, 
ultrasonic  transducers  are  limited  in  average  power-handling  ability  due 
to  heating  effects  and  are  limited  in  peak-power  input  capability  due  to 
voltage-breakdown  effects.  An  additional  limitation  of  conventional  pulse- 
echo  systems  is  that  they  offer  no  improvement  in  the  S/N  ratio  of  the 
received  signal. 

The  inherent  difficulties  with  the  conventional  pulse-echo  system  mentioned 

above  were  very  successfully  overcome  in  a  laboratory  system  developed  by 
3 

Ferguson,  et^  al.  This  system  transmits  a  long  burst  of  ultrasonic  noise 
and  performs  a  correlation  of  the  received  echoes  with  a  delayed  version 
of  the  transmitted  burst.  In  this  system  the  resolution  does  not  depend 
upon  the  transmitted  pulse  width,  as  in  conventional  pulse-echo  systems, 
but  upon  the  bandwidth  of  the  transmitted  noise.  In  addition,  an  improved 
S/N  ratio  of  the  received  echo  can  be  achieved  (at  the  expense  of  a  longer 
data-acquisltion  time)  by  scanning  the  signal  slowly. 

From  a  practical  point  of  view,  the  system  just  described  has  a  serious 
drawback.  The  water  delay  line  which  was  used  consisted  of  two  matched 
ultrasonic  transducers  in  a  water  bath,  and  delay  was  varied  using  a 
micrometer  drive.  Even  if  this  delay  line  were  automated  and  compute 
controlled,  it  would  be  bulky  and  require  two  closely  matched  transducers. 
These  two  transducers,  along  with  the  transmitting  and  receiving  trans¬ 
ducers,  contribute  to  the  overall  system  response,  making  it  difficult 
to  isolate  the  desired  response  due  to  the  flaw  along. 

An  ultrasonic  random-signal-correlation  system  which  does  not  employ  a 
delay  line  and  which  may  be  used  with  a  single  transducer  has  been 
designed  and  built.  Section  IV  describes  this  Ultrasonic  Pseudorandom- 
Signal-Correlation  System. 
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Section  V  contains  conclusions  on  the  progress  reported  herein.  Suggestions 
are  also  made  for  further  efforts  which  could  be  initiated  utilizing  the 
results  reported  here. 
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SECTION  II 

CW  ULTRASONIC  SPECTROSCOPY 


Current  Air  Force  requirements  include  the  ability  to  detect  delamination 
and  degradation  of  adhesive  bonds  or  sealants  in  multilayered,  adhesively 
bonded  aircraft  structures.  For  nondestructive  evaluation,  ultrasonic 
examination  provides  the  only  access  to  the  inner  layers  of  such  sandwich¬ 
like  structures.  Because  the  ultrasonic  signal  will  be  affected  by  each 
material  layer  through  which  it  travels,  some  means  must  be  utilized  to 
extract  information  concerning  the  layer  of  interest.  Ultrasonic  spectros¬ 
copy  provides  such  a  method  of  isolating  the  contribution  of  each  material 

layer  to  the  total  frequency  spectrum  of  the  multilayered  structure.  Ultra- 
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sonic  spectroscopy  has  been  proposed  as  a  means  for  determining  grain-size, 
for  measuring  plate  thicknesses  or  sound  velocities  for  detecting  delami¬ 
nations  in  bonded  structures,^  for  classifying  bond  strength, ^  and  for 

8 

measuring  attentuation. 

A  method  has  been  developed  which  makes  use  of  cw  ultrasonic  spectroscopy 

in  the  determination  of  the  acoustic  parameters  of  an  unknown  material 

used  singly  or  as  the  inner  layer  of  a  three-layer  sandwich  structure. 

This  method  utilizes  a  theoretical  model  of  layered  structures  developed 
1-2  9-10 

by  Brekhovskikh  ’  which  incorporates  a  functional  dependence  upon 
the  density,  thickness,  and  longitudinal  sound  velocity  as  well  as  the 
quadratic,  linear,  and  constant  frequency-dependent  attenuation  coeffi¬ 
cients  for  each  layer  of  the  structure.  A  computer  fit  of  the  theore¬ 
tical  model  of  the  experimental  transmission-coefficient  spectrum  prov' 
a  resultant  set  of  acoustic  parameters  which  may  be  used  to  indicate 
material  degradation  such  as  microcrazing  or  to  evaluate  proposed  mode1'', 
of  tensile  or  shear-bond  strength.  The  cw  spectroscopy  results  are 
included  as  Appendix  A  in  a  paper  entitled,  "Computer  Determination 
of  Acoustic  Parameters  for  Multilayered  Structure  Using  cw  Spectroscopy". 

Table  I  in  Appendix  A  gives  an  overall  summary  of  the  procedure  used  to 
extract  the  acoustic  parameters  and  Table  II  shows  the  data  flow  for  the 
computer  program  for  the  layered  structures.  Appendix  B  contains  the 
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listing  of  the  fit  program  used  to  extract  the  physical  parameters  of  the 
inner  layer. 
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SECTION  III 


MEASUREMENT  OF  CRACKS  USING  THE  DOPPLER  EFFECT 

There  has  been  considerable  interest  in  the  nondestructive  detection  and 
measurement  of  radial  cracks  under  bolt  holes  of  wing  fasteners.  For 
detection  of  these  cracks,  conventional  pulse-echo  ultrasonic  techniques 
are  usually  employed.  In  general,  methods  for  measurement  of  the  size 
of  cracks  have  depended  upon  echo  duration  from  the  crack  or  upon  the 
attenuation  of  the  ultrasound  wave  across  the  crack.  It  will  be  demon¬ 
strated  in  this  section  that  for  special  geometries  such  as  the  radial 
cracks  emanating  from  the  bolt  hole  at  the  faying  surface  on  the  inner 
layer  of  an  adhesively  bonded  structure,  crack  size  can  be  measured 
using  the  Doppler  effect.  A  sample  having  a  radial  crack  is  mounted 
on  a  rotating  platform  in  a  water  bath  and  Interrogated  using  a  train 
of  tone  bursts  from  a  stationary,  fo9used  transducer.  The  Doppler-shif t 
spectrum  of  the  backscattered  ultrasonic  beam  from  the  rotating  crack  is 
recorded.  The  maximum  Doppler  shift  is  then  related  to  the  radial  size 
of  the  crack.  Paper  describing  the  Doppler  studies  of  cracks  is 
included  as  Appendix  C. 
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SECTION  IV 

PSEUDORANDOM-SIGNAL-CORRELATION  SYSTEM 

Detailed  information  on  a  pseudorandom-slgnal-correlatlon  system  designed 
and  built  under  this  program  (Including  a  logic  description) ,  on  the 
experimental  efforts  with  the  system,  and  on  the  experimental  results 
can  be  found  in  the  papers  which  appear  as  Appendix  D  of  this  report. 

This  section  contains  a  brief  discussion  of  a  cw  pseudorandom-signal- 
correlation  system  for  which  the  preliminary  design  is  now  complete. 

It  is  expected  that  future  work  on  this  system  will  give  new  insight 
into  ultrasonic  pseudorandom-signal  correlation. 

CW  PSEUDORANDOM-SIGNAL-CORRELATION  SYSTEM 
(DISCUSSION  OF  A  PRELIMINARY  DESIGN) 

The  pseudorandom-signal-correlation  system  discussed  in  detail  in  Appendix  D 
may  be  characterized  as  a  "burst-type"  system.  Since  the  bursts  are  sub¬ 
sequences  of  a  longer  sequence,  they  do  not  have  the  desirable  correlation 
properties  of  a  continuously  transmitted  complete  m-sequence.  Correlation 
of  these  truncated  m-sequences  produces  range  sidelobes.  In  the  burst-type 
system,  the  magnitude  of  the  range  sidelobes  is  reduced  by  an  averaging 
technique  which  consists  of  transmitting  many  bursts  per  correlation  point. 
The  increased  signal-acquisition  time  inherent  in  this  averaging  technique 
is  the  price  which  must  be  paid  for  the  advantages  of  a  burst-type  system. 
These  advantages  Include  the  capability  of  using  a  single  transducer  for 
both  transmitting  and  receiving  and  the  fine  control  over  burst  length 
which  is  Inherent  in  the  burst-generation  method  used. 

In  applications  where  two- transducer  operation  is  feasible,  a  pseudorandom- 
signal-correlation  system  which  continuously  transmits  a  complete 
m-sequence — a  "cw"  system — offers  certain  advantages  over  the  burst-type 
system.  The  auto-correlation  of  a  complete  N-blt  m-sequence  produces 
a  periodic  waveform  which  consists  of  triangular  pulses  spaced  N  bits 
apart  and  a  non-zero  constant  level  elsewhere.  For  a  sequence  having 
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the  binary  values  of  +1  and  -1,  the  pulses  will  be  N  units  high  and  the 
constant  level  will  be  -1  units.  (The  pulse  will  be  two  bits  wide  at  its 
base,  regardless  of  sequence  length.)  These  auto-correlation  properties 
imply  that  range  sidelobes  are  not  produced  in  a  cw  system  and  that  for 
any  sequence  length  which  does  not  produce  range  ambiguity,  the  cw  system 
can  transmit  more  energy  to  the  target  than  the  burst-type  system.  (Range 
ambiguity  results  when  all  discernible  signals  are  not  received  before  the 
m-sequence  begins  to  repeat  itself.) 

The  preliminary  design  of  the  cw  pseudorandom-signal-correlation  system  was 
completed  during  this  program.  This  system  is  intended  for  use  in  experi¬ 
mental  investigations  to  determine  how  well  a  practical  system  of  this  type 
would  conform  to  the  theoretical  ideal  performance  and  how  such  a  system 
could  be  applied  to  actual  NDE  problems. 

The  cw  pseudorandom-signal-correlation  system  will  be  a  flexible  laboratory 
instrument  incorporating  the  following  major  features: 

(1)  Capable  of  correlation  and  auto-correlation  using  N-bit 

N 

m-sequences  (of  length  2  -  1,  N  being  a  positive  non¬ 

zero  integer)  up  to  N  =  20. 

(2)  Possessing  two  basic  correlation  modes: 

(a)  "True"  continuous  correlation  where  the  correlation 
process  continues  until  a  halt  command  is  given. 

(b)  Delayed  correlation  over  a  preset  number  of  correla¬ 
tion  points.  In  this  case,  transmission  starts  on 
command,  but  the  reference  sequence  is  gated  off 
for  a  preset  delay  interval.  At  the  end  of  this 
delay,  both  the  reference  and  transmitter  sequence 
generators  run  for  a  preset  number  of  correlation 
points  after  which  the  entire  cycle  is  repeated. 
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In  either  of  the  above  correlation  modes,  a  correlation  point 
Is  defined  as  R  periods  of  the  m-sequence.  R  Is  a  presettable 
Integer  and  Is  analogous  to  "repetitions  per  correlation  point" 
In  the  burst-type  system.  The  repetitions  are  not  used  in  the 
cw  system  to  reduce  sldelobes  but  to  control  the  correlation 
"speed."  Just  as  in  the  burst-type  system,  correlation  speed 
determines  the  bandwidth  of  the  correlator  output. 

The  cw  system  algorithmic  design  is  essentially  complete.  Since  this  Is 
the  crucial  part  of  the  design,  any  anticipated  changes  (except  for  very 
minor  ones)  in  system  operation  should  be  made  before  proceeding.  The 
algorithm  is  now  being  reviewed  with  this  in  mind.  Future  work  will 
involve  completing  the  review  and  translating  the  algorithmic-state 
machine  chart  into  logic  diagrams  and  then  into  hardware. 


SECTION  V 

CONCLUSIONS  AND  FUTURE  WORK 


CONCLUSIONS 

cw  Spectroscopy 

A  swep t- frequency ,  continuous-wave  ultrasound  technique  for  determining 
longitudinal-mode  acoustic  parameters  in  multi-layered  structures  has 
been  developed.  A  computer  program  performs  a  nonlinear  curve  fit  of 
theroretical  to  experimental  frequency  spectra.  A  theoretical  model  of 
layered  structures  derived  by  Brekhovskikti  and  programmed  in  Fortran  is 
used  for  this  purpose.  A  computer  FIT  program  iteratively  computes  the 
wave  speed  and  frequency-dependent  attenuation  within  a  given  layer. 
Although  all  of  the  outer-layer-parameter  values  must  be  known,  charac¬ 
terization  of  the  inner  layer  requires  knowledge  only  of  the  thickness 
and  density.  ' 

Three-layer  aluminum/sealant /aluminum  sandwiches,  a  one- layer  aluminum 
plate,  and  several  one-layer  epoxy-resin  plates  have  been  characterized 
using  this  technique.  The  resultant  theoretical  frequency  spectra 
nearly  duplicated  the  experimentally  determined  spectra  in  most  cases. 
Figures  3-5  in  Appendix  A  (dashed  curves)  illustrate  the  FIT  program 
output  for  one  of  the  three- layer  samples. 

Comparison  of  the  epoxy-resin  results  with  independent  tone-burst  atten.m- 
tion  results  suggests  that  the  computed  values  for  attenuation  are  wit!  in 
±0.5  nepers/cm  of  the  actual  values  at  10  MHz  (broadband  10-MHz  trans¬ 
ducers  were  employed).  Two  different  thicknesses  of  a  standard  commercial 
sealant  (RTV)  used  as  the  inner  layer  of  three-layer  aluminum  plate 
sandwiches  were  evaluated.  The  computed  values  for  attenuation  per 
centimeter  varied  by  only  0.3  nepers. 
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Computed  wave  speeds  were  within  32  of  those  measured  for  the  epoxy-resin 
samples  by  the  tone  burst  technique.  The  computed  wave  speeds  for  the 
sealant  used  were  up  to  9 X  higher  than  the  measured  values.  This  dis¬ 
crepancy  is  probably  due  to  error  in  measuring  the  thickness  of  the 
inner  sealant  layer. 

A  major  result  of  cw  spectroscopy  was  the  development  of  a  versatile 
computer  program  for  determining  acoustic  parameters  for  a  variety  of 
layered  structures.  This  program  permits  the  calculation  of  the  ultra¬ 
sonic  attenuation  of  a  given  layer  in  a  three-layer  sandwich  structure 
directly  from  the  experimental  frequency  spectrum  of  that  structure. 

It  is  hoped  that  this  unique  ability  will  aid  in  the  detection  of 
sealant  and  adhesive  bond  degradation  in  layered  structures. 

The  following  conclusions  have  been  reached  as  a  result  of  this  study: 

•  The  knowledge  of  either  the  wave  speed  or  thickness  of 
inner  layer  demands  some  a  priori  knowledge  of  the 
properties  of  this  layer.  Hence,  this  technique  cannot 
be  applied  to  an  unknown  sealant. 

•  The  attenuation  as  determined  by  the  above  technique 
for  the  single-layer  epoxy-resin  samples  increased 
with  aging  of  the  samples.  This  was  confirmed  inde¬ 
pendently  by  the  tone-burst  method.  It  could  be 
inferred  from  this  that  the  periodic  determination 
of  the  attenuation  for  the  inner  layer  of  a  multi¬ 
layer  structure  may  give  an  indication  of  the 
degradation  or  the  mechanical  Integrity  of  the 
bond. 

•  It  proved  to  be  impractical  to  utilize  real  Air  Force 
sealants  in  the  present  study  due  to  their  high  attenua¬ 
tion.  A  commercial,  silicon-based  sealant  (RTV)  was 
substituted  because  of  its  ease  of  use  and  much  lower 
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attenuation.  Since  the  aluminum/sealant/alumlnum  sandwiches 
were  immersed  in  a  water  bath,  the  aluminum-water  impedance 
mismatch  created  another  source  of  signal  loss.  The  use  of 
impedance-matching  techniques — particular  when  real  sealants 
are  to  be  characterized — could  prove  to  be  necessary  in 
overcoming  this  problem. 

•  The  choice  of  the  Simplex  Method^  to  perform  the  nonlinear 
curve  fit  proved  to  be  very  advantageous.  The  method  is 
ideally  suited  for  the  complex  and  (nearly)  intractable 
Brekovskikh  Model.  Mo  additional  subroutines  are  required 
to  invert  matrices  or  to  determine  derivations.  As  cur¬ 
rently  implemented,  the  program  can  handle  up  to  100  param¬ 
eters  with  the  user  having  freedom  to  vary  or  fix  those  as 
desired.  This  flexibility  was  invaluable  in  evaluating  the 
technique. 

Measurement  of  Crack  Size 

As  shown  in  Section  III,  it  is  Indeed  possible  to  measure  the  radial  size 
of  the  cracks  under  bolt  holes  both  in  single-  and  multi-layer  mechanically 
fastened  aircraft  wing  structures.  The  accuracy  of  the  present  method  is 
better  than  20%.  In  the  presence  of  multiple  radial  cracks  or  nonradial 
cracks,  the  method  can  measure  only  the  maximum  radial  extent.  However, 
this  is  sufficient  since  with  knowledge  of  the  radial  extension  of  the 
crack,  the  bolt  hole  can  be  reamed  with  an  over-size  reamer  and  a  crack 
can  be  repaired  before  it  reaches  a  critical  size. 

Since  the  size  of  the  radial  crack  is  computed  from  the  frequency  shift  and 
not  from  the  amplitude  of  the  backscattered  ultrasound,  the  method  is 
applicable  for  measuring  radial  cracks  in  the  third  layer  in  a  multiple- 
layer  aircraft  structure.  The  method  can  be  used  in  the  field  by  employing 
a  rotating  transducer.  The  sound  is  coupled  to  the  teat  sample  through 
water  in  a  boot  capping  the  transducer.  Results  were  similar  to  those 
obtained  by  rotating  the  disk  in  water. 


-  L  ^ 
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Ultrasonic  Pseudorandom-Signal-Correlation  System 


The  experimental  characterization  of  the  Ultrasonic  Pseudorandom-Signal- 
Correlation  System  (UPSC  system)  has  shown  that  such  a  system  can  produce 
results  comparable  to  those  obtained  with  a  random  signal-correlation 
system.  The  UPSC  system  has  better  range  resolution  than  a  random-signal 
system  which  makes  use  of  a  narrow-band  variable  delay  line.  However,  a 
high  system  clock  rate  is  required  to  generate  and  control  the  wide  band¬ 
width  pseudo-noise  bursts  needed  for  good  range  resolution.  Furthermore, 
the  signal-to-noise  ratio  improvement  gained  with  the  UPSC  system  is 
obtained  at  the  expense  of  long  signal-acquisition  times. 

FUTURE  WORK 

Research  should  be  initiated  on  the  special  case  of  a  three-layer  structure 
where  two  matching  pieces  of  metal  are  separated  by  a  thin  liquid  layer. 

This  type  of  structure  may  improve  efficiency  in  nondestructive  inspection 
of  cast,  "near-net-shape”  components  such  as  turbine  disks  which  contain 
topologically  complex  features  which  are  difficult  to  scan  in  a  typical 
water- tank  inspection  device.  It  is  hoped  that  a  matching  part  to  the 
turbine  disk,  for  instance,  would  limit  the  diffraction  and  focusing  effects 
of  the  water-metal  interface  by  decreasing  the  path  length  of  ultrasound  in 
water.  In  the  final  limit  of  an  infinitesimally  thin  water  layer,  the  com¬ 
ponent  to  be  inspected  should  be  seen  only  as  an  interior  region  of  a 
larger,  but  topologically  simpler,  structure. 

Initially  a  comparison  should  be  made  between  theoretically  predicted  losses 
due  to  thin  water  interfaces  at  various  angles  and  thicknesses  and  losses 
observed  experimentally.  Figure  1  illustrates  the  transmission  and  reflec¬ 
tion  coefficients  of  a  5-mil-thick  aluminum/water/alumlnum  interface. 

A  proof-of-principal  device  to  test  the  matching-shape  hypothesis  has  been 
constructed  and  used  with  encouraging  results.  Figure  2  illustrates  this 
device  which  consists  of  a  cylindrical  well  within  an  aluminum  cube.  An 
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aluminum  rod  may  be  Inserted  Into  the  well  with  a  glycerin-filled,  10-mil 
clearance  between  the  rod  and  well  wall.  An  off-center,  1/8-in.  hole  was 
drilled  into  the  rod  to  provide  a  simulated  defect.  By  scanning  the  outer 
surfaces  of  the  aluminum  cube,  this  discontinuity  could  be  detected;  the 
ultrasonic  return  echo  was  noted  to  shift  in  time  as  the  rod  was  rotated 
in  the  well. 


DISCONTINUITY 


Fig.  2.  Experimental  Setup 


Another  future  effort  which  should  be  undertaken  is  automation  of  the  mea¬ 
surement  of  the  radial  size  crack.  Radial  bolt-hole  cracks  can  be  detected 
TM 

by  the  Autoscan  (I,  II,  and  III),  an  instrument  which  can  center  the 
ultrasonic  beam  on  the  bolt  hole  rather  conveniently.  A  tone  burst  rather 
than  pulses  can  be  used  as  the  incident  beam.  The  Doppler  spectra  from 
the  crack  could  be  stored  or  processed  by  a  microprocessor  to  measure  the 
maximum  Doppler  frequency.  This  frequency,  fQ,  then  would  be  fed  into  the 
alogrlthm  for  calculating  the  crack  size,  which  could  be  displayed  in 
real  time.  However,  such  a  study  may  fall  within  the  domain  of 
Manufacturing-Technology  programs. 
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A  high  clock  rate  for  the  UPSC  system  could  be  obtained  by  using  emitter- 
coupled  logic  (ECL) ,  while  the  signal  acquisition  time  could  be  improved 
by  using  parallel  processing  with  multiple  correlators.  These  two  system 
improvements  might  be  the  subject  of  future  work  on  the  UPSC  system. 

Future  work  on  pseudorandom  signal-correlation  systems  should  also  include 
the  completion  of  the  design*  fabrication,  and  experimental  evaluation 
of  the  cw  pseudorandom  signal  correlation  system. 
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APPENDIX  A 

CW  SPECTROSCOPY  STUDIES 
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COMPUTER  DETERMINATION  OF  ACOUSTIC  PARAMETERS  FOR 
MULTILAYERED  STRUCTURES  USING  CW  ULTRASONIC  SPECTROSCOPY 


J.  S.  Leffler  and  N.  K.  Batra 
Systems  Research  Laboratories,  Inc. 
Dayton,  OH 


INTRODUCTION 

The  Increasing  role  of  multilayered,  adhesively  bonded  aircraft  structures 
requires  a  corresponding  advance  in  nondestructive  evaluation  techniques.  In  par¬ 
ticular,  the  Inaccessibility  of  the  Inner  layers  of  such  sandwlch-llke  structures 
renders  an  Isolated,  layer-by-layer  examination  of  the  structure  Impossible.  Thus, 
any  ultrasonic  examination  of  the  inner  layers  of  the  component  In  question  must  take 
into  account  Interference  due  to  the  outer  layers.  Characterization  of  inner-layer 
physical  properties  can  be  accomplished  by  means  of  ultrasonic  spectroscopy  which. 

In  effect.  Isolates  the  contribution  of  each  material  layer  to  the  total  frequency 
spectrum  of  the  multilayered  structure. 

This  paper  concerns  the  utilization  of  ultrasonic  spectroscopy  to  determine  the 
acoustic  parameters  of  an  "unknown"  material  used  singly  or  as  the  Inner  layer  of  a 
three-layer  sandwich  structure.  The  resultant  set  of  parameters  may  be  used  to  Indi¬ 
cate  material  degradation  such  as  microcrazing  or  to  evaluate  proposed  models  of  ten¬ 
sile  or  shear  bond  strength. 


THEORETICAL  CALCULATIONS 

Determination  of  the  acoustic  parameters  la  accomplished  by  means  of  a  computer 
program  which  performs  a  nonlinear  curve  fit  of  a  theoretical  frequency  spectrum  to 
the  experimental  frequency  spectrum  for  the  structure  In  question.  The  theoretical 
model  used  was  derived  by  Brekhovskikh*.  A  subroutine  for  implementing  this  model 
calculates  the  transmission-coefficient  frequency  spectrum  for  normal-incidence 
ultrasound  given  the  density,  thickness,  and  longitudinal  sound  velocity  as  well 
as  the  quadratic,  linear,  and  constant  frequency-dependent  attenuation  coefficients 
of  each  layer  of  the  structure. 

The  nonlinear  curve-fit  program  varies  the  unknown  inner-layer  parameters  in 
such  a  way  as  to  minimize  the  differences  between  the  theoretical  and  experimental 
frequency  spectra.  This  process  Is  Iterative  in  nature.  All  outer-layer  acoustic 
parameters  and  at  least  one  inner-layer  acoustic  parameter  (thickness  or  wavespeed) 
are  "fixed”  (i.e.,  set  equal  to  some  known  constant).  The  remaining  inner-layer 
parameters  are  initially  estimated.  A  theoretical  frequency  spectrum  is  then 
generated  by  the  computer  model,  based  upon  the  initial  parameter  set,  and  compared 
to  the  original  experimental  frequency  spectrum.  This  comparison  yields  a  modified 
chl-squared  value  which  is  related  to  the  "closeness-of-f it"  of  the  two  spectra. 

The  nonlinear  fit  program  makes  use  of  this  chi-squared  value  as  a  guide  In  generat¬ 
ing  a  new  parameter  set  which  will  better  fit  the  experimental  results.  The  computer 
model  utilizes  the  new  parameter  set  In  calculating  another  theoretical  frequency 
spectrum,  and  the  iteration  continues  until  the  best  chi-squared  fit  Is  achieved. 

The  acoustic  parameters  for  the  inner  layer  of  the  structure  are  then  given  by  the 
last  Iterated  parameter  set. 

-  ■ 

Work  supported  under  U.S.  Air  Force  Contract  No;  F33615-77-C-5022. 
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EXPERIMENTAL  PROCEDURE 


The  experimental  set-up  Is  Illustrated  In  Figure  1.  A  digital  synthesizer  Is 
swept  In  frequency  (0-20  MHz)  by  a  1024-channel  signal  averager.  The  output  sine 
wave  Is  fed  to  an  Immersion  transducer  and  to  the  reference  Input  of  a  gain/phase 
meter.  The  Incident  ultrasound  is  coupled  through  water  to  the  sample  at  normal 
Incidence  setting  up  standing-wave  resonances  at  certain  frequencies.  The  ultra¬ 
sound  transmitted  through  the  sample  Is  received  by  a  second  transducer  and  com¬ 
pared  to  the  reference  signal  by  the  gain/phase  meter.  An  analog  output  which  can 
selectively  represent  either  the  relative  attenuation  (In  dB)  or  the  phase  change 
(In  degrees)  of  the  transmitted  versus  the  Incident  signal  Is  digitized  by  the  signal 
averager  and  stored  on  magnetic  tape. 

Figure  2  shows  a  three-layer  sample  spectra  which  has  been  deconvolved  to  elimi¬ 
nate  the  effects  of  the  particular  transducers  or  equipment  used.  The  peaks  corre¬ 
spond  to  standing-wave  resonances  Induced  in  the  three  material  layers. 

RESULTS 

Three-layer  structures  (alumlnum/sealant/alumlnum)  (see  Figures  3-5)  and  several 
one-layer  epoxy-resin  plates  (Figures  6-7)  have  been  characterized  using  the  computer 
program.  Each  figure  contains  a  table  giving  the  complete  acoustic  parameter  set  for 
that  sample.  T,  p,  and  C  represent  thickness  (cm),  density  (gm/cc),  and  longitudinal 
wave  velocity,  respectively.  The  total  frequency-dependent  longitudinal  attentuatlon 
Is  given  by 


a(u>)  “  C2'w2  +  +  C0  , 

where  C2,  Cl,  and  C0  are  the  quadratic,  linear,  and  constant  attenuation  coefficients, 
respectively.  (The  sheAr-mode  parameters,  B,  B2,  Bl,  and  B0  In  the  figures  were  not 
used  In  this  longitudinal  wave  study  and  were  set  equal  to  zero.)  Since  the  flrst- 
and  third-layer  aluminum  plates  were  identical,  parameters  are  specified  for  the 
first  layer  only.  In  the  lower  section  of  each  figure,  FMIN  and  FMAX  denote  the 
frequency  range  of  interest  and  WN(2)  is  a  scale  factor  normally  set  equal  to  one. 
Theoretical  and  experimental  spectra  are  signified  by  dashed  and  solid  lines,  respec¬ 
tively. 


Three-Layer  Structures 

Figure  3  shows  the  theoretical  frequency  spectrum  derived  from  the  initial  parame¬ 
ter  estimation  for  a  particular  three-layer  structure.  The  theoretical  spectrum 
based  upon  computer-derived  acoustic  parameters  for  that  same  structure  is  shown  in 
Figure  4.  The  value  of  chi-squared  goes  from  16.1  in  Figure  3  to  2.21  in  Figure  4, 
which  reflects  a  much  better  agreement  between  the  theoretical  and  experimental  spec¬ 
tra  in  the  latter  plot.  Figure  5  shows  computer-determined  acoustic  parameters  for 
a  similar  sandwich  structure  but  one  having  a  "double-thick"  inner  layer  (sealant). 

The  constant  attenuation  coefficient,  C0,  computed  for  the  sealant  (same  type  in 
both  samples)  was  2.03  nepers/cm  in  Figure  5  and  1.89  nepers/cm  in  Figure  4.  This 
close  agreement  tends  to  validate  the  procedure  for  determining  inner-layer  parame¬ 
ters. 


Single-Layer  Structures 

Single-layer  epoxy-resin  samples  which  had  been  hydrothermally  spiked  to  create 
scattering  centers  were  studied  to  determine  the  effect  of  sample  degradation  upon 
the  attenuation  of  ultrasound.  Figures  6  and  7  show  the  computer-determined 
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attenuation  parameter*  (C2,  Cl,  and  C0)  for  both  unaged  and  aged  (hydrothermally , 

318  daya)  caaea.  Notice  the  quadratic  frequency  dependence  of  the  attenuation  in 
both  casea.  A  comparison  of  typical  nonlinear-f lt-program  reaulta  with  those  from 
an  Independent  analysis  using  a  tone-burst  method  Is  shown  In  Table  1.  It  ahould 
be  pointed  out  that  the  attenuation  Increases  with  aging. 

SUMMARY 

Continuous-wave  ultrasonic  spectroscopy  Indeed  provides  a  tool  for  the  charac¬ 
terization  of  the  inner  layer  of  multilayered  structures.  Through  the  use  of  a  com¬ 
puter  curve-fit  technique  and  a  suitable  theoretical  model,  acoustic  parameters  can 
be  determined  from  analysis  of  experimental  frequency  spectra.  The  determined  parame¬ 
ters  can  be  used  to  characterize  the  layered  structure. 
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TAPE  TRANSPORT 


ANALOG  A/D  CONVERTER 
INPUT  I  AND  STORAGE 


Q=P 


SWEEP  CLOCK  JUTTL 


GAIN/PHASE 

METER  REFERENCE  SIGNAL 


DIGITAL 

SYNTHESIZER 


RECEIVER 

INPUT _ e 

RECEIVER  ^ 
TRANSDUCER 


RECEIVER  PKTRANSMITTER 
TRANSDUCER  ^TRANSDUCER 
W  ^ADHESIVE 
BOND 
SAMPLE 

FIGURE  1 — CW  ultrasonic  spectroscopy  experimental  arrangement. 


FIGURE  2 — Deconvolved  experimental  frequency  spectrum  for  three-layer  aluminum/ 
•ealant /aluminum  structure. 


t 


FIGURE 

spiked" 


7— Computer-determined  acoustic  parameters  for. single^ layer  "hydrothermally 
epoxy  resin. 
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PROGRAM  FIT (INPUT, OUTPUT, TAPE6,DEBUG=0UT PUT, PLFILE*0) 
COMHON/FITl/NPTS*NPAR,YD(105l)  ,  WOdOEu)  ,  PZ  (IOC)  ,  TO  (1050) 
COMMON  /FITB/  FM IN , FM AX, I T ITLE  (5) 

COMMON  /FITC/  IPO ( 10  3 ) , I°SMCM 
COMMON  /FITO/  9(2) ,02(2) , 31(2) , 30<E) 

COMMON/ST 2PIT/NV ,NTRACE,CHISQ»MASK(100) , 

1  X (100), DEL  TAX (10C) , DEL  MIN (10C) , NF 

COMMON  /1EST/  9X ( 100 ) ,NF9£ST , NRP3ST , 8CHI5D 
COMMON  /CONST5/  NPLN, NPlSH, NPINT.NEVAL 
EXTERNAL  FCN 
IPSMCH  r  1 

PLOT  EXPERIMENTAL  curve? 

NPLWS  =  0  ;  PLOT  THEORETICAL  L  EXPERIMENTAL 
NPLSW  =  2  5  PLOT  ONLY  THEORETICAL 
NPLSH  =  0 

LIMIT  *  OF  EVALUATIONS 
NFVAL=99 

PLOT  EVERY  50  EVALUATIONS 
NPINT  =  50 

SET  UP  FOR  BEST  FIT 
9CHISQ  =  1.0E20 
NF  «  0 

NORMAL  INCIDENCE  ONLY 
DO  7  I  =  1,2 
9(1)  =  0.0 
92(1)  -  0.0 
91(1)  =0.0 
90(1)  =  0.0 
7  CONTINUE 


THIS  PROGRAM  WILL  ACCEPT  EXPERIMENTAL  OR  THEORETICAL  DATA  IN 
FILE  T APEb  ALONG  WITH  KNOWN  AND  "GUESS TI MATEO"  PARAMETERS  FOR 
A  THREE  LAYER  METAL-90N0-.NET AL  SYSTEM  AND  ATTEMPT  TO  FIT  THE 
UNKNOWN  PARAMETERS  TO  THE  OATA. 


FOR  THE  THREE  LAYER  SYSTEM,  12  PARAMETERS  ARE  INPUTTED  AND 
STORED  IN  THE  PZ  ARRAY  BY  THE  INOATA  SUBROUTINE  AS  FOLLOWS) 

THICKNESS  OF  FIRST  LAYER  (CM). 

THICKNESS  OF  SECOND  (MIDDLE)  LAYER  (CM). 
VELOCITY  OF  SOUND  IN  FIRST  LAVER  (CM/SEC). 
VELOCITY  OF  SOUND  IN  SECONO  LAYER  (CM/SEC). 
RELATIVE  OFNSITY  OF  FIRST  LAVER. 

RELATIVE  DENSITY  OF  SECOND  LAYER. 

OMEGA  SQUARED  LONGITUDINAL  ATTENUATION 


PZ( 

1) 

r 

T  ( 1 ) 

PZ( 

2) 

= 

T  (2  ) 

PZ( 

3) 

r 

C  ( 1 ) 

PZ( 

4) 

s 

C  ( 2) 

PZ( 

5) 

s 

R(l> 

PZ( 

6) 

r 

R(  2 ) 

PZ( 

7) 

r 

A2  ( 1  > 

rUlS  FASX  IS  H’JAIiXYT 

.  .  70  DO  3 
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COEFFICIENT  FOR  FIRST  LAYER. 

PZ(  8)  =A2(21  -  OHEGA  SQUARED  LONGITUDINAL  ATTENUATION 
COEFFICIENT  FOR  SECOND  LAYER. 

PZ<  91  sAl(l)  -  OMEGA  LONGITUDINAL  ATTENUATION  COEFFICIENT 
FOR  FIRST  LAYER. 

PZ(IJ)  = A1 ( 2)  -  OMEGA  LONGITUDINAL  ATTENUATION  COEFFICIENT 
FOR  SECONO  LAYER. 

PZ (111  = AO  ( 1 1  -  CONSTANT  LONGITUDINAL  ATTENUATION  TERM  FOR 
FIRST  LAYER. 

PZ(121  = AO (1)  -  CONSTANT  LONGITUOINAl  ATTENUATION  TERM  FOR 
SECONO  LAYER. 

*  THE  THIRD  LAYER  IS  ASSUMED  TO  BE  IDENTICAL  TO  THE  FIRST. 


MASK(I1  INDICATES  WHETHER  OR  NOT  THE  I'TH  PARAMETER  IS  KNOWN. 
IF  MASK ( I 1 = I.  THE  I'TH  PARAMETER  IS  DETERMINED  AND  MAY  NOT  9E 
VARIED  IN  ORDER  TO  03TAIN  A  FIT.  IF  MASK(I)=0,  THE  I'TH 
PARAMETER  IS  UNKNOWN  ANO  REPRESENTS  A  "GUESSTIMATE”  WHICH  MAY 
BE  VARIED  WHILE  OBTAINING  A  FIT  TO  THE  DATA. 


THE  YD  ARRAY  CONTAINS  THE  REAL  PART  OF  THE  COMPLEX 
TRANSMISSION  ^UNCTION  ARRAY  READ  IN  FROM  FILE  TAPE6  IN 
SUBROUTINE  INDATA.  THIS  FILE  IS  PRODUCED  FROM  EXPERIMENTAL  DATA 
BY  THE  OATA  PROGRAM  OR  FROM  A  THEORETICAL  MODEL  3Y  THE  LAYER 
PROGRAM.  EACH  CHOICE  OF  PARAMETERS  LEADS  TO  A  THEORETICAL 
COMPLEX  TRANSMISSION  FUNCTION  WHICH  IS  COMPAREO  TO  YO  IN  OROER 
TO  SURMISE  THE  CLOSENESS  OF  THE  FIT.  THIS  IS  OONE  BY  THE 
FCN  SUBROUTINE. 


THE  WD  ARRAY  IS  USED  AS  A  WEIGHTING  FUNCTION  TO  FIT  THE  DATA 
WITH  REGARD  TO  THE  TRANSDUCER  RESPONSE.  MORE  IMPORTANCE  IS 
GIVEN  TO  HATCHING  THE  THEORETICALLY  GENERATED  TRANSMISSION 
FUNCTION  WITH  THE  EXPERIMENTAL  TRANSMISSION  FUNCTION  AT  THE 
FREQUENCIES  WHERE  THE  TRANSDUCER  RESPONSE  IS  THE  GREATEST. 


THE  SEARCH  SUBROUTINE  PERFORMS  A  NON-LINEAR  FUNCTION 
MINIMIZATION  3Y  THE  USE  OF  A  SIMPLEX  METHOD.  IT  VARIES  THE 
UNKNOWN  PARAMETERS  AND  CALLS  SUBROUTINE  FCN  TO  DETERMINE 
THE  SUCCESS  OF  THE  LATEST  ITERATION  BY  USE  OF  A  CHI-SQUAREO 
TEST. 


VARIABLES 


NAME  TYPE  FORMAT  DEFINITION 


Ml»  TM2  IS  BtST  QUALITY  MUfctt**** 

PROM  99?  T  W3 
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AO 

REAL  ARRAY 

E10.2 

CONSTANT  ATT.  TERM  IN  LAYER 

A1 

REAL  ARRAY 

E10.2 

• 

W  ATT.  COEFFICIENT  IN  LAYER 

A2 

REAL  ARRAY 

E10.2 

- 

W**2  ATT.  COEFFICIENT  IN  LAYER 

3CHISQ 

real 

m 

BEST  FIT  CHI-SQUAREO  VALUE 

BX 

REAL  ARRAY 

(130) 

- 

BEST  FIT  PARAMETER  VALUES 

C 

REAL  ARRAY 

E10.2 

- 

VELOCITY  OF  SOUNO  IN  LAYER  (CM/S) 

CHISQ 

REAL 

- 

CHI-SQUAREO 

0 

COMPLEX  ARRAY 

£15.5 

- 

TEMPORARY  STORAGE, TRANS.  FUNCTION 

OELMIN 

REAL  ARRAr 

(100) 

• 

SMALLEST  STEP  SIZE  FOR  X(I) 

OELTAF 

real 

FREE 

- 

FREQUENCY  INCREMENT  FOR  DATA  SET 

OELTAX 

real  ARRAY 

(103) 

• 

INITIAL  STEP  SIZE  FOR  X(I) 

OUM 

REAL  array 

- 

DUMMY  FILLER  FOR  COMMON  STATEMENT 

0UM1-8 

real  araays 

- 

DUMMY  FILLER  FOR  COMMON  STATEMENT 

EYE 

complex  const. 

• 

=SQRT (-1) 

FMAX 

REAL 

£10.2 

- 

MAXIMUM  FREQUENCY  FOR  VALID  FIT 

FMIN 

REAL 

£19.2 

- 

MINIMUM  FREQUENCY  FOR  VALIO  FIT 

FZERO 

REAL 

FREE 

- 

STARTING  FREQUENCY  FOR  DATA  SET 

I 

INTEGER 

- 

DO  LOOP  PARAMETER 

IPO 

INTEGER  ARRAY 

(100) 

- 

PEAK  (S.W.  RESONANCE)  INDICES 

IPSWCH 

INTEGER 

- 

IF  =1, UTILIZES  DATA  PEAKS  IN  FIT 

ITITLE 

INTEGER  ARRAY 

4A10 

- 

CONTAINS  PLOT  TITLE 

MASK 

INTEGER  ARRAY 

(100) 

• 

FIXED  VARIABLE  MASK 

NEVAL 

INTEGER 

- 

MAX  #  OF  EVALUATIONS 

NF 

INTEGER 

- 

CURRENT  NUMBER  OF  EVALUATIONS 

NFBEST 

NTEGER 

- 

BEST  FIT  *  OF  ITERATIONS 

NPAR 

INTEGER 

- 

#  OF  PARAMETERS 

NPAR1 

INTEGER 

« 

=NPAR-1,  DUMMY  VARIA9LE 

NPINT 

INTEGER 

• 

PLOT  INTERVAL  ( #EVAL 'S/PLOT) 

NPLN 

INTEGER 

• 

PLOT  NUMBER 

NPLSW 

INTEGER 

- 

=2, PLOT  ONLY  THEORETICAL, 0=9OTH 

NPTS 

INTEGER 

FREE 

- 

NUMBER  OF  POINTS  IN  DATA  SET 

NRP3ST 

INTEGER 

• 

#  REJECT  PEAKS(FCN) (NOT  YET  IMPL. 

NTRACE 

INTEGER 

- 

=0 , no  trace:  =i, trace  search 

NV 

INTEGER 

• 

f  OF  parameters  to  be  varied 

NAME 

type 

FORMAT 

PI 

real 

CONSTANT 

PZ 

real 

ARRAY 

(100) 

R 

real 

ARRAY 

£10.2 

T 

REAL 

ARRAY 

E  10.  2 

TO 

REAL 

ARRAY 

(1050) 

TI 

REAL 

F  R*-E 

WO 

REAL 

ARRAY 

(100) 

X 

X  JUMP 

REAL 

REAL 

ARRAY 

(100) 

YO 

REAL 

ARRAY 

(100) 

Z1 

real 

CONSTANT 

ZN1 

real 

CONSTANT 

DEFINITION 


=3.14159 

contains  initial  parameters 

RELATIVE  DENSITY (RHO)  OF  LAYER 
THICKNESS  OF  LAYER  (CM) 

THEOR.  DATA  POINTS  AT  EACH  FREQ. 
ANGLE  OF  INCIDENCE  (OEG.) 
WEIGHTING  FUNCTION  -  T.  RESPONSE 
CONTAINS  3EST-FIT  COMPUTED  PARAMS 
DUMMY  VARIABLE  =  107. 

REAL  (O)-RLAL  PART  OF  TRANS.  FNC , 
IMPEDANCE  OF  FIRST  LIQUID  MEDIA 
IMPEDANCE  OF  LAST  LIQUID  MEOIA 
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INPUT  OATA  CAROS 


CARO  #1 

-  FIRST  LAYER  (6E10.2) 


CC  1-10 
CC  11-20 
CC  21-30 
CC  31-40 
CC  41-50 
CC  51-60 


R  I  RELATIVE  DENSITY 
T  «  THICKNESS  (CM) 

C  t  VELOCITY  OF  SOUND  (CM/SEC) 

62  I  OMIGA-SQUAREO  LONGITUDINAL  6TTEN.  COL  FF ICIENT 
61  l  OMEGA  LONGITUDINAL  ATTENUATION  COEFFICIENT 
60  I  CONSTANT  ATTENUATION  COEFFICIENT 


CARD  »2 

-  SECOND  (MIODLE)  LAYER 

SAME  FORMAT  AS  FIRST  LAYER  PARAMETERS 

CARO  #3 

-  SECOND  LAYER  MASK  (8011) 

CC  1  t  =1  IF  T  FIXED, 0  (OR  BLANK)  IF  NOT  FIXEO. 

CC  2  I  =1  IF  C  FIXED, 0  (OR  BLANK)  IF  NOT  FIXED. 

CC  3  I  =1  IF  R  FIXEO, 0  (OR  BLANK)  IF  NOT  FIXEO. 

CC  4  »  =1  IF  A2  FIXED, 0  (OR  eLANK)  IF  NOT  FIXEO. 

CC  5  I  =1  IF  A 1  FIXFO, C  (OR  BLANK)  IF  NOT  FIXED. 

CC  6  I  =1  IF  AO  FIXED, 0  (OR  BLANK)  IF  NOT  FIXEO. 

CC  11-20  -  FMIN  :  MINIMUM  FREQUENCY  FOR  FIT ,OEFAULT=S.O 
CC  21-30  -  FMAX  I  MAXIMUM  FREQUENCY  FOR  FIT , DEFAULT =2 0. 0E06 
CC  31-70  -  ITITLEJ  PLOT  TITLE,  h 0  CHARACTERS 
CCSO  I  =1  IF  WEIGHTING  FUNCTION  NOT  NECESSARY 


INPUT  FILE  TAPE6  (  FROM  OATA  OR  LAYER  PROGRAM) 


4110  -  (TRN( I) , 1=1, 4) 5  EXPERIMENTAL  OATA  TAPE  RECORD  NUMBERS 
FREE  -  NPTS  ?  NUMBER  OF  DATA  POINTS  IN  OATA  SET 
FREE  -  FZERO  5  START  FREQUENCY  (HZ) 

FREE  -  OELTAFJ  FREQUENCY  INCREMFNT  (HZ) 

FREE  -  TI  5  ANGLE  OF  INCIDENCE  ( DEG, ) 

2E15.5  -  (0(1) ,1=1, NPTS) ?  COMPLEX  TRANSMISSION  FUNCTION 
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MAIN'  PROGRAM - * - 

CALL  PLOTTER  INITIALIZATION  ROUTINE 

CALL  COMPRS 

NPLN  =  0 

PRINT  4 

PRINT  6 

GET  OATA 

CALL  INOATA 

SET  UP  FOR  00  LOOPS 
NPAR1  =  NPAR  -  1 

LOCATE  OATA  PEAKS 
CALL  PEAK(Y0,NPTS,IP0,1) 

TOP  OF  FORM 
PRINT  4 
PRINT  5 

PRINT  OUT  INITIAL  PARAMETERS 
PRINT  11 

PRINT  13, (ITITLE(I) ,  1=1,4) 

PRINT  10,  (  1,PZ(I),  1=1, N PARI, 2) , (2,PZ(I) ,1=2, NPAR, 2) 

IF  (  MASK (2)  .EO.  1  )  PRINT  22 

IF  (  MASK (4)  .£Q.  1  )  PRINT  24 

IF  (  MAS K(6)  .EQ.  1  )  PRINT  26 

IF  (  MASK (8)  .EO.  1  )  PRINT  28 

IF  (MASK ( 10)  .EQ.  1)  °RINT  30 

IF  (MASK  (12)  .EQ.  1)  PRINT  32 

COMPUTE  STEP  SIZE  1  ~  '*  *-  '  . . 

00  1  1=1, NPAR 
XJUMP=0. 1 

OELTAX(I)=XJUMP*PZ(I) 

IF (OELTAX(I)  .EQ.  0.0)  OELTAX  (I) =0. 01 
HATCH  OUT  FOR  A2(2) 

IF  (  I  .EQ.  8  )  DELTAX(I) =1.0E-14 
DELMIN(I)=XJUMP*XJUMP*OELTAX(I) 

1  CONTINUE 

TOP  OF  FORM 
PRINT  4 

COMPUTE  CORRECT  PARAMETERS 
PLOT  FIRST  GUESS 

CALL  fcn  . . . 

CALL  PICTR(l) 

CALL  SMPLEX (FCN) 

PLOT  COMPUTED  FIT 
CALL  FCN 
CALL  PICTR(l) 

TOP  OF  FORM 
PRINT  4 
PRINT  5 


'■-a  “•  ' 
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PRINT  OUT  INITIAL-  PARAMETERS 
PRINT  11 

PRINT  13, (ITITLE (I) ,  1*1,4) 

PRINT  13,  (  1 , P  Z ( I » ,  I=1,NPAR1,2> , (2,PZ(I> ,I=2,NPAR,2) 

PRINT  5 

PRINT  OUT  BEST  COMPUTED  PARAMETERS 
PRINT  12 

PRINT  14,  3CHISR ,  NFSEST 

PRINT  13,  ( 1 , BX ( I ) ,  1  =  1 » NPAR1 ,2) , (2 , 9X ( I ) , I =2,NPAR,2) 

IF  (  MAS X ( 2)  .EQ.  1  )  PRINT  22 

IF  (  MASK ( 4 )  .EQ.  1  >  PRINT  24 

IF  (  MASK ( 6)  .EO.  1  )  PRINT  26 

IF  <  MASK ( 8)  .EQ.  1  >  PRINT  2« 

IF  (MAS<(10)  .EQ.  1)  PRINT  30 
IF  (MASK ( 12)  .EQ.  1)  PRINT  32 
PRINT  5 
PRINT  3, NF 
CALL  OOMEPL 

FORMAT  STATEMENTS 

2  FORMAT (10X,6E15.5) 

3  FORMAT <10X, ’NUMBER  OF  FUNCTIONAL  EVALUATIONS  =  *,I5) 

4  FORMA T ( 1H1) 

5  FORMAT (/////) 

6  FORMAT (*-• ,13X,  *FIT  PROGRAM  EXECUTION*) 

1C  FORMAT  <*(>♦, 10X,  *T(*,I1,*>  =  *,  E15 . 5 , 5X  ,  *C  (*  ,  II ,  *)  =  *,E15.5*5X, 

1  *RHO ( * , 1 1 , *)  =  * , £15, S ,/*  *,10X,*A2(*,I1,*)  =  *  ,C15.5,5X, 

2  *A1(*,I1,*)  =  *,E15.5,5X,*A0(*,I1,*>  =  *,E15.5,5X, 

3  *HN(*»I1»*)  =  *,F5.2) 

11  FORMAT  (*-  INITIAL  PARAMETERS  *) 

12  FORMAT  (*-  COMPUTSO  PARAMETERS*) 

13  FORMAT (*3*,13X,4A10,//> 

1L  FORMAT  <*0*,SX,*CHISi  =  *,F12.3,*,  EVALUATION  NUMBER  =  *,I5,//> 

22  FORMAT  (*0*,10X,  *FIX  THICKNESS*) 

24  FORMAT  (*3*,10X,  *FI<  VELOCITY*) 

26  FORMAT  ( *  0* , 10 X , *  FIX  OENSITY*) 

28  FORMAT  ( *  0* , 13  X , *FIX  A2*) 

30  FORMAT  (* 0* , It  X , *FIX  Al*) 

32  FORMAT  { *  3* , 10X , *  FIX  A0*) 

STOP 

ENO 

SUBROUTINE  SMPLEX (FUNK) 


FUNCTION  MINIMIZATION  BY  THE  SIMPLEX  METHOO  SMPLX 

(NELOER  AND  MEAD,  THE  COMPUTER  JOURNAL,  JANUARY  1965)  SMPLX  4 

COMMANOMENf  1  ....  THOU  SHALT  NOT  OMIT  THINE  "EXTERNAL"  CARO.  SMPLX  8 


*****  SEARCH  COMMON  •*♦**♦ 
COMMON/ST £PI T/NV , NTRACE, CH ISQ, MASK ( 100) » 


OOOOOOO  O  OOOOOOOOOOOOOOOOOOO  o  o  o 


1  X(100),OELTAX(10C) ,OELMIN(10P),NF 

COMMON  /C0NST6/  NPLN, NPLSW , NPINT ,LIM 

**•*•  INTERNAL  DIMENSIONS  ***** 

OIMENSION  CHI(2i> ,Z(21,2C) ,ZBAR<20> ,ZSTAR(20)  SHPLX  12 


THE  SEARCH  PARAMETERS  ARE  OEFINEO  AS  FOLLOWS 

NV  THE  NUMBER  OF  PARAMETERS  TO  BE  VARIED 

NT RACE  =  1  PRINT  OUT  MAP  OF  SEARCH  SEQUENCE 

=  0  NO  MAP 

X.  PARAMETER  ARRAY 

MASK  =0  SEARCH  ON  PARAMETER  I 

=1  HOcO  PARAMETER  I  FIXED 
OELTAX  INITIAL  STEP  SIZE  FOR  VARIABLE  XCI) 

OELMIN  SMALLEST  STEP  SIZE  ALLOWED  FOR  VARIABLE  X<I> 


CHISQ  RETURNED  FUNCTION  VALUE 


ALPHA=1.0  SMPLX  14 

BETA=0.5  SMPLX  15 

6AMMA=2. 0  SMPLX  16 

HUGE*1. 0  £37 
SIGNIF=1.E5 

NVP*NV*1  SMPLX  20 

LIM  IS  »  OF  EVALUATIONS  PERHITTEO  BEFORE  TERMINATION. 

00  10  J=1,NV  SHPLX  21 

10  Z(NV*1, J)=X( J)  SMPLX  22 

CALL  FUNK  SMPLX  23 

NF=1  SMPLX  24 

CHI(NV+1)=CHISQ  SMPLX  25 


CALCULATE  INITIAL  PCI)  AND  Y(I>. 


SMPLX  29 


00  30  J-l.NV 
00  20  K=1,NV 

20  Z(J|K)SX(K) 

SIGN  =  1.0 

23  IF (MASK  t  J)  .EQ.  1)  GO  TO  21 
Z(J,J)  =  Z ( J, J) ♦SIGN’OELTAX(J) 

21  XS=X(J) 

X(J)=Z(J, J) 

CALL  FUNK 
NF=NF»1 

IF  (  NF  .GE.  LIM)  GO  TO  600 

xcj)«xs 

IF(CHI(NV*1)-CHISQ)  22.30.30 

22  IF (SIGN  .EQ.  -2.0)  GO  TO  30 


SMPLX  31 
SMPLX  32 
SMPLX  3S 


SMPLX  36 
SMPLX  37 
SMPLX  38 

SMPLX  39 


..  MAX.W  r*tX-> 


A  *■'£“ 
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•  SIGN  *  -2.0 
GO  TO  23 
30  CHI ( J) =CHISQ 


SMPLX  40 


DETERMINE  H  AND  l. 


SHPLX  42 


40  JH=1 

SHPLX 

44 

JL  *1 

SHPLX 

45 

OO  80  J=2,NVP 

SMPLX 

46 

IF(CHI(J)-CHI(JH) ) 60, 60, 50 

SMPLX 

47 

50  JH= J 

SMPLX 

46 

60  IF(CHIU>-CHI(JL>  >70,80,80 

SMPLX 

45 

70  JL* J 

SMPLX 

50 

80  CONTINUE 

SHPLX 

*  * 

51 

CALCULATE  P9AR. 

SMPLX 

54 

90  DO  120  J=1,NV  SMPLX  56 

ZBAR ( J) -9 . 0  ‘  SMPLX  57 

00  HO  K=1,NVP  SMPLX  56 

IF(K-JH) 100,110,100  SMPLX  59 

10C  Z9AR(J>=2BAR(J)*Z(K,J)  SMPLX  6C 

11C  CONTINUE  SMPLX  61 

20  Z9AR(J)=Z9AR(J> /FLOAT  (NV)  SMPLX 


ATTEMPT  A  REFLECTION.  SMPLX  64 

FORM  P*  •  SMPLX  65 


DO  130  J=i , Ntf 

SMPLX 

67 

IF (MASK ( J)  .E9.  1>  GO  TO  110 

X(J>= ( 1.0+ ALPHA >»Z3AR(J)-ALPHA*Z<JH, J> 

SMPLX 

6  o 

ZST  AR ( J>  =  X ( J> 

SMPLX 

6 

CALL  FUNK 

SMOLr 

7 1 

NF*NF  *1 

SMPLX 

71 

IF  (  NF  .GE.  LIM)  GO  TO  600 

CHISTA  =CHISQ 

SMPLX 

72 

IF <CHIS3-C HI UL>  >14  0,180,180 

SMPLX 

73 

THE  REFLECTION  SUCCEEOEO.  ATTEMPT  AN  EXPANSION.  SMPLX  75 

FORM  P**  .  SMPLX  76 


140  IF(NTRACE) 150,160,150 


SHPLX  70 
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150  PRINT  43,  CHISQ 

43  FORMAT (1JX8HCHISQ  =  £15.  5, 10X2CHREFLLCTI3N  SUCC-'EOEO 
PRINT  44,  C* < J> ,J=1,NV> 

44  FORMAT t 5 X8HX<I> ....7(10(1X112.5)1) 

PRINT  155 

155  FORMAT ( 1H  > 

160  00  170  J=l,NY 

IFCMASK(J)  ,EQ.  1)  50  TO  170 
X(J)=GAMMA*X(J) •Ml. 0 -GAMMA) *ZBAR{ J) 

170  CONTINUE 
CALL  FUNK 
NF=NF*1 

IF  (  NF  .G£.  LIM)  GO  TO  600 
IF(CHISl-CrtKJL)  >  175,210,210 
175  IFtCHISa-CHISTA)  4C0, 210,210 


SMPLX  60 

SMPLX  62 

SHPLX  64 
SHPLX  65 


SMPLX  67 
SHPLX  66 


THE  REFLECTION  FAILEO. 

•  • 

SMPLX  91 

160  OC  200  J=1,NVP 

SMPLX  93 

IF(J-JH)190, 200,100 

SMPLX  94 

190  IF(CMISl-CHKJ)  )  230,200, 200 

SMPLX  95 

20C  CONTINUE 

SMPLX  96 

IF (CHISQ-CHI(JH) >250,270,270 

SMPLX  97 

210  IF( NT RACE) 220,230,220 

SMPLX  99 

220  PRINT  46 

46  FORMAT (43X1EHEXP ANSIQN  FAILEO  ) 

SMPLX 101 

REPLACE  P(JH>  BY  P*  • 

SHPLX103 

230  OO  240  Jsl,NY 

SMPLX10  5 

24C  Z(JH,J)=ZSTAR(J) 

SMPLX1C6 

CHI (JH)=CHISrA 

GO  TO  440 

SMPLX1C 6 

250  OO  260  J=1  ,NV 

SMPLX 110 

260  Z ( JH, J) =  X  ( J) 

SMPLX11 1 

CHI (JH)=CHIS1 

|  f  f 

SMPLX112 

ATTEMPT  A  CONTRACTION. 

SMPLX114 

FORM  P**  . 

»  ♦  t 

SHPLX115 

270  OO  280  J=1,NV  SMPLX117 

IF(MAS<(J)  .El.  1)  GO  TO  280 
X( J)=BETA*Z(JH, J)M1. 0-3ETA) *Z0AR(J) 

260  CONTINUE 


38 


non  oooo  o  oooooooo 


CALL  FUNK  SMPLX119 

NF*NF  *1  SMPLX120 

IF  (  NF  .  GE.  LIM)  GO  TO  600 

IF(CHISQ-CHI < JH> >290,330,330  SMPLX121 

290  IF ( NT RACE >  3f  0 , 42? ,  303  SHPLX122 

3CC  IF(CHISQ-CHI(JL)  >310,320,320  SMPLX123 

310  °RI  NT  333,  CHISQ 

303  FOPMAT ( 10  XBHCHISQ  *  £15. 3 , 10X21HCONTRACTION  SUCCEEOEO  )  SHPLX125 

PRINT  44,  (X(J>  ,  J  =  1,NV> 

PRINT  155 

GO  TO  420  SNPLX126 

320  PRINT  49 

49  FORMAT (43X21HOONTRACTION  SUCCEEOEO  >  SHPLX130 

GO  TO  420  SHPLX131 


THE  CONTRACTION  FAILED.  SMPLX133 

REPLACE  ALL  PtI)  3Y  (PC  I) *P ( JL> > /2.0  SMPLX134 


33C  IF(NTRACE) 340,350,340 
340  PRINT  56 

56  FORMAT (43X1BHCONTRACT ION  FAILED  > 
350  JS-JL 

00  390  J=1 ,NYP 
IF(J-JL)360,390,360 
360  OO  370  <=1,NV 

Z< J,K>=3.5*<Z< J,K>  +Z( JL, K>  > 

370  X(K)-Z(J,K) 

CALL  FUNK 
NF=NF  *1 

'  IF  <  NF  •  GE.  LIM)  GO  TO  600 
IF(CHISQ-CHKJL)  >  300,390,390 
360  JS*J 

390  CHI ( J) -CHISQ 
JL  =  JS 
GO  TO  440 


SHPLX136 

SHPLX136 

SMPLX139 

SMPLX140 

SMPLX141 

SMPLX142 

SMPLX143 

SMPLX144 

SMPLX145 

SMPLX146 

SMPLX147 

SMPLX146 

SMPLX149 

SMP'  X  15  0 
SMPLX15J 


400  IF( NT  RACE >410, 42 0,410  ShPLX 153 

410  PRINT  52,  CHISQ 

52  FORMAT (10X6HCHISQ  =  E15.  9 , 13X19HEXPANSION  SUCCEEOEO  >  SMPLX1 55 

PRINT  44,  (X(J> , J=1,NV> 

PRINT  155 


REPLACE  P(JM>  BY  P*  OR  P**  • 


SMPLX160 


420  OO  430  J- 1, NV 
430  Z ( JH, J) =X ( J) 
CHI(JH)=CHISQ 
C 


SMPLX161 

SMPLX162 

SMPLX163 
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TEST  FOR  CONVERGENCE. 


SMPLX165 


440  DO  470  J=i,NV  SMPLX167 

IF (MASK ( J )  .EQ.  1)  GO  TO  470 

ZMAX=-HUGE  SMPLX166 

ZMIN=HUGE  SMPLX169 

00  450  K=i,N VP  SMPLX170 

ZMAX=AM4X1(ZMAX,Z(K,  J>)  SHPLX 

re  ZMIN=AMIN1(ZMIN,Z(K, J))  SHPLX 

IF  ( SIGNI FMZMAX-ZMIN)- AMI  N1(ABS(ZM  AX  >,A3S(ZMIN)>  >470,470,460  SHPLX  - 

60  IF ( ZM AX" ZMIN- &9S ( OELMIN (J) ) >470,470,475 

470  CONTINUE  SHPLX175 


GO  TO  480 

475  JH=1 
JL  =  1 

DO  479  J=2,NVP 

IF(CHI(J)-CHI(JH> >  477,477,476 

476  JH-J 

477  IF(CHI(J)-CHIIJL> >  478,479,479 

478  JL= J 

479  CONTINUE 

IF(SIGNIF*(CHI( JH)-CHI(JL) )-CHI(JL>  >  480,480,90 


480  IF (CHI(JH)-CHI(JL) >490,500,500 
490  JL=JH 
500  CHISQ=CHI ( JL) 

90  510  J=1,NV 
510  X(J)-Z( JL, J) 

RFTURN 
OUT  OF  TIME 
600  CONTINUE 

PRINT  510,  NF,  LIH 

610  FORMAT (*3FUNCTlONAL  ITERATION  LIMIT  REACHED,  NF  =*,I5> 

RETURN 
END 

SUBROUTINE  FCN 

COMMON/FI Tl/NPT S.NPAR, YD ( 10 >0) ,ND(1050) , °Z (100)  ,  TO (1050) 

COMMON  /FITB/  FMIN ,FM AX , I T ITLE  (5 > 

COMMON  /FITC/  IPD ( 100) , IPSWCH 

COMMON/ST£PIT/NV,NTRACE,CHISO, MASK ( 100 > ,T(2) ,C(2>  ,R(2>  , 

1  A2 (2) , A1 ( 2) , 4  3 ( 2 )  ,  Mtl  (  2)  ,  9UM  (  86) ,  DEL  TAX  ( 130  ,OELMIN(100>  »NF 

COMMON  /BEST/  BX(100) , NF3EST , NRP3ST , BCHIS J 
COMMON/CONST1/Z1 , ZN1, PI , EYE 
COMMON/CONS  T*»/OELTAF,FZERO 
COMMON  /C0NST6/  NPLN, N? LS W , NPINT ,NE VAL 
COMMON  /DISP/  DSP  ( 2, J ) 

COMPLEX  K(2),CC(2) ,SINP(2) , COSP (2) , EYE 

COMPLEX  0231,3232,023,0321,0322,032,0221,0222,022,0331,0332,033 
DIMENSION  119(100) 

COMPLEX  V31 , VB2 , VB 
COMPLEX  K1 
REAL  K2 

NPINT  IS  HOW  OFTEN  TO  PLOT.  I.E.,  EVERY  50  EVALUATIONS,  ETC.. 

THIS  SUBROUTINE  IS  USED  TO  DETERMINE  THE  ADEQUACY  OF  THE  FIT 


SMPLX177 
SMPLX178 
SMPLX179 
SMPLX18  4 
SMPLX185 
SMPLX188 


SMPLX 
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ACHIEVED  BY  THE  CURRENT  PARAMETER  VALUES.  IT  PRODUCES  A 
CHI-SQUAREO  ESTIMATE  OF  TfclE  DIFFERENCE  BETWEEN  A  THEORETICAL 
MOOEL  UTILIZING  THE  CURRENT  PARAMETERS  ANO  THE  EXPERIMENTAL  DATA 
VALUES,  IT  IS  A  MODIFIED  mXLAYER“  SUBROUTINE  WHERE  THE  THIRD 
LAYER  IS  ASSUMED  TO  BE  IDENTICAL  TO  THE  FIRST. 

CHISQ=0.0 

INIT  PEAK  DETECTION 
N  =  1 

OON'T  PRINT  PEAK  MATCHES 
MPSWCH  a  0 

NEGATIVE  PARAMETERS  UNPHYSICAL 

IF  {  T (2 )  « LT,  0)  GO  TO  4 

IF  (  C(2)  .LT.  0)  GO  TO  4 

IF  (  R(2)  .LT.  0)  GO  TO  4 

IF  (  AO ( 2 )  .LT.  0)  GO  TO  4 

OO  1  I=1,NPTS 

F=< FLOAT (I)*DELTAF+FZERO) 

OMEGA=2.9*PI*F 
OO  6  J=l,2 

ATT=A2(J)*OMS5A*OMEGA«-Al(J)  *OM£GA>A0 (J) 

K(J)=C MPLX( OH EGA/C  (J) ,-ATT) 

KEEP  SQUARE-ROOT  ARG.  FROM  BLOWING  UP 
IF  C  CABS (  K(J>  )  .LT.  1.0E-99  )  PRINT  ll.F.J 
11  FORMAT  (*  K(J)  TOO  SHALL,  FREQ  =  *,E15.5,  *,J  =  *,I2) 

CC( J)=C3QRT( (OMEGA*OMEGA)/ (<( J)*K(J) )> 

K1  =  K(J)*T(J) 

KEEP  FUNCTION  FROM  BLOWING  UP  -  IMAGINARY  PART  MAG.  >300. 

IF  ( ABS ( AIMAG (K 1) )  .GT.  300.0)  PRINT  10,F,J 
10  FORMAT  (♦  K1  BLOW  UP,  FREQ  =  *,E15.5,*,  J=  *,I2) 

IF  (A0S(AIMAG(K1>  )  .GT.  300.0  )  GO  TO  4 
SlNPtJ)  =  CSIN(Kl) 

COSP(J)  =  CCOS(Kl) 

6  CONTINUE 

0231= (1.  0/ (R(2)*CC<2> ) ) * SINP ( 2) *C0SP <1) 

Q232= (1.  0/ (R(l>  *CC(1) ) ) *SINP(1)*C03P<2) 

Q23=-EYEMQ231«-Q232> 

Q321=R(2) *CC(2) *SINP(2) *COS?(l) 

Q322=R(1) *CC(1)*SINP(1)*C0SP(2) 

Q32=-EYE» (Q321+Q322) 

Q221=Q331=C03P(2)*C0SP(1) 

Q222= ( <R(1)*CC (1 ) I / (R(2) *CC  (2) ) ) *SINP(2) *SINP<1) 

Q332= ( (R(2)*CC(2) >/<R(l)*CC(l) ) )*SINP(2)*SINP<1> 

Q22=Q221-Q222 

Q33=Q331-Q332 

V81=C03P(1)M  <EYE*ZN1)/(R(1)*CC(1) )  )*SINP(1) 

VB2  =  ZN1*C0SP(1> ♦EYE*R(1) *  CO  ( 1 ) *SINPtl> 

VB=VB1*{Q32-Z1*Q33)-V32* (Q22-Z1*Q23> 

Y  =  CABS ( (-?• 0*Z1)/V3) 

STORE  THEORETICAL  RESULTS 
NORMALIZE  THEORETICAL 

Y  =  WN ( 2 )  * Y 
TO ( I)  =  Y 

FMIN  IS  MINIMUM  FIT  FREQUENCY 
FMAX  IS  MAXIMJM  FIT  FREQUENCY 
IF  TF  .LT.  FMIN)  GO  TO  1 
IF  (F  .GT.  FMAX)  GO  TO  1 
RES=YD(I)-Y 
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CHISQ=CMISQMRSS*RES*HO(I)) 

1  CONTINUE  *  '/ 

GO  TO  5 

C  UNPHYSICAL  PARAMETER 

4  CHISQ  =  1E1C 

5  CONTINUE 

•C  COMPARE  PEAKS? 

CHFACT  =  CHISO 
PHEIGH  s  0.0 

MATCH  »  0  i 

IF  (IPSWCH  .NE.  1)  GO  TO  29 
C  GET  THEORETICAL  PEAKS 

CALL  PEAK(TO,NPTS,ITD,0> 

N=1 
H  »  1 

C  SEE  IF  PEAKS  MATCH 

DO  21  H  =  1.100 

IF  (IABS  (IPD(N)-ITD(M))  . LT .  10)  MATCH*1 
IF  (MATCH  .EO.  1)  PRINT  31 
IF  (  MATCH  • E3.  1)  GO  TO  22 

31  FORMAT  (*  THEORETICAL  l  EXPERIMENTAL  PEAKS  MATCH*) 

21  CONTINUE 

C  NO  HATCH 

GO  TO  29 

C  MATCH  -  COMPUTE  PEAK  AMPLITUDE  OIFF 

22  00  24  I  =  N,  100 

C  END  LOOP  IF  OUT  OF  PEAKS 

IF  (  IPO ( I)  .GE.  NPTS)  GO  TO  29 
F  =  FLOAT  (IPKI)  >*OELTAF  ♦  FZERO 
C  DON'T  COMPARE  IF  OUT  OF  RANGE 

IF  (  F  .LT.  FMIN)  GO  TO  24 

IF  (  F  .GT.  FMAX)  GO  TO  24 

C  SEE  IF  EXTRA  PEAK 

26  CONTINUE 

IF  (M*l  .GT.  100)  GO  TO  23 

11  =  IA3S  (IP0(I)-IT0(M«-1)> 

12  =  IA3S  (I°D ( I ) -IT3 (M) ) 

IF  (  II  • LT.  12)  M  =  M*1 
IF  (  II  .LT.  12)  GO  TO  26 

IF  (  12  .GT.  10)  PRINT  34,  IPO(I) 

34  FORMAT  (*  REJECT  DATA  PEAK  INDEX)*, I4> 

IF  (  12  .GT.  10)  GO  TO  24 

23  CONTINUE 

C  PRINT  PEAK  MATCH-UPS  IF  MPSHCH  =  1 

IF  (MPSHCH  • EQ.  0)  GO  TO  27 
PRINT  33,  IPD(I),  ITO (M) 

33  FORMAT ( *  MATCH  OATA  PEAK  INDEX)*, 14,*  WITH  THEORETICAL «*, 14) 

27  CONTINUE 

PPES  =  YDdPO(I)  )-TO(ITO(M)) 

PHEIGH  =  PHEIGH  MPRES*PRES) 

M  e  H*1 

IF  (  M  .GT.  100)  GO  TO  25 

24  CONTINUE 
29  CONTINUE 

PHEIGH  »  PHf I5H*100.0 

25  CHISO  *  CHISO  ♦  PHEIGH 
PRINT  32,  PHEIGH 

32  FORMAT (•  PHEIGMT  =  *,  £15. 5) 
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PRINT  2,  NF*  CHISQ 

2  FORMAT (20 X, ’EVALUATION  NUMBER  *,I5,*  CHI  SQUAREO  *  ’,E13.5> 
PP1NT  3, tT(I),C(I) ,R(I) ,42(I> , A1 (I ) , AC (I)  ,  MN ( I) , 1*1, 2» 

3  FORMAT (11X,7£15«5) 

«»LOT  EVERT  NPINT  EVALUATIONS* 

NEW  (JESTFIT?  • 

IF  (CHISO  *LT*  9CHIS3)  CALL  BFIT 
IF  (  NF  • £Q«  0)  GO  TO  15 
X  *  FLOAT ( NF) /NPINT 

IF  (  X-(MF/NPINT)  *EQ.  0*0)  CALL  PICTR(l) 

15  RETURN 
ENO 

SUBROUTINE  INOATA 
INTEGER  TRN(4» 

COMPLEX  EYE* 0( 1950) 

COMMON/FIT1/NPTS,NPAR,YO(105C) ,WO(1050) ,PZ(100) ,TD(1050>. 

COMMON  /FITB/  FMIN , FM AX, IT ITLE <5> 

COMMONEST EPIT/NV ,NTR4CE, CHISQ, MASK (ICC) *T  ( 2) ,C(2) ,R(2) » 

1  A2  (  2)  •  41  (  2  >  *  AO  (  2 )  *  WN  ( 3)  ,  0UH2  (  3)  ,  0UM3(  3>  »0UM<*(3)  * 

2  0UM5(3>  *  DU (6 (3)  ,  DUM7 ( 3) , OUM9 (67) ,OELT AX (100 ) , 

3  OELMIN(IGO) ,NF 
COMMON/CONST1/Z1 , ZN1, PI,£Y£ 

COMMON/CONST4/OELTAF,FZERO 

THIS  SUBROUTINE  INPUTS  THE  EXPERIMENTAL  OATA  AND  TRANSDUCER 
RESPONSE  FROM  FILE  TAPE6*  AND  ALSO  INPUTS  INITIAL  LAYER  PARAMETERS 
FROM  OATA  CARDS. 


INITIALIZE  NUMBER  OF  PARAMETERS 
NPAR  =  14 

INITIALIZE  NUMBER  OF  PARAMETERS  NOT  VARIED 
NV=NPAR 


TURN  OFF  TRACE  OF  SEARCH  PROCEDURE 
NTRACE*0 

INITIALIZE  CONSTANTS 
LYE=CMPLX(0.  3,1.0) 

PI=3, 141592654 

IMPEOANCE  AT  NORMAL  INCIDENCE  -  FIRST  MEOIUM 
Zl=l. 4BE05 

IMPEOANCE  AT  NORMAL  INCIDENCE  -  LAST  MEOIUM 
ZN1S1.40E05 

READ  OATA  SET. 

REA0(6*1I (TRN(I) *  1  =  1*4) 

1  FORMAT (4110) 

PRINT  10 

10  FORMAT  (*-*, 10X, ’DATA  SET*! 

PRINT  11,  TRN(l) , T RN ( 3 ) 

11  FORMAT  ( * 0*, 10X ,* OAT 4  RECORO  «** 13 » • ,I3,*)*  I 
PRINT  12,  TRN(2) ,  TRN(4» 

12  FORMAT  (*0’,10X,  ’TRANSDUCER  RC.CORO  (’,13,’,’,  13,’)’  > 
READ(6,’)NPTS,FZ£R0, OELT  AF,TI 


J 
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1 


RFAO  COMPLEX  TRANSMISSION  FUNCTION 
RE  A  0  (  &  1 5  )  ( D(  I  )•,  1  =  1  ,Nr*TS) 

5  F0RMAT(2E 15.51 

REAO  PARAMETERS 

RE  AO  2,  (T(I),C(I),R(I>,A2(I)  ,A1(I) , AO (I) ,  MN(I),I=1»2> 

2  FORMAT ( 7E 10 . 2) 

STORE  PARAMETERS  IN  PZ  ARRAY 

PZ(1)*T(1) 

PZ(2)=T(2) 

PZ(3)=C(1> 

PZ(4>=C(2» 

PZ (5) =R(1) 

PZ  (6)  =R(  2) 

PZ (7) =A2  tl) 

PZ ( 8) =A2 (2) 

'  PZ  ( 9 ) =A1 ( 1) 

PZ (10 )=A1 (2) 

PZ(11 ) =A0 (1) 

PZ ( 12) =A0 ( 2) 

SFT  UP  NORMALIZING  PARAMETER 
PZ(13)  =  WN ( 1 )  =  0.0 
DEFAULT  VALUE  OF  HN(2)  =  1 
IF  (WN<? )  .EQ.  0.0)  HN( 2 )  =  1.0 
PZ ( 14  >  =  WN  ( 2) 

USE  ONLY  REAL  PART  OF  COMPLEX  TRANSMISSION  FUNCTION 
OO  7  1=1 i NPTS 
YD ( I) =  CA3S(0( I) ) 

7  CONTINUE 

INITIALIZE  WEIGHTING  FUNCTION 
00  4  I  =1,NPTS 
WO ( I )  =1 
4  CONTINUE 

FIX  FIRST  LAYER  PARAMETERS 

MAS  K(  1 ) =MASK ( 3) =M ASK( 5)  =  M ASK  ( 7) =MASK (9) =  MASK (11 ) =MASK (1 3) =1 

READ  SECONO  LAYER  MASK, FM IN , FMAX , ITITLE  l  WEIGHTING  FLAG 

RE  AO  6, (MASK ( I)  , 1=2, MPAR, 2)  ,  FM  IN , FMAX , ( ITI TLE (I ) ,1  =  1,4) , IWFLAG 

6  FORMAT  (7I1,3X,2E10.2,<*A10,9X,I1) 

SET  UP  DEFAULT  VALUE  FOR  FMAX 

IF  (FMAX  .LE.  0.0)  FMAX=20.0E06 
ENO  TITLE 
ITITLE (5 ) =  10H$ 

PRINT  17,  FMIN,  FMAX 

17  FORMAT (•  * , 10X , *  FMIN  =  *,E15.5,*  FMAX  =  *,E15.5> 

SEE  IF  TRANS0U3ER  RESPONSE  CURVE  TO  BE  USED  AS  WEIGHTING  FUNCTION. 
IF  (IWFLAG  .E3.  1)  PRINT  15 
15  FORMAT  (*-*,10X, ’WEIGHTING  FUNCTION  =  1*) 

IF  (  IWFLAG  .EQ.  1)  GO  TO  8 


c 


RFAO  TRANSDUCER  RESPONSE  CURVE 
*rAD(6,l> (TRN(I) ,1=1,41 
RFA0(6,*)NPTS,F7£R0,DELTAF,TI 
REA0(6, 5) (OCX) »I=1,NPTS> 

C 

c 

PRINT  IS,  TRN<1> ,TRN(3> 

16  FORMAT(*-*»19X, ’WEIGHTING  FUNCTION  RECORD  C* , 13 , 13, *> *  ) 
C  USE  ABSOLUTE  MAGNITUDE  FOR  WEIGHTING  FUNCTION 
00  3  1=1, NPTS 
W0(I)=CA3S(D(I)) 

3  CONTINUE 

c 

6  RETURN 
ENO 

SUBROUTINE  PICTR(LOOP) 

COMMON  /CONSTL/  BELT AF, FZERO 

COMMON/FIT1/NPTS,NPAR,VO(1C50) , WO (1950)  ,  PZ ( 10 0) , TO C 1050) 
COMMON  /FITE/  FLO, FHI , IT  I TLE ( 5) 

DIMENSION  X  ( 1950 )  ,  Y(1C50>,  Y2(105Q) 

COMMON  /FITO/  3 ( 2 1  , 02 (2) , 3 1(2) , 30 ( 2) 

C  BEST  FIT  PARAMETERS 

COMMON  /STEPIT/  NV , NT R AC" , SCHISQ, MASK( 19 0)  , 

1  XX  (IOC)  .OSLTAXUOC)  ,DELMlN(ldO>  ,SNF 

COMMON  /BEST/  THICK( 2) , C ( 2 ) ,ROW (2) , C2 (21  ,C1 (2) ,C0 (2) , WN (2) , 

1  DUM(b6) ,NT,NRP3ST,CHISQ 

DIMENSION  LBLX(IO) ,LBLY(19) , LBLH (10 ) ,LST RG ( 10) , LPAK ( 150 ) 
COMPLEX  VSV,OSO 
COMPLEX  0(1024) , V ( 1024) 

COMMON  /CONST5/  NPLN,NPLSW,NPINT,NEVAL 

NLAYER=2 

TI  *  0.0 

RTO=57.29577951 

LBLX( 1) = 10HFRE QUENCY 

LBL X ( 2) s 1 OHI N  MHZ. S 

LBLH( 1) =10H  T  ( 

LBLH(2)=10HR>  C 
LBLH ( 3) =10H  9 

LBLH ( 4) = 10H  C2 

LBLH(5) =10H  Cl 
LBL  H(6) =10H  CO 
LBLH( 7) * 10H  92 
L9LH( S) *10HB1  B 

LBLH(9)=10HO$ 

XL=10.0 

YL=5.0 

XSF=XL/19.0 

YSF=YL/5.0 

HITE=  0. 1 2*YSF 

WIDTH=HITE*6. 0/7.0 

XLEFT=0. 25*XSF 

YBTM=YL+0.50*YSF 

XLC-0.0 

XRC=95. 3*WIOTH’2.0*XLEFT 
YBC*VBTH-0.25*YSF 

YTC=YBC^0. 5C*YSF+1.35*HITE* (FLOAT (NLAYER)^l.O) 

C  SET  UP  HORIZONTAL  SPACING  FOR  PARAMETERS  IN  BOX. 

XJUHP=9»0* WIDTH 
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XJO=XLEFT 

XJ1*XJ0»7.C*WI9TH. 

XJ2=X Jl*5.  0* WIDTH 
XJ3=XJ2+X JUMP 
XJ4=X J3*X JUMP 
XJ5=XJ4+XJUMP 
XJ6  =  X J5  +  X  JUMP 
XJ7=XJ5+X JUMP 
XJB=XJ7*XJUMP 
XJ9=XJ8*XJUMP 
PRINT  1,  LOOP 
1  FORMAT  (*OLOOP  *  *,  15) 

GO  TO (6, 7 ,8, 9, 10, 11, 12, 13, 21) LOOP 

6  CONTINUE 

00  5  1=1, NPTS 

X(I)  =  (FLOAT  (I)  *DELTAF«-FZERO)  /l*  0E06 
THEORETICAL  3ATA 
Y( I )  =  T0(I) 

EXPERIMENTAL 
Y2(I)  =  YO (I) 

5  CONTINUE 

LBLY (1)=10HTRANSMISSI 
LBLY(2)=10HON  COEFFIC 
LBLY (3)=10HIENTS 
GO  TO  100 

7  CONTINUE 

00  14  Is 1 i NPTS 

X(I)= (FLOAT (I) *DELTAF+FZSRO)/1.0E06 
Y(I)  =  CABS(V(D) 

14  CONTINUE 

LBLY (1)=10HPEFLECT ION 
LBLY ( 2) =10 H  COEFFICIE 
LBLY(3)=10HIENTS 
GO  TO  100 

8  CONTINUE 

00  15  1=1, NPTS 

X(I)= (FLOAT (I) *DELTAF+FZERO)/1.0E06 
YCI)=RTO*  ATA-HAIMAGO(I)  > /REAL  (0  (I ) )  ) 

IF(REALO(I)).LT.O.O.  AND.  AIM  AG  (9  ( I )  ) .  GT.  9.  0  )  Y  ( I )  =Y  (  X) +180. 0 
IF  (REAL  (0(D)  .LT.0.0.  AND.  AII1AG  (0  ( I ) )  .LT  .  0.  0  )  Y  ( I )  =Y  ( I) -1 80. 0 

15  CONTINUE 

LBLYC1) =10HPHASE  OF  D 
LBLY(2)=10H  IN  OEG.S 
GO  TO  100 

9  CONTINUE 

00  16  1=1, NPTS 

X(I>= (FLOAT <I)*DELTAF*FZERO> /1.0E06 
Y<I)=RT9*AT#N(AIMAG(VU))/REAL(Y(I))) 

IF ( RE  AL ( \M  I ) ) .LT.0.0. A. 40. A I  MAG  ( V ( I ) ) . GT . 0 .  C ) Y ( I ) = Y ( I ) *1 80 . 0 
IF(REAL(Y(I) ).LT.0.C. ANO. A IMAG  ( V ( I ) ) . LT .  9.  0) Y (I ) =Y ( I) -180, 0 

16  CONTINUE 

LRL  Y  ( 1)  =10HPr)ASE  OF  V 
LBLY ( 2)  =  1  OH  IN  OEG.S 
GO  TO  100 
lb  CONTINUE 

00  17  1=1, NPTS 

X(I>= (FLOATCI)*OELTAF»F2ERO)/1.0E06 
Y(I)  =  REAL  (0(D) 


17  CONTINUE 
LBLY(1)=10HR£AL  COMPO 
LBLY(2)=10HNENT  OF  OJ 
GO  TO  100 

11  CONTINUE 

00  18  1=1, NPTS 

X(I)= (FLOAT(I)*DELTAF+FZERO>/1.0E06 
Y(I)=AIMAG(D(I)) 

18  CONTINUE 
LBLY(1)=10HIMAGINARY 
L3LY(2)=10HCOMPONENT 
LBLY ( 3) =10HOF  0$ 

GO  TO  100 

12  CONTINUE 

DO  19  I=1»NPTS 

X(I)  = (FLOAT (I) *OELTAF*FZSRO)/1.0E06 
Y ( I )  =  REAL ( V (I >  > 

19  CONTINUE 
LBLY(1)=10HR£AL  COMPO 
LBLY (2) =1 OHNENT  OF  tf$ 

GO  TO  100 

13  CONTINUE 

00  20  1  =  1 , NPTS 

X(I>=(FLOAT(I)*OELTAF+FZERO>/1.0EO& 

Y(I)  =  AIMAG(V(I)  ) 

20  CONTINUE 

L9LY ( 1) =10HIM AGINARY 
LBLY(2)=10HCOMPONENT 
L3LY(3)=10HOF  V$ 

GO  TO  100 

21  CONTINUE 

00  22  1=1, NPTS 

X(I)=(FLOAT(I)*OELTAF*FZERO)/i.OE06 
VSV=V(I) *CONJG  (V ( I ) ) 

DS0=0(I)*C0NJG(0(I)» 

Y(I)=REAL(VSV*DSO> 

22  CONTINUE 
LBLY(1)=10HFNERGY  CON 
L9LY(2)  =  13HSERVATION  " 

LBLY(3>  =  10H(YS\/*OSO)S 

100  CONTINUE 

FMIN=FZERO/1.0E06 

FMAX= (FLOAT (NPTS >*DELTAF+FZERO)/1.0E06 
FSTEP=(FMAX-FHIN)/10.  0 
NPLN=NPLN*1 
PRINT  2,  NPLN 

2  FORMAT  (*0  CALL  3GNPL,  NPLN  =  *,I5> 

CALL  BGNPL (NPLN) 

XFO=(11.0-XD*Q.75 
YFO=(9.53-YU*0.29 
CALL  PHYSOR(XFO, YFO) 

CALL  XTICKS ( 10 ) 

CALL  YTICKS(4) 

CALL  TITLE (1 H  , -1 ,LBLX, IP  0 ,LBL Y, in 0 , XL ,  YL) 
IF(TI.5T.0,C)CALL  M£SSAG(”CM  ACOUSTIC  SPECTRUM  THETA 
12.25*XSF,7.1*VSF) 

IF  (TI.GT.O.O)CALL  REALNO  (TI,  1,“A9UT~,-A9UT"> 
IF(TI.GT.0.0)CALL  MESSAGC*  OEG.f,  100, ”A3UTM, “ABUT-) 
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IF(TI.EQ. 0.0>CALL  MESSAG ( “CM  ACOUSTIC  SPECTRUM  NORMAL  INCIOENC 
1ES",100,2.25*XSF,7.1*YSF) 

C  PRINTI  EVALUATION  V,  CHISQ 

CALL  MES3AG  ("EVALUATION  NO.  S”» 10 0 , 2.  25*XSF, 6. 8* YSF) 

CALL  INT  NO (NF, "A  OUT”, "ABU  T"> 

CALL  MES3AG ("  CHISQ=  S", 10C, "ABUT", "AOUT") 

CALL  REALNO(CHISO, 103, "AOUT”, “A3UT") 

CALL  MESSAG(ITITLE,1)G,2.  25*XSF,6.4*YSF) 

CALL  MESSAG(“FMIN  =  S“, 1 0 0 , 0 . 0 , -. 9*YSF> 

CALL  REAlNO  (FLO,  13 3 ,"A JUT", "A3UT") 

CALL  MESSAG("  , FMAX  =$“, 13 0  ,  "A0UT"f"A0UT“> 

CALL  REALNO  (FHI , 103 , "A  OUT", "ABUT") 

CALL  MESSAG  <"  ,WN(2>  =  3“ , 1G 0 ,"ABUT","A3UT“> 

CALL  REALNO  ( MN ( 2) ,103,"A3UT”»"A9UT“) 

CALL  8LNK1(XLC,XRC,YJC,YTC,2> 

IF (LOOP. EQ. 3. OR. LOOP. EQ. 4 >  CALL  GRAF ( FMIN, FSTEP, FMAX ,-180. 0 ,60. 0, 18 
10.0) 

IF  (LOOP.  GT.4. ANO.LOOP.LT. 9) CALL  GRAF (FMI N, FSTEP, FMAX ,-l . 0 ,0 . 2, 1. 0) 
IF  (LOOP  .EQ.  1  .OR.  LOOP  . EQ.  2  .OR.  LOOP  .EQ.  9) 

1  CALL  GRAF  (FMIN, FSTEP, FMAX, 3.0, "SCALE", 1.0) 

IF (LOOP, GT.2.AN0.L003.LT. 9) CALL  XSTAXS (FMIN , FSTEP, XL , 1H  ,1,0.0,YL/ 
12.0) 

CALL  MIXALF("L/CGREEK”) 

CALL  HEIGHT (HITE) 

OYM=1.35*HITE 

YMl=YBTM*FLOAT (NLAYER)’DYM 
CALL  RESET ("BLNK1") 

CALL  MESSAG(LBLH,100,XLEFT,YML) 

00  3  I=1,NLAYER 

C  SKIP  INITIAL  LAYER  IF  ZERO  THICKNESS 
IF  (  THICK ( I )  .EG.  0.0)  GO  TO  3 
YML=YML-OYM 

CALL  REALNO (THICK (I) ,105, X JO, YML) 

CALL  REAL  NO (ROM (I) ,1Q3,XJ1, "A9UT") 

CALL  REALNO (C (I ) , 10 3 , XJ2 , "ABUT”) 

CALL  REALNO (3(1)  ,  1 G3, X J3 , "ABUT") 

CALL  REALNO(C2 (I)  ,  -2 , X J4 , “ABUT") 

CALL  REALNO (Cl ( I )  ,  -2,XJ5,"A9UT") 

CALL  REALNO (CO ( I ) ,  -2 ,XJ6 ,"ABUT") 

CALL  RSALNO( 32(1),  - 2, X J7 , " ABUT") 

CALL  REALNOOKI),  -2  »XJ3  ,“AOUT") 

CALL  REALNOOO  (I)  ,  -2,XJ9,"ABJT") 

3  CONTINUE 

C  INDICATE  "FIXED"  PARAMETERS  IN  BOX 
YML  =  YML-OYM 
W=MIOTH 

IF(MASK(2).EQ.1>  CALL  MESSAG  ("FIXED*",  100,  XJ04-W,  YML) 

IF(MASKCj)  .EQ.l)  CALL  MESSAGC’FXOS"  ,193,XJ1  +  H,YML) 

IF(MASK(4) .EQ.l)  CALL  MES S AG  (“FIXE J 3”, 13 0  ,  X J2tW , YML) 

IF(MASK(3) .EQ.l)  CALL  MESS AG ( "FI XED5" , 13 0 , X J4»W , YML ) 

IF  (MASK (10) . EQ.l ) CALL  MES SAG ( "FI XEO S“, 10 0  ,  X J5+H , YML ) 
IF(MAS<(12).EQ. 1 ) CALL  MESSAG ("FIXEOJ", 100,  XJ6  +  M, YML) 

CALL  SIMPLX 

CALL  BASALF  ("L/CSCRIPT") 

CALL  MIXALF  ("SCRIPT") 

CALL  MESSAG  (“(D)ASHZD  =  (T)HEORETICAL,  (S)OLIO  *  (E) XPERIMENTALS" 
1  ,100,2.25*XSF,5.00*YSF) 

CALL  CARTOG 
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o  o  nnonn 


CALI  RESET  (“BASALF") 

IF  (NPLSW  .GT.  1)  GO  TO  \ 

C  PLOT  ACTUAL  DATA  -  SOLID  LINE 
CALL  CURVEU,  Y2,NPTS,0) 

U  CONTINUE 

C  PLOT  THEORETICAL  CURVE  -  SOLID  LINE 
CALL  OASH 

CALL  C'JRVE(X,Y,NPTS,0) 

CALL  RESET ("DASH“» 

CALL  ENOPL1NPLN) 

RETURN 
ENO 

SUBROUTINE  PEAK(Y,NPTS,IPD,NPRINT) 

THIS  SUBROUTINE  DETERMINES  THE  LOCATION  OF  DATA  PEAKS 
IT  STORES  THE  PEAK  INDICES  IN  IPO  UP  TO  lOt  PEAKS. 

IF  NPRINT  =  1,  PEAK  INFORMATION  IS  PRINTED 

DIMENSION  Y(NPTS),  I»0(100) 

COMMON  /CONST<»/  OElTAF.FZcRO 
INITIALIZE 

NFLAG  INDICATES  POSITIVE  SLOPE  IF  *  1 
YMAX  -0.0 
YMIN=  1.0 
NFLAG  =  0 
J  -  i 

IPO(l)  =  NPTS»1 
00  11  I  =  l.NPTS 
C  REJECT  INITIAL  GARBAGE 


IF 

(  Y ( I )  .GT. 

0.95) 

GO  TO 

10 

IF 

(Y(I)  .GT. 

YMAX 

.AND. 

NFLAG  .EQ. 

1) 

YMAX  = 

Y(I) 

FOUND  PEAK 

IF 

(  Y  ( I)  .LT. 

YMAX 

.ANO. 

NFLAG  .EQ. 

1) 

GO  TO 

5 

IF 

CYCI)  .LT. 

YMIN 

.AND. 

NFLAG  .EQ. 

0) 

YMIN  * 

V<I> 

C  POSITIVE  SLOPE? 

IF  (  Y ( I )  .EQ.  YMIN)  GO  TO  10 

C  SEE  IF  STILL  ON  DOWNTREND 

ki  =  m 
K2  *  1*5 

IF  (K2  .GT.  NPTS)  GO  TO  10 
C  SCAN  NEXT  10  POINTS 

DO  3  K  =  K1.K2 

IF  (  YMIN  .GT.  Y  (K) )  GO  TO  10 
3  CONTINUE 
NFL AG=1 
GO  TO  10 

C  FOUNO  PEAK  -  STORE  INOEX 

5  CONTINUE 

C  LOOK-AHEAD.  MAKE  SURE  YMAX  *  PEAK 

Kl  a  I 
K2  =  I  *-20 

IF  (K2  .GT.  NPTS)  K2  a  NPTS 
C  SCAN  NEXT  20  POINTS 

00  7  K  =  K1,K2 

C  YMAX  REALLY  PEAK  OR  GLITCH? 

IF  (  YMAX  .LT.  Y (K>  )  GO  TO  10 

7  CONTINUE 

C  YMAX  REALLY  PEAK 
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oooooooo 


IF  (J  .Lit  100)  IPO(J)*  1-1 
J  =  J*1 

IF ( J  .LE.  100)  IPO(J)  *  NPTS*1 
IF  (NPRIMT  .LT.  1)  GO  TO  25 
F  =  <FLOAT(I-l)*OELTAF  ♦  FZERO) /l. 0 E06 
PRINT  23,  F,  YMAX,  IOO(J-I) 

20  FORMAT  (•  FOUNO  PEAK,  FREQUENCY  =  *^0.3,*  MHZ,  AMPLITUOE  «*,F0.3, 
1  •,  INOEX  =  •,  15) 

25  CONTINUE 
YMIN  =  Y ( I) 

NFL  AG  =  0 
YMAX  =  0.0 

10  CONTINUE 

IF  (NPRINT  .EQ.  2)  PRINT  8 , 1 , YMIN, YM AX , Y (I ) ,NFL AG 
0  FORMAT  <*  I=*,I5,*,YMIN=  *,El5.5,*,YMAX=  • , E15. 5  ,*  ,  Y»* ,E15.5 , 

1  *,  NFLAG  =  *,  13) 

11  CONTINUE 
RETURN 
ENO 

SU9R0UTINE  B" IT 

COMMON 7STEPIT /NV , NTRACE , CH ISQ , MASK (100)  , 

1  X (100) , OELTAX  (IOC) »Q£LMIN(100)  ,NF 

COMMON  /BEST/  BX(100) ,NFBEST , NRPOST , BCHIS3 

THIS  SUBROUTINE  STORES  THE  PARAMETERS  RELATED  TO  THE  BEST  FIT 


BX  -  PARAMETER  ARRAY 

NFBEST  -  #  OF  THE  FUNCTIONAL  EVALUATION 
NRPBST  -  #  OF  DATA  PEAKS  REJECTED 
9CHISQ  -  BEST  CHI  SQUAREO 


00  10  I  =  1 ,  NV 
BX(I)  =  X(I) 

1C  CONTINUE 
NFBEST  =  NF 
8CHISQ  =  CHISQ 
RETURN 
END 

0.635. 64327E+ 06 
f • 3127  1.000E05 

1011101  3.0E06 


2.70  0.3 

1,07  1.64 

15. 0 E 0 6 A L/ ADHESIVE/ AL  SANDWICH-T6R035- VARY  C 
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MEASUREMENT  OF  THE  RADIAL-BOLT-HOLE  CRACK 
SIZE  USING  THE  DOPPLER- SHIFT  TECHNIQUE 


by 


N.  K.  Batra 


ABSTRACT 


The  Doppler  effect  can  be  used  to  estimate  the  size  of  radial-bolt-hole  cracks 
in  mechanically  fastened  aircraft  structures.  A  sample  having  a  radial  crack 
is  mounted  on  a  rotating  platform  in  a  water  bath  and  Interrogated  using  a 
train  of  tone  bursts  from  a  stationary,  focused  transducer.  The  Doppler-shift 
spectrum  of  the  backscattered  ultrasonic  beam  from  the  rotating  crack  is 
recorded.  The  maximum  Doppler  shift  is  then  related  to  the  radial  size  of 
the  crack.  The  merits  and  limitations  of  this  technique  are  discussed. 
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MEASUREMENT  OF  THE  RADIAL-BOLT-HOLE  CRACK 
SIZE  USING  THE  DOPPLER- SHIFT  TECHNIQUE 


by 


N.  K. 


* 

Batra 


I.  INTRODUCTION 

A  aajor  problem  in  extending  the  operational  life  of  aircraft  Is  the  detection 
and  measurement  of  the  size  of  fatigue  cracks  In  highly  stressed  areas  of  the 
airframe.  The  cracks  that  present  the  most  serious  threat  to  flight  safety 
are  those  which  grow  from  the  bolt  holes  In  mechanically  fastened  structures. 
These  flaws  are  among  the  most  difficult  to  detect  for  several  reasons.  First, 
these  fatigue  cracks  most  often  occur  under  the  heads  of  counter-sunk  bolts 
at  the  base  of  the  bolt  holes.  Second,  the  use  of  new  high-strength  structural 
alloys  has  necessitated  the  detection  and  size  measurement  of  very  small  flaws 
in  these  crack-intolerant  materials.  At  present  the  inspection  procedure  for 
detecting  these  defects  consists  of  removing  the  fastener  and  inserting  an 
eddy-current  probe  Into  the  hole.  Since  the  average  cost  of  this  Inspection 
is  approximately  $150. 00/bolt  and  large  cargo  aircraft  may  contain  20,000 
critical  fasteners,  the  economic  incentive  to  develop  semi-automatic  Instru¬ 
mentation  Is  readily  apparent. 

^Systems  Research  Laboratories,  Inc.,  2800  Indian  Ripple  Road,  Dayton,  OH  45440. 
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Although  several  commercial  lnstrumenta  have  been  developed  to  detect  these 
flaws,  they  do  not  have  the  capability  to  measure  the  size  of  the  flaw.  How¬ 
ever,  crack  size  has  been  estimated  using  a  number  of  ultrasonic  techniques^ 
involving  amplitude, ^  time, ^  and  acoustic-spectrum  methods. Methods 
of  size  estimation  which  make  use  of  the  amplitude  of  the  reflected  sound  wave 
are  inaccurate  for  many  reasons,  the  principal  one  being  the  differing  reflec¬ 
tivities  of  fatigue  cracks.  Such  methods  also  disregard  several  other  impor¬ 
tant  factors  such  as  crack  orientation,  shape,  and  quality  of  the  acoustic 

(3) 

coupling — all  of  which  affect  measurement  accuracy.  Timing  methods '  '  have 
been  used  to  estimate  crack  dimensions  from  the  time  delay  of  a  signal  return¬ 
ing  from  a  crack  relative  to  the  excitation  pulse,  or  some  known  scattering 
point  on  the  fastener  hole.  These  methods  are  cumbersome,  and  the  best  results 
are  obtained  for  very  large  cracks  such  as  those  found  in  welded  steel  struc¬ 
tures. 

Because  of  the  circular  symmetry  of  the  bolt-hole,  a  method  based  upon  the 
Doppler  Effect  has  been  developed  to  measure  radial  crack  size.  This  method 
does  not  depend  upon  the  amplitude  or  time  delay  of  the  reflected  signal  from 
the  crack.  The  maximum  dimension  of  the  crack  is  proportional  to  the  maximum 
Doppler  shift  of  the  backscattered  signal  from  the  crack.  The  method  is 
potentially  more  accurate  than  those  used  previously. 


II.  THEORETICAL  DEVELOPMENT 

The  classical  Doppler  effect  refers  to  the  change  in  frequency  which  occurs 
when  a  sound  wave  is  backscattered  (reflected)  from  a  discontinuity  moving 
in  the  medium  of  the  Incident  wave.  A  Doppler  shift  of  the  sound  wave  also 
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occurs  when  the  veve  enter*  e  medium  which  1*  In  motion  relative  to  Its  source. 
The  magnitude  of  the  frequency  shift  can  be  used  to  estimate  the  size  of  a 
crack  or  reflector  in  the  moving  medium.  In  this  section  the  equation  used 
to  estimate  crack  size  is  developed. 

The  schematic  representation  of  the  experiment  Is  shown  in  Fig.  1.  An  ultra¬ 
sonic  transducer  is  in  a  water  bath  along  with  an  aluminum  plate  rotating 
with  angular  velocity  £1  In  the  x^,  x2  plane.  The  wave-nunber  vector.  It,  of 
the  incident  longitudinal  wave  is  shown  in  the  x^,  x^  plane.  Several  points 
concerning  this  physicsl  arrangement  should  be  noted.  First,  a  Doppler  shift 
in  the  frequency  of  the  incident  sound  wave  occurs  both  when  the  wave  enters 
and  exits  the  sample  at  the  water /aluminum  boundary.  Second,  since  the  angle 
of  incidence  of  the  beam  Is  above  the  first  critical  angle  (-  13°  for  a  water/ 
aluminum  Interface),  only  a  vertically  polarized  shear  wave  Is  transmitted 
Into  the  sample.  Finally,  the  incident  beam  Is  backscattered  to  the  trans¬ 
ducer  only  when  the  reflector  is  nearly  at  right  angles  to  the  Incident-wave 
vector  k.  Therefore,  the  plane  of  polarization  of  the  sound  beam  Is  nearly 
perpendicular  to  the  centripetal-acceleration  vector  during  the  backseat terlng 
process.  For  this  reason,  it  can  be  assumed  that  any  Doppler  shift  measured 
as  a  result  of  backscatterlng  is  independent  of  this  acceleration.  The  magni¬ 
tude  of  the  Doppler  shift  of  a  backscattered  tone  burst  from  a  rotating  medium 
can  be  derived  by  performing  a  simple  coordinate  transformation.  In  Fig.  1(b) 
the  position  vector  r*  at  time  t,  in  the  rotating  coordinate  frame,  is  repre¬ 
sented  by 


r*  -  R  cos  (0  -  ftOx^  +  R  sin  (0  -  Ot)x2 


(1) 
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where  R  Is  the  magnitude  of  vector  jr*  and  x^  and  Xj  are  unit  vecotra  along  the 
end  Xj  **is»  respectively.  In  the  rotating  coordinate  system,  the  soli¬ 
tude  Is  given  by 


A(r*,u’)  -  A  exp  j(tc  *r*  -  w't)  (2) 

o 

where  j  -  and  u>'  is  the  circular  frequency  of  the  sound  wave  in  the  rotat¬ 
ing  frame.  The  transformation  to  the  laboratory  coordinate  frame  is  given  by 


where  the  dot  refers  to  a  time  derivative.  To  the  observer  in  the  laboratory, 
the  wave  is  described  by 


A(r,u)  -  AQ  exp  j  [£  •  r  -  (£  *  r*  +  a»')t]  (4) 

If  the  frequency  of  the  incident  beam  in  the  moving  system,  (o’,  is  related  to 
the  observed  frequency,  u>,  in  the  laboratory  system,  the  following  relation¬ 
ship  is  obtained: 

e 

us  *  a) 1  +  k  •  r* 

If  i  is  the  angle  of  incidence  in  the  x^  x3  plane,  the  wave  vector  k  can  be 
written  as 


£  -  |£|  sin  i  +  x3  |£|  cos  i 
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Using  Eqs.  (1),  (5),  and  (6),  the  expression  for  the  Doppler-shifted  frequency 
can  be  written  as 

u  ■  u'  -  |£|  R ft  sin  1  8 In  (©  -  nt)  (7) 

where  |tc|  -  u/c  and  c  Is  the  sound  speed  in  the  Incident  nedlum.  A  Doppler 
shift  occurs  when  the  wave  enters  the  rotating  media  and  again  when  the  wave 
backscattered  from  the  crack  exits  the  disk.  The  wave  vectors  of  the  Incident 
and  backscattered  waves  are  £  and  -tc,  respectively.  Therefore,  the  Doppler- 
shifted  angular  frequency,  w",  of  the  backscattered  wave  in  the  laboratory 
reference  frame  is  given  by 

w"  ■  id'  +  |£|  Rfl  sin  i  sin  (0  -  nt)  (8) 

Hence,  the  Doppler-shifted  angular  frequency  in  the  laboratory  reference  frame 
is  given  by 

Up  ■  w"  -  oi  ■  2  |tt|  R  n  sin  i  sin  (0  -  nt)  (9) 

The  net  change  in  frequency  of  an  incident  wave  after  backscattering  by  a  point 
reflector  at  radial  distance  R  in  a  rotating  medium  is  given  by 

fD  "  fo  8ln  1  8ln  (0  *  flt)  (10) 

where  T  is  the  period  of  rotation  of  the  rotating  reference  frame,  and  is  defined 
by  T  -  2ir/n,  and  is  the  frequency  of  the  incident  sound  wave.  Since  back- 
scattering  occurs  only  in  a  small  angular  range  where  0  -  tit  *  90*,  Eq.  (10)  may 
be  simplified  to 


4irR  , 

cT  o 


sin  1 


(ID 


A  Doppler- frequency  shift  f^  will  occur  for  every  point  at  radial  distance  R 
along  a  crack  or  defect  for  which  acoustic  energy  is  backscattered  from  the 
rotating  coordinate  system.  The  maximum  value  of  this  shift  should,  therefore, 
yield  the  maximum  radial  dimension  of  the  defect,  provided  the  beam  Is  back- 
scattered  from  the  tip  of  the  crack. 


III.  EXPERIMENTAL  PROCEDURE 

In  these  experiments  21  aluminum  (6061-T6)  disks  6.35  mm  thick  and  8.9  cm  In 
dlam.  having  6. 35-nm-dlam.  central  holes  were  used.  Six  disks  contained  radial 
electro-discharge-machined  (EDM)  slots  extending  radially  from  the  central  hole. 
The  slots  were  from  1.02  to  2.84  mm  in  radial  length  and  their  shape  approxi¬ 
mated  a  quarter  circle,  as  shown  schematically  In  Fig.  1.  One  disk  contained 
two  radial  EDM  slots,  180°  apart.  Four  disks  contained  three  square  EDM  slots, 
each  of  a  different  size  and  placed  at  different  radial  distances  from  the 
center  of  the  specimen.  Nine  disks  contained  radial  fatigue  cracks.  These 
cracks  were  obtained  by  placing  an  EDM  slot  in  the  center  of  each  plate,  cycli¬ 
cally  loading  the  plate  in  three-point  bending  to  grow  a  fatigue  crack  in  the 
slot,  and  then  drilling  a  hole  in  the  plate  such  that  the  crack  alone  extended 
from  the  base  of  the  hole.  The  fatigue  cracks  were  all  approximately  quarter- 
circle  shape.  Figure  2  Is  a  micrograph  of'  such  a  crack. 

To  test  the  viability  of  this  technique  in  detecting  flaws  in  the  second  layer 
of  multilayered  structures,  several  samples  were  constructed  which  consisted 
of  two  aluminum  plates  bonded  together  with  a  silicon  (RTV) ^  rubber  that 
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simulated  the  sealant  used  In  aircraft  construction.  As  shown  In  Fig.  3, 
the  second  aluminum  contained  the  crack  or  EDM  slot.  Care  was  taken  during 
bonding  of  the  two  plates  to  ensure  that  the  RTV  did  not  enter  the  slot. 

The  silicon- rubber  Inner  layer  was  nominally  0.30-mm  thick. 

From  the  theoretical  development  it  Is  apparent  that  either  the  transducer 
or  the  sample  containing  the  crack  may  be  rotated.  For  mechanical  convenience, 
the  experiments  were  conducted  using  a  rotating  sample  with  a  stationary  trans¬ 
ducer  as  shown  In  Figs.  1  and  3. 

Each  sample  was  mounted  on  a  platform  which  was  rotatable  at  a  constant  angular 
velocity  (2.  The  platform  was  placed  in  a  water  bath  which  served  as  a  couplant. 
A  focused  transducer  was  used  to  transmit  a  10-MHz  tone  burst,  typically  35  iiaec 
in  duration,  at  a  repetition  rate  of  10  KHz.  This  ultrasonic  beam  was  Incident 
at  an  angle,  1,  which  ranged  from  14  to  28°  with  respect  to  the  normal  to  the 
disk.  At  this  angle  of  incidence,  the  longitudinal  wave  mode  converted  to 
a  shear  wave  at  the  water/aluminum  interface.  The  position  of  the  transducer 
relative  to  the  hole  was  adjusted  in  such  a  way  that  the  refracted  shear  wave 
propagated  tangentially  to  the  circumference  of  the  hole.  Care  was  taken  to 
ensure  that  the  beam  did  not  intersect  the  hole  for  small  cracks  and  that  it 
covered  the  tip  of  the  crack  for  large  cracks. 

The  instrumentation  used  in  these  experiments  is  shown  schematically  in  Fig.  4. 
The  received  signal  containing  the  Doppler  shift  was  heterodyned  with  the  trans¬ 
mitted  signal  in  a  double-balanced  mixer.  The  sum-frequency  band  from  the 
mixer  was  rejected  by  a  low-pass  filter  (1-100  Hz) ,  while  the  difference  band 
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was  fed  Into  a  Hewlett-Packard  HP  3580A  low-frequency  spectrum  analyzer  and 
recorded  on  an  X-Y  plotter.  A  100-sec/dlv.  window  having  a  5-Hz  frequency 
span/dlvlslon  at  a  bandwidth  resolution  of  1  Hz  was  used.  This  sampling  rate 
was  adequate  to  avoid  aliasing  of  signals  band-limited  to  frequencies  less 
than  50  Hz.  Data  were  also  taken  with  an  HP  3582A  real-time  dual-channel 
spectrum  analyzer  which  has  better  bandwidth  resolution.  Both  analyzers 
were  used  in  the  free-running  mode  with  .the  RMS  average  of  the  spectrum 
taken  on  the  HP  3582A.  Since  the  backscattered  signal  from  the  crack  is 
transient  (occurs  only  once  during  the  period  of  rotation  of  the  disk) , 
a  uniform  filter  was  used  on  the  HP  3582A  to  achieve  optimum  results. 


IV.  RESULTS  AHD  DISCUSSION 

The  Doppler  shift  in  the  backscattered  signal  from  the  crack  is  linearly 
dependent  [Eq.  (11)]  upon  the  frequency  of  the  incident  wave,  fo>  and  the 
angular  speed,  0,  of  the  disk.  To  ensure  that  the  backscattered  signal 
was  due  to  the  Doppler  effect,  the  frequency  shift  corresponding  to  the 
maximum  amplitude  of  the  backscattered  signal  was  measured  as  a  function 
of  the  frequency  of  the  incident  wave  and  the  angular  speed  of  the  disk. 

In  both  cases  a  linear  relationship  resulted,  as  shown  in  Figs.  5(a)  and 

(b). 

The  scan  of  a  disk  without  a  notch  is  shown  in  Fig.  6.  The  spectrum  is 
flat  except  'for  the  slight  hump  around  30  Hz  caused  by  a  slight  eccentric¬ 
ity  in  alignment  of  the  transducer  and  the  central  hole  of  the  aluminum 
disk  which  results  in  scattering  from  the  hole  for  a  small  sector  in  the 
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revolution  of  the  disk.  Since  this  signal  also  has  a  Doppler  shift,  it  nay 
be  Interpreted  as  resulting  from  a  crack,  thus  illustrating  the  need  for 
careful  alignment  of  the  transducer,  especially  for  small  cracks. 

A  typical  output  for  the  Doppler  experiment  is  shown  in  Fig.  7.  It  can  be 
seen  that  the  rotating  crack,  with  its  radial  distribution  of  velocities, 
does  Indeed  produce  a  broadened  Doppler  peak.  The  spikes  in  the  envelop 
of  the  spectrum  are  due  to  the  slow  sampling  rate  of  the  HP  3580A  swept 
analyzer.  The  maximum  frequency  of  the  Doppler  shift  (the  point  at  which 
the  peak  enters  the  noise  level)  is  29.0  Hz.  If  Eq.  (11)  and  the  experi¬ 
mental  parameters  are  used,  this  frequency  corresponds  to  a  reflector 
having  a  maximum  radial  dimension  of  4.53  mm.  With  an  a  priori  knowledge 
of  the  hole  size  (3.18  mm),  the  maximum  slot  dimension  is  calculated  to 
be  1.35  mm.  The  actual' size  (1.52  mm)  was  determined  by  optical  micros¬ 
copy.  Figure  8  shows  the  Doppler  spectrum  from  a  disk  containing  a  0.94-mm 
EDH  slot  taken  with  a  HP  3582A  operated  in  the  RMS  mode.  The  maximum 
Doppler-shift  frequency  for  this  spectrum  occurs  at  43  Hz,  yielding  a 
maximum  radial  dimension  for  the  slot  of  4.16  mm.  A  calculated  radial 
size  o£  0.98  mm  is  in  excellent  agreement  with  the  actual  size  of  0.94  mm. 

When  the  Doppler-shift  technique  is  used  to  determine  flaw-size,  it  is 
necessary  to  ensure  that  the  focal  spot  of  the  ultrasonic  beam  will  cover 
the  tip  of  the  crack.  To  accomplish  this  it  is  necessary  to  scan  the 
crack  moving  the  focal  spot  radially  outward  (along  x^)  from  the  point 
of  tangency  to  the  hole  to  the  tip  of  the  crack.  When  the  maximum  fre¬ 
quency  of  the  spectra  ceases  to  increase,  the  focal  spot  covers  the 
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maximum  extent  of  the  flaw.  Figure  9  shows  the  Doppler  spectrum  taken  at  four 

.3 

different  locations  on  a  disk  containing  a  2.84-mm  EDM  slot.  As  expected,  the 
Doppler  spectrum  attained  a  maximum  frequency  of  20  Hz.  The  amplitude  of  the 
spectrum  began  to  decrease  after  this  point  since  less  of  the  ultrasonic  energy 
In  the  focal  spot  was  reflected  from  the  flaw.  Consequently,  It  Is  concluded 
that  so  long  as  the  focal  spot  of  the  beam  Intersects  the  tip  of  the  flaw  and 
does  not  intersect  the  hole,  the  beam  need  not  be  tangential  to  the  circumfer¬ 
ence  of  the  hole. 

When  a  real-time  analyzer  was  used.  It  was  possible  to  detect  multiple  radial 
cracks  in  the  time  domain.  However,  with  a  swept  analyzer,  the  size  of  the 
multiple  radial  slots  could  be  detected  and  measured  simultaneously  at  dif¬ 
ferent  angular  positions  because  the  Doppler  spectra  of  the  flaws  are  super¬ 
imposed.  From  a  maximum  Doppler  shift  from  each  peak,  the  maximum  radial 
dimension  for  each  slot  could  be  calculated.  With  a  real-time  spectrum 
analyzer,  the  spectrum  for  each  flaw  could  be  obtained  by  time  gating. 

Figure  10  shows  the  Doppler  spectrum  obtained  from  a  sample  containing  a 
fatigue  crack.  The  spectrum  Is  complex  since  fatigue  cracks  are  diffuse 
scatterers.  When  the  procedure  outlined  above  was  used,  the  crack  size  was 
estimated  to  be  1.39  mm.  Microscopic  examination  revealed  the  radial  dimen¬ 
sion  to  be  1.40  mm.  The  data  collected  on  all  disks  are  shown  in  Fig.  11. 

It  should  be  noted  that  the  data  for  each  flaw  size  were  taken  from  several 
disks  having  identical  size  slots  or  cracks,  with  at  least  three  replicates 
of  the  spectra  being  taken  on  each  disk.  The  agreement  between  the  acoustic 
and  optical  size  estimates  Is  quite  good. 
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pata  were  also  taken  on  the  three- layered  samples  containing  radial  EDM  slots 
at  the  upper  surface  of  the  lower  disk.  One  difficulty  in  taking  data  on  these 
samples  was  that  the  frequency  of  the  Incident  wave  had  to  be  such  that  no 
standing  waves  were  generated  in  the  RTV  inner  layer.  Also,  because  of  the 
high  attenuation  of  the  incident  signal  in  the  RTV  layer,  a  power  amplifier 
was  inserted  into  the  output  circuit  to  achieve  the  desired  signal- to-nolse 
ratio.  Figure  12  gives  the  experimental  results  from  these  measurements. 

Again  the  agreement  between  the  optical  and  acoustic  results  is  quite  good. 

Results  in  Figs.  11  and  12  show  that  it  is  possible  to  measure  the  radial  size 
of  cracks  under  the  holes  in  single  and  multilayered  structures.  Measurement 
accuracy  is  estimated  to  be  within  20Z.  It  should  be  mentioned  that  the  most 
frequent  error  encountered  was  failure  to  achieve  the  required  tangency  to  the 
central  hole  for  .email  flaws.  In  one  case,  the  measurement  error  was  caused 
by  a  slight  error  in  the  position  of  the  hole — its  center  was  not  at  the  rota¬ 
tional  center  of  the  plate.  Since  the  reflection  from  an  eccentrically  rotat¬ 
ing  hole  is  identical  to  that  of  a  crack,  it  may  appear  that  a  crack  is  present 
when,  in  fact,  it  is  not. 

Several  curved  fatigue  cracks  were  also  studied.  As  expected  for  the  nonradial 
cracks,  only  the  maximum  radial  dimension  of  the  crack  could  be  measured. 

V.  CONCLUSIONS 

It  has  been  demonstrated  that  the  radial  size  of  cracks  which  extend  radially  from 
bolt  holes  in  both  single  and  multilayer  mechanically  fastened  structures  can  be 
measured  using  the  Doppler-shlf t  technique.  The  technique  is  free  from  the  many 
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uncertainties  associated  with  else  measurement  by  amplitude  and  timing  methods. 
However,  the  method  Is  susceptible  to  external  vibration  and  also  requires 
adjustment  of  the  direction  of  the  incident  beam  to  avoid  its  Intersecting 
the  hole  during  rotation.  In  practice,  this  method  could  provide  a  reliable 
estimate  of  the  maximum  radial  extent  of  cracks  which  emlnate  from  bolt  holes 
in  aircraft  structures. 
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FIGURE  CAPTIONS 


(a)  Schematic  representation  of  rotating  aluminum  disk  in  water 
bath,  k  is  wave  vector  of  incident  longitudinal  wave. 

(b)  Experimental  arrangement  of  ultrasonic  beam  and  rotating 
crack.  Angle  of  incidence,  1,  is  chosen  such  that  only 
shear  wave  propagates  in  aluminum  disk. 

Optical  micrograph  of  a  typical  fatigue  crack  used  to  obtain 
the  Doppler  data.  The  bolt  hole  is  in  the  lower  left-hand 
corner  of  the  micrograph. 

Experimental  arrangement  for  Doppler  experiments  on  layered 
samples. 

Block  diagram  of  electronic  set-up  for  Doppler-shlft  studies 
of  radial  cracks. 

(a)  Doppler-frequency  shift  as  a  function  of  frequency  of 
incident  beam.  Note  linear  dependence  of  Doppler- 
shift  peak  upon  frequency  of  incident  beam. 

(b)  Doppler-shift  frequency  f^  vs  angular  speed  of  crack. 

Note  linear  dependence. 

Doppler-shlft  spectra  taken  on  disk  containing  no  flaws 
or  slot.  Low-amplitude  hump  around  30  Hz  is  due  to  a  slight 
eccentricity  of  beam  with  respect  to  the  hole.  fQ  ■  11.65  MHz, 
T  ■  A. 5  sec/rev,  i  ■  23°,  r  •  3.18  on. 
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Doppler-shlft  spectra  for  radial  EDM  notch.  Low-amplitude 

hump  under  main  peak  is  caused  by  eccentricity  in  beam 

alignment.  fQ  ■  11.65  MHz,  T  ■  4.515  sec/rev,  i  ■  17°, 

|f  \  ■  29  Hz,  c  ■  1.48  *  10^  cm/sec,  r  ■  3.18  mm, 

\  7  max 
a  -  1.35  mm. 

Doppler-shlft  spectra  for  0.94-mm  radial  slot  taken  with  real¬ 
time  analyzer.  fQ  ■  12.76  MHz,  T  “  3.24  sec/rev,  i  ■  18*, 

^f^  -  43  Hz,  r  -  3.18  mm,  a  -  0.98  mm. 

Doppler-shlft  spectra  taken  on  disk  with  slot  greater  than 
ultrasound  focal  spot  size.  Note  that  maximum  Doppler 
frequency  (20  Hz)  does  not  change  as  beam  is  moved  across 
rotating  slot  but  the  amplitude  of  Doppler  peak  drops  once 
the  focal  spot  does  not  cover  the  slot.  Spectra  D  is  for 

the  case  where  beam  does  not  Intersect  slot,  f  >12.76  MHz, 

o 

T  ■  10.0  sec/rev,  i  *  18°,  ^f^  *  20  Hz,  r  ■  3.18  mm, 

a  "  2.79  mm. 

Doppler-shift  spectra  for  real  fatigue  crack,  f  ■  12.76  MHz, 
T  ■’ 3.40  sec/rev,  i  ■  18°,  ^  •  45  Hz,  r  -  3.18  on, 

a  ■  1.39  mm. 

Plot  of  crack  size  as  measured  in  Doppler  studies  as  function 
of  actual  crack  size  for  single- layer  samples.  The  following 
notation  is  used:  o  Round  EDM  slot 
*  Square  EDM  slot 
+  Fatigue  crack 

Plot  of  crack  size  as  measured  in  Doppler  studies  as  function 
of  actual  crack  size  for  multi-layered  samples. 
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DETECTION  AND  SIZE  ESTIMATION  OF  RADIAL 
BOLT-HOLE  CRACKS  IN  LAYERED  STRUCTURES 

by 

N.  K.  Batra*,  R.  L.  Crane+,  and  P.  A.  Benadum* 


A  major  problem  In  aircraft  maintenance  Is  the  detection  and  measurement  of  ra¬ 
dial  cracks  emanating  from  bolt  holes  in  layered  structure.  Frequently,  these  cracks 
are  in  the  second  or  hidden  l&yer  of  such  structures.  While  they  can  be  detected  using 
conventional  pulse-echo  techniques,  the  measurement  of  flaw  size  using  these  methods  is 
difficult  and  subject  to  a  number  of  inherent  errors.  This  paper  briefly  describes  a 
new  technique  for  size  measurement  which  makes  use  of  the  Doppler  effect  to  relate  the 
change  in  reflected  signal  frequency  to  flaw  size  (1). 

t 

EXPERIMENTAL  PROCEDURE 

Figure  1  shows  schematically  the  physical  arrangement  used  in  the  experiments. 

The  sample  consisted  of  two  aluminum  plates  bonded  together  with  a  silicon  (RTV)**  rub¬ 
ber  that  simulated  the  sealant  used  in  aircraft  construction.  The  aluminum  plates  were 
0.635-cm  (0.250-in.)  thick  and  contained  a  0.318-cm  (0. 125-in.)  radius  central  hole.  A 
radial  electro-discharge-machined  (EDM)  slot,  which  simulated  a  fatigue  crack,  was 
placed  in  the  second  or  hidden  aluminum  plate.  Care  was  taken  during  bonding  of  the 
two  plates  to  ensure  that  the  RTV  rubber  would  not  enter  the  slot.  As  seen  in  Figure  1 
the  slot  had  the  shape  of  a  quarter  circle.  The  silicon-rubber  innerlayer  was  nominally 
0.30-mm  (0.012-in.)  thick.  This  sandwich  sample  was  mounted  on  a  rotatable  platen  and 
immersed  in  a  water  bath.  A  0.635-cm  (0.250-in.)  diameter  focused  transducer  was  used 
to  send  a  3.5  ps  long  tone  burst  of  7  MHz  acoustic  energy  into  the  sample  at  a  repeti¬ 
tion  rate  of  10  KHz.  The  angle  of  incidence  was  held  constant  at  18  degrees.  Since 
this  angle  is  above  the  first  critical  angle  of  the  water /aluminum  interface,  a  verti¬ 
cally  polarized  shear  wave  was  generated  in  the  sample.  This  combination  of  angle  of 
incidence  and  frequency  was  chosen  to  provide  the  largest  signal  return  (reflection) 
from  the  simulated  crack  in  the  specimen.  In  most  cases  the  ultrasonic  beam  was  posi¬ 
tioned  such  that  its  focal  spot  was  tangential  to  the  hole  in  the  second  layer.  In  the 
case  where  the  focal  spot  of  the  transducer  was  smaller  than  the  flaw  size,  the  ultra¬ 
sonic  beam  was  positioned  in  such  a  way  that  focal  spot  covered  the  crack  tip.  The  re¬ 
flected  signal  was  received  by  the  transducer,  mixed  with  the  rf  carrier  and  filtered, 
and  the  output  was  displayed  on  a  Hewlett-Packard  3582  spectrum  analyzer.  A  schematic 
diagram  of  the  electronic  set-up  is  shown  in  Figure  2. 
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THEORY 


The  Doppler  effect  typically  refers  to  the  change  In  frequency  which  occurs  when  a 
sound  berm  Is  reflected  from  a  discontinuity  that  Is  moving  in  the  medium  of  an  Inci¬ 
dent  wave.  A  Doppler  shift  can  also  occur  when  the  sound  wave  enters  a  medium  that  is 
in  motion  relative  to  the  emitter.  This  latter  case  Is  applicable  to  the  experiment 
described  schematically  in  Figure  1.  Here  the  ultrasonic  wave  experiences  a  frequency 
shift  at  the  water/aluminum  interface  both  upon  entering  and  exiting  the  sample.  This 
Doppler  shift  has  a  range  of  frequencies  as  a  result  of  the  radial  distribution  of  velo¬ 
cities  across  the  beam  spot  on  the  upper  surface  of  the  sample.  As  a  result,  the  ex¬ 
tent  of  reflection  of  the  signal  by  the  crack  determines  the  extent  of  the  Doppler 
shift.  It  should  be  noted  that  the  reflection  of  the  signal  occurs  over  a  very  small 
angular  range  because  of  the  physical  arrangement  of  the  experiment  and  the  small  focal 
spot  size  of  the  transducer.  With  these  considerations  in  mind,  it  is  possible  to  de¬ 
rive  a  mathematical  relationship  relating  the  Doppler  shift  to  the  crack  dimensions  (1). 
The  frequency  shift  of  an  ultrasonic  wave  that  enters  a  rotating  media  at  radial  dis¬ 
tance  R  and  is  reflected  back  into  a  non-moving  media  is  given  by, 

f  „  l»lo_R  8ln  t  (1) 

D  CT 

where  fj)  is  the  Doppler  frequency,  f0  is  the  incident  frequency,  c  is  the  wave  speed  of 
sound  in  water,  t  is  the  period  of  rotation  of  the  sample,  and  i  is  the  angle  of  inci¬ 
dence  of  the  ultrasonic  beam  in  water.  This  equation  can  easily  be  used  to  establish 
the  maximum  extent  of  the  crack.  With  a  priori  knowledge  of  the  hole  6ize,  the  crack 
dimension  can  be  determined. 


RESULTS  AND  DISCUSSION 

The  data  from  a  typical  experiment  is  shown  in  Figure  3.  The  time  record  of  the 
wave  reflected  from  the  EDM  slot  is  shown  in  Fig.  3A,  while  the  spectrum  of  the  signal 
is  shown  in  Fig.  3B.  The  bandwidth  of  the  Doppler  spectrum  extends  from  approximately 
16  to  35  Hz.  Using  Equation  (1)  and  the  data  shown  in  Figure  3,  this  corresponds  to 
radial  reflector  dimensions  of  2.6  mm  (0.1  in.)  to  5.9  mm  (0.23  in.).  Assuming  the  ul¬ 
trasonic  focal  spot  completely  covers  the  slot,  then  these  data  indicate  a  radial  EDM 
slot  dimension  of  3.3  mm  (0.13  in.).  The  actual  dimensions  of  3.175  mm(0.125  in.)  for 
the  radius  of  the  hole  and  3.0  nun  (0.12  in.)  for  the  slot  correspond  favorably  to  the 
predicted  sizes.  Note  that  it  is  the  maximum  frequency  of  the  Doppler  shift  or  the 
maximum  radial  slot  dimension  that  is  important  rather  than  the  absolute  distances.  For 
actual  measurements  the  hole  size  is  known  a  priori,  eliminating  the  uncertainty  in  the 
measurement  of  one  parameter.  Figure  4 A  shows  the  time  record  of  data  taken  on  a  sample 
containing  two  radial  EDM  slots  180°  apart,  one  being  1.24  mm  and  the  other  1.32  mn. 
Figure  4B  is  an  expanded  region  of  the  time  record  which  shows  that  a  slot  is  in  view 
for  only  0.185  seconds.  Figure  4C  shows  the  Doppler  spectrum  from  one  of  the  slots.  It 
was  possible  using  the  HP3582  spectrum  analyzer  to  detect  multiple  radial  crack*  using 
the  time  mode;  here,  however,  it  was  not  possible  to  estimate  the  size  of  a  slot  unique¬ 
ly  when  two  or  more  such  slots  were  present.  However,  it  would  be  possible  to  measure 
the  size  of  one  slot  uniquely  if  the  motion  of  the  platform  could  be  synced  with  the 
spectrum  analyzer. 

Table  I  contains  a  summary  of  the  data  taken  thus  far  for  sandwich  samples  con¬ 
taining  both  EDM  slots  and  fatigue  cracks.  It  is  apparent  that  the  Doppler  technique 
is  quite  capable  of  yielding  a  reasonable  estimate  of  flaw  size. in  layered  or  sandwich 
aircraft  structure.  An  experimental  bound  on  the  error  in  these  measurements  has  been 
found  to  be  about  20  percent.  While  it  is  possible  to  achieve  much  better  accuracy  if 
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TABLE  I.  Summary  of  Experimental  Results 


Sample  Measured  Crack 

No.  Size  (Optical) 


Estimated  Crack 
Size  (Doppler) 


mm 

in. 

nm 

in. 

1 

1.32 

0.052 

1.65 

0.065 

2 

1.52 

0.060 

1.57 

0.062 

3 

1.79 

0.070 

1.93 

0.076 

4 

2.03 

0.080 

1.91 

0.075 

5 

2.84 

0.112 

2.72 

0.107 

great  care  is  taken,  such  effort  does  not  seem  to  be  warranted  at  this  time.  It  should 
be  mentioned  that  the  most  frequent  error  encountered  was  the  failure  to  achieve  the 
required  beam  tangency  to  the  central  hole.  In  one  case  the  measurement  error  was 
caused  by  a  slight  error  in  the  position  of  the  hole,  such  that  its  center  was  not  at 
the  rotational  center  of  the  plate.  Since  the  reflection  from  an  eccentrically  rotat¬ 
ing  hole  is  quite  similar  to  that  of  a  crack,  it  may  appear  that  a  crack  is  present 
when,  in  fact,  it  is  not.  Also,  with  the  existing  experimental  set-up,  alignment  of 
the  ultrasonic  beam  with  the  central  hole  is  quite  difficult.  Again  a  reflection  from 
the  hole  can  mimic  a  crack. 


Another  possible  difficulty — one  that  plagues  all  ultrasonic  inspections  of  layere 
aircraft  structure — is  the  total  reflection  of  ultrasonic  energy  from  a  layer  that  is 
(n+l)/4  wavelengths  thick.  This  problem  is  well  known  to  those  using  wide-band  ultra¬ 


sonic  techniques.  In  the  Doppler  experiment,  the  problem  can  become  acute  because  of 


the  use  of  a  tone  burst  or  single-frequency  ultrasonic  energy.  When  this  problem  was 


encountered  in  the  present  study,  the  frequency  of  the  incident  wave  was  changed. 


CONCLUSIONS 

It  has  been  shown  that  the  Doppler  effect  may  be  used  to  estimate  the  size  of  ra¬ 
dial  crack-like  flaws  in  multilayer  aircraft  structures.  The  technique  is  capable  of 
predicting  crack  size  with  an  accuracy  of  better  than  ?0  percent.  While  this  technique 
may  not  be  used  to  detect  cracks  emanating  from  bolt  holes,  it  does  represent  a  possibl 
method  of  estimating  defect  size.  Also  it  does  not  suffer  from  many  of  the  inherent 
errors  that  limit  amplitude  methods  for  measuring  crack  size. 


REFERENCES 

1.  Batra,  N.  K.  and  Crane,  R.  L.  ,  "Measurement  of  Radial  Bolt  Hole  Crack  Size  Using 
Doppler-Shif t  Techniques",  to  be  published. 


8U 


iW 

i* 


TRANSDUCER 


FIGURE  1 — Experimental  arrangement  of  Doppler 
experiments. 
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FIGURE  2 — Block  diagram  of  electronic  set-up 
used  in  Doppler  experiments. 
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FIGURE  3A — Tine  record  output 

for  e  typical  Doppler 
experiment. 


FIGURE  3B — Fast-Fourler  analysis 
of  time  record  of 
Fig.  3 A.  Spike  at 
45  Hz  Is  reference 
marker. 
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FIGURE  AA — Time  record  output  for 
sandwich  specimen  con¬ 
taining  two  EDM  slots. 


cHhNNEL  h 

LOM-PflS§  F  ILTEREUi-T  I  MET  FUNCT  I  OHj>' 

•  ' *  t>s  '  •:  ..-A  .•».  *  = M 


..  s  m  ' 

~V*>  ,  ,  J  v- 

•  csV"  . 

..  ,  -w  s  • 

'.  A  -  v  ’  s  ■/*»?»•  .4*'  .  '  t  . 


.  1/  :•* 

;  *•  V.  ' '  ^v*. Xv.  >*  j 

I**..  ■  •  ►  .  ;  ■*  '•*'•**&  ■  -’-v 

W*  «£ 


.  "  ,4  •  V!$r  -  ■ 

HLTERNflTE  SflNPLES  'flFTER  FILTER-' 


FIGURE  AB — Tima  expanded  version 
of  one  of  the  echoes 
In  Fig.  AA. 


FIGURE  AC — Spectral  analysis  of 
one  of  the  slots 
indicated  in  Fig.  A  A. 
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An  Ultrasonic  Pseudorandom-Signal-Correlation  System 

CHARLES  M.  ELIAS,  Student  Member  IEEE 

ABSTRACT 

A  working  ultrasonic  pseudorandom  signal-correlation  system  is  described  which, 
unlike  ultrasonic  random  signal-correlation  systems,  does  not  require  an  acous¬ 
tic  delay  line.  Elimination  of  the  delay  line  allows  faster  data  acquisition 
and  better  range  resolution.  The  system  uses  two  identical  shift-register- 
type  generators  to  produce  pseudo-noise  bursts  which  are  subsequences  of  a 
65  535-bit  complementary  m-sequence.  One  generator  produces  the  transmitted 
bursts  while  the  other  generates  identical  reference  bursts  which  start  at  a 
variable  correlation  delay  time  after  the  transmitted  bursts.  The  reference 
bursts  are  cross-correlated  with  the  received  echoes  to  obtain  the  approximate 
Impulse  response  of  the  transducer /specimen  system  under  test.  Range  side- 
lobes  are  reduced  by  transmitting  and  correlating  many  bursts  at  a  given  corre¬ 
lation  delay  before  incrementing  the  delay.  Signal-to-sidelobe  ratios  of 
greater  than  47  dB  have  been  obtained  using  this  method. 

Limitations  of  the  system  due  to  sampling  constraints  and  the  pseudo-noise 
power  spectrum  are  discussed,  and  the  system  design  and  implementation  are 
outlined.  Results  of  experimental  characterization  of  the  system  show  that 
the  pseudorandom  signal-correlation  system  has  approximately  the  same  range 
resolution  as  a  conventional  pulse-echo  system  but  can  yield  a  significant 
Increase  in  slgnal-to-nolse  ratio. 
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An  Ultrasonic  Pseudorandom-Signal-Correlation  System 

CHARLES  M.  ELIAS,  Student  Member  IEEE 

INTRODUCTION 

Random-signal  correlation  techniques  are  well  known  In  the  radar-systems  field 
as  pulse-compression  methods  [1].  In  a  pulse-compression  system,  the  pulse 
transmitted  has  a  wide  bandwidth  (for  hlgh-range  resolution)  and  long  duration 
(for  high  velocity  resolution  and  high  transmitted  energy) .  The  received 
echoes  are  processed  to  obtain  short  pulses  [2]. 

In  the  ultrasonics  area  a  pulse-compression  system  which  makes  use  of  random 
signal-correlation  techniques  for  flaw  detection  has  been  described  by 
Furgason,  et  al.  [3]  and  evaluated  by  Bilgutay,  et  al.  [4].  This  system  was 
shown  to  provide  a  much  improved  signal- to-noise  ratio  and  a  much  lower  peak- 
to-average  power  ratio  as  compared  to  conventional  pulse-echo  methods.  In 
the  correlation  method  used,  random  noise  bursts  are  transmitted,  therefore, 
a  "reference"  noise  burst  that  is  a  replicn  of  the  transmitted  burst,  but  is 
delayed  t  seconds  relative  to  the  transmitted  burst.  Is  required,  t  must  be 

This  research  is  sponsored  by  the  United  States  Air  Force  Materials  Laboratory 
under  Contract  F33615-77-C-5022. 

The  author  is  with  Systems  Research  Laboratories,  Inc.,  2800  Indian  Ripple 
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variable  over  the  time  Interval  In  which  echoea  of  Intereat  will  be  received. 
That  la,  the  correlation  process  as  performed  by  the  system  Is  a  finite  Inter¬ 
val  approximation  of  the  cross-correlation  function  of  the  transialtted  burst 
f(t)  and  the  received  burst  g(t).  This  approximation,  designated  C(t),  Is 
mathematically  given  by 

.  T  t  t/2 

Z(X)  -  i  Jo  f(t-T)g(t)dt  i  £  J  f(t-T)g(t)dt 

The  reference  bursts  are  analogous  to  f(t-T>  in  (1),  and  the  averaging  and 
Integration  are  simulated  with  a  low-pass  filter.  However,  since  each  trans¬ 
mitted  burst  is,  by  definition,  unique,  the  signal  from  a  noise  source  Is 
split  into  two  channels,  one  channel  producing  the  transmitted  bursts  and 
the  other,  which  contains  a  variable  delay  line,  producing  the  reference  bursts. 
The  variable  delay  line  consists  of  two  ultrasonic  transducers  in  a  water 
bath,  and  the  delay  is  varied  by  changing  the  separation  between  the  trans¬ 
ducers.  Such  a  delay  line  is  bulky  and  has  a  slow  mechanical  scan.  Further¬ 
more,  the  presence  of  the  two  transducers  in  the  reference  signal  path 

significantly  decreases  the  path  bandwidth  and  degrades  the  range  resolu- 

* 

don  of  the  system  [5]. 

The  present  paper  describes  an  ultrasonic  correlation  system  in  which 
pseudorandom  binary  noise  generators  produce  the  transmitted  and  reference 
noise  bursts  and,  therefore,  no  delay  line  Is  required.  Results  of  experi¬ 
mental  characterization  of  the  system  are  given.  Similar  work,  using 

phase-modulation  techniques,  is  being  done  by  Chapelon,  et  al.  [6]  at  the 
* 

Although  not  explicitly  mentioned,  this  degradation  is  implicit  In  Eq.  (34) 
of  this  reference. 


(1) 
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Instltut  National  de  la  Sante  et  de  la  Recherche  Medlcale  In  Bran,  France. 

At  the  time  of  this  writing,  however,  inaufficlent  detailed  information  is 
available  for  a  comparison  with  the  preaent  work.  Pseudorandom  binary  coding 
has  also  been  used  by  Waag,  et  al. ,  in  an  ultrasonic  Doppler  system  for  non- 
invaslve  cardiac  chamber  flow  measurement  [7]. 

THEORY 

Pulse  compression  is  obtained  in  the  ultrasonic  pseudorandom  signal-correlation 
(UPSC)  system  by  deriving  the  Impulse  response  of  the  "system",  which  is  com¬ 
posed  of  the  ultrasonic  transducer (s)  and  the  specimen,  from  a  cross  correla¬ 
tion  of  the  received  echo  and  the  reference  bursts.  (This  assumes  that  the 
transmitting  amplifier,  receiving  amplifier,  and  correlator  have  no  signifi¬ 
cant  effect  upon  the  impulse  response.) 

It  is  known  from  systems  theory  that  if  "white"  noise  is  input  to  a  linear 
system,  then  the  cross  correlation  of  the  input  noise  and  the  system  output 
will  be  proportional  to  the  system  Impulse  response.  Although  the  theory 
specifies  white  noise,  which  has  a  constant  power  density  spectrum  for  all 
frequencies  and  an  autocorrelation  function  which  is  an  impulse,  a  good 
approximation  of  system  Impulse  response  can  be  obtained  using  band-limited 
noise  bursts.  This  was  demonstrated  in  the  random  signal-correlation  system 
mentioned  earlier  [3,4].  The  UPSC  system  utilizes  pseudorandom  binary  noise 
bursts  to  provide  a  good  approximation  to  the  Impulse  response. 

L  l 


A  pseudorandom  binary  sequence  or  pseudo-noise  sequence  is  a  periodic  sequence 
of  bits  having  certain  properties  associated  with  randomness  [8].  A  l 

length  or  m-sequence  Is  a  pseudo-noise  sequence  which  can  be  obtained  as  the 
serially  docked  output  of  a  binary  linear  shift  register  when  specific  feed¬ 
back  connections  are  made  by  a  modulo-2  adder.  An  m-sequence  has  a  period  or 
length  p  ■  2n  -  1  bits,  with  n  being  the  number  of  stages  In  the  shift  regis¬ 
ter  [8,9].  The  autocorrelation  function  of  an  m-sequence  {a^}  is  defined 
»•  [8], 


P 

C(t)  -  Z 

k 


*k*k  +  t 

1 


t  -  0,±1,±2,... 


(2) 


For  {a^}  having  binary  values  + 1  and  -  1  , 

p  if  t  -  qp,  q  -  0,1,2,... 

C(t)  -  (3) 

-1  if  T  y  qp 

This  -relationship  can  be  obtained  from  the  known  properties  of  m-sequences  and 
the  result  from  Ref.  8  for  m-sequences  having  binary  values  of  0  and  1.  From 
Eq.  (3)  it  can  be  seen  that  the  autocorrelation  of  pseudo-noise  is  not  a  single 
impulse  at  r  ■  0,  as  in  the  case  of  white  noise,  but  a  triangular  pulse  which 
has  height  equal  to  p,  width  at  its  base  of  two  units,  and  period  p  units  on 
the  T  axis.  If,  however,  an  m-sequence  is  considered  as  a  burst  of  length  p 
bits  or  If  It  is  broken  up  into  short  bursts,  the  autocorrelation  function  will 
have  other  peaks  or  "range  sldelobes"  In  addleion  to  the  main  peak  [1]. 


93 


Another  important  property  of  an  m-sequence  Is  its  power  spectrum  which  has  nulls 
at  Integer  multiples  of  the  clock  frequency  f£  and  line  spacing  inversely  pro¬ 
portional  to  the  sequence  length  p  [8,10].  Thus,  by  proper  choice  of  p  and  fc, 
the  power  spectrum  of  the  pseudo-noise  can,  in  principle,  be  made  to  approximate 
that  of  white  noise  sufficiently  for  a  given  application. 

GENERAL  DESCRIPTION  OF  THE  SYSTEM 

Figure  1  is  a  block  diagram  of  the  experimental  UPSC  system  used  in  this  investi¬ 
gation.  This  system  utilizes  two  pseudorandom  binary  noise  generators — PNG  #1 
and  PNG  #2 — to  generate  the  required  pseudorandom  binary  sequences.  The  PNG's 
employ  shift  registers  having  exclusive  NOR  feedback  to  generate  the  pseudo¬ 
noise  sequence.  PNG's  #1  and  #2  use  identical  shift  registers  having  the  same 
feedback  circuitry  to  produce  an  m-sequence  of  length  p.  The  feedback  connec¬ 
tions  determine  a  unique,  periodic  m-sequence;  therefore,  if  initial  conditions 
are  the  same,  identical  PNG's  will  produce  identical  pseudo-noise  subsequences 
(or  noise  bursts)  when  enabled  for  N  (N  <  p)  clock  periods.  This  is  precisely 
how  PNG '8  #1  and  #2  generate  the  transmitted  and  reference  noise  bursts  in  the 
experimental  system.  The  timing  diagram  of  Fig.  2(a)  shows  several  transmitter 
cycles,  where  Q  is  the  time  between  successive  transmitted  bursts  or  the  period 

of  a  transmitter  cycle,  x  is  the  initial  reference  burst  delay,  and  t  is  the 

6 

time  to  the  received  burst  (assuming  a  fixed  target). 

System  operation  proceeds  as  follows:  both  PNG's  are  set  to  the  same 
state,  and  the  burst  length  is  set  to  N  clock  periods.  PNG  #1  generates  the 
initial  transmitted  burst.  After  a  delay  of  t  clock  periods,  PNG  #2  generates 


a  reference  buret,  t  +  At  clock  periods  after  PNG  #1  produces  the  next  trans¬ 
acted  burst.  PNG  #2  produces  the  second  reference  burst.  As  this  process  con¬ 
tinues.  with  the  reference  burst  delay  Increasing  In  Increments  of  At,  the 
analog  multiplier  output  Is  the  Instantaneous  product  of  the  reference  burst 
and  the  received  signal  which,  after  low-pass  filtering  and  amplification, 
becomes  the  correlator  output.  The  correlator  output  Is  Illustrated  In 
Pig.  2(b).  Note  that  the  time  between  successive  points  on  the  correlated 
output  Is  Q  +  At;  that  Is,  one  new  point  of  the  correlated  output  Is  genera¬ 
ted  for  each  transmitted  burst  (or,  more  precisely,  for  each  time  that  t  is 
Incremented  by  At) .  A  At  change  on  Fig.  2 (a)  ("real  time")  corresponds  to  a 
Q  +  At  change  In  the  correlated  output  time  frame  ("correlation  time").  Since 
Q  »  At,  then  At/Q  yields  the  ratio  of  "real  time"  to  "correlator  time".  In 
other  words,  with  the  system  operating  as  described,  if  the  correlator  output 
Is  observed  on  an  oscilloscope,  the  scope  time  base  setting  would  be  multiplied 
by  1/Q  to  obtain  "real  time".  Or,  conversely,  a  response  which  would  take  t 
seconds  In  "real  time"  would  be  traced  out  in  Q*  t  seconds  at  the  correlator 
tfUtput . 

In  the  above  description  of  the  system  operation,  one  very  Important  feature 
of  this  system  has  been  omitted  for  the  sake  of  simplicity.  For  reasons 
which  will  be  subsequently  explained,  it  Is  not  desirable  to  increment- t  by 
At  In  every  transmitter  cycle.  In  the  actual  system,  the  correlation  delay 
t  remains  unchanged  for  I  transmitter  cycles  prior  to  each  Increase  of  At. 

Thl«  Mans  that  I  subsequences  are  transmitted  and  I  reference  subsequences 
are  generated  by  PNG  #2  at  some  fixed  delay  t  prior  to  being  Increased  by  At 
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and  the  process  being  repeated.  Nov  the  acale  factor  1/Q  becomes  1/QI  since 
the  previous  description  assumed  I  to  be  unity.  The  reason  for  the  use  of 
I  "repetitions"  per  correlation  point  Is  to  reduce  range  sldelobes.  The 
"algtial-to-sldelobe  ratio"  Is  defined  as  the  ratio  of  the  amplitude  of  the 
main  correlation  peak  to  that  of  the  largest  range  sldelobe.  The  slgnal- 
to-sldelobe  ratio  might  be  Increased  In  the  following  ways: 

1.  Continuous  Operation  of  the  Transmitter  PNG  (cw  Operation) 

This  mode  Is  very  attractive  for  certain  applications  since  it  produces 
no  range  sldelobe  peaks.  Uhlle  the  maximum  range  of  a  cw  system  is 
limited  by  the  periodicity  of  the  m-sequence  autocorrelation  function 
[1],  this  limitation  can  be  overcome  by  using  a  very  long  m-sequence. 

The  most  serious  limitation  of  the  cw  mode  Is  that  its  use  virtually 
precludes  single-transducer  operation.  For  this  reason  the  cw  mode 
was  not  uaed  In  the  UPSC  system. 

2.  Use  of  a  Complete  m-Sequence  as  a  Burst 

In  this  case,  the  signal-to-sidelobe  ratio  should  be  on  the  order  of  /p 
where  p  Is  the  sequence  period  [1].  But,  for  a  system  where  many  different 
burst  lengths  are  required,  this  alternative  presents  design  complications. 
The  shift  register  feedback  logic  would  be  quite  complex  since  a  different 
feedback  configuration  Is  required  for  each  change  In  sequence  length. 

Such  complex  logic  would  not  only  tend  to  limit  the  highest  clock  frequency 
of  the  pseudo-noise  generators,  but  It  would  also  nullify  two  primary  advan¬ 
tages  of  using  shift  registers  to  generate  the  pseudo-noise  sequences: 
simplicity  and  ease  of  implementation. 
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Use  of  I  Repetitions  wr  Correlation  Point 


The  process  of  transmitting  end  correlating  I  bursts  before  changing  the 
correlation  delay  r  is  somewhat  analogous  to  the  operation  of  the  randoa 
signal  correlation  system  described  earlier  [3].  However,  to  the  best  of 
the  author's  knowledge,  this  technique  has  not  been  used  previously  in  a 
pseudo-noise  system.  Use  of  this  method  greatly  simplifies  control  of  the 
burst  length.  This  was  considered  to  be  especially  important  in  an  experi¬ 
mental  system  where  as  many  system  parameters  as  possible  should  be  easily 
controllable.  Of  course,  an  increase  in  I  increases  the  data-acquisltlon 
time,  but  the  trade-off  of  data-acquisltlon  time  for  signal- to-noise  ratio 
enhancement  is  an  inherent  property  of  this  type  of  system.  Slow  range 
scanning  (large  I)  allows  a  lower  low-pass  filter  cutoff  frequency  which 
provides  improved  signal- to-noise  ratio  and  simultaneous  improvement  in 
the  signal- to-sidelobe  ratio.  That  is,  the  ratio  of.  signal- to-noise  power 
at  the  system  output  to  that  at  the  receiving  amplifier  can  be  expressed 
as 


BV 


in 


BW 


out 


(4) 


where  BW1q  and  BW  are,  respectively,  the  bandwidth  at  the  system  input 
(usually  transducer  limited)  and  the  system  output  (set  by  the  low-pass- 
filter  cutoff  frequency)  [3]. 


It  Is  apparent  from  the  foregoing  deacrlptlon  of  the  system  operation  that  the 
CPSC  system  Is,  In  effect,  sampling  an  Impulse  response  waveform  at  i.  sampling 
frequency  of  1/At  Hz  (recall  that  a  new  sample  datum  Is  available  each  time  t 
Is  Incremented  by  At) .  This  sampling  frequency  must  be  greater  than  twice  that 
of  the  highest  spectral  component  of  significant  amplitude  which  is  present  In 
the  Impulse  response  of  the  transducer/specimen  system  under  test.  Another 
sampling  effect  which  occurs  in  the  UPSC  system  correlation  process  can  be 
explained  by  considering  the  autocorrelation  of  two  pseudo-noise  bursts  when 
these  bursts  have  some  arbitrary  phase  relationship.  Refer  to  Fig.  2(a)  and 
suppose  that  t  is  repeatedly  Incremented  by  At  and  the  product  of  the  two 
bursts  Is  obtained  each  time  t  Is  Incremented.  This  can  be  considered  a 
sampling  process  with  the  sampling  frequency  again,  being  1/At.  In  this  case, 
however,  1/At  should  be  twice  the  pseudo-noise  clock  frequency  in  order  to 
sample  twice  per  bit  time.  This  sampling  rate  will  be  adequate  in  practice 
because  the  received  bursts  are  band-limited  by  the  transducer/specimen  sys¬ 
tem  (see  Fig.  2(c)  for  an  example  of  an  actual  received  burst). 


SYSTEM  DESIGN  AND  IMPLEMENTATION 


A  simplified  block  diagram  of  the  logic  subsystem  of  the  UPSC  system  Is  shown 
In  Fig.  3.  The  subsystem  is  totally  synchronous,  and  all  timing  functions  are 
performed  by  presettable  synchronous  counters  which  are,  in  turn,  controlled 
by  the  primary  control  logic.  Storage  registers  hold  numbers  which  are  used  as 
values  for  counters  or  directly  by  the  primary  control  logic.  These  numbers 
are:  burst  length.  Initial  reference  burst  delay  (Initial  value  of  t)  maximum 
reference  burst  delay  (final  value  of  t),  number  of  repetitions  per  correla¬ 
tion  point  (I),  and  transmitter  cycle  length  (Q).  The  storage  registers  are 
selected  and  loaded  using  front-panel  thumbwheel  switches.  Each  pseudo-noise 
generator  uses  a  16-stage  shift  register  with  exclusive  NOR  feedback  from 
Stages  1,  3,  12,  and  16.  This  configuration  will  produce  a  65  535-blt 
complementary  m-sequence  (this  sequence  having  properties  nearly  identical 
to  those  of  an  m-sequence).  The  complementary  sequence  was  chosen  because 
It  permits  use  of  the  all-zero  state,  a  forbidden  state  in  the  m-sequence. 

This  means  that  initialization  of  the  generators  can  be  performed  for  most 
Integrated-circuit  shift  registers  simply  by  activating  the  clear  input  [9]. 

The  system  master  clock  period  is  the  minimum  increment  by  which  any  time 
interval  or  delay  can  be  changed.  Due  to  the  sampling  considerations  dis¬ 
cussed  earlier,  PNG's  #1  and  #2  are  run  at  one-half  the  master  clock  fre¬ 
quency. 
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The  logic  In  this  system  was  Implemented  using  (primarily)  Schottky-damped 
transistor-transistor  logic  (TIL)  Integrated  circuits  with  wire-wrapped  inter¬ 
connections.  The  highest  master  clock  frequency  obtained  with  error-free 
operation  was  25  MHz.  No  effort  has  been  made  to  Increase  system  speed  beyond 
this  point;  however,  the  logic  speeds  in  this  system  do  not  approach  state- 
of-the-art  speeds.  Uith  faster  logic  and  careful  design,  a  much  faster  system 
could  be  built.  A  typical  experimental  set-up  for  the  UPSC  system  is  shown 
in  Fig.  4.  A  Philips  PM  3265  dual-trace  multiplier  oscilloscope  is  used  as 
a  variable-gain  broad-band  amplifier  and  analog  multiplier.  The  received  sig¬ 
nal  and  the  reference  bursts  are  fed  to  the  vertical  inputs  of  the  oscilloscope. 
The  oscilloscope  is  set  in  the  multiply  mode  to  permit  the  product  of  the  two 
vertical  amplifier  outputs  to  be  available  at  the  Y  OUT  jack.  This  output  is 
fed  to  an  active  low-pass  filter  having  variable  cutoff  frequency.  The  low- 
pass-filter  output  is  the  correlated  signal  and  can  be  observed  on  a  second 
oscilloscope  as  shown. 

The  logic  subsystem  supplies  synchronization  signals  for  both  oscilloscopes. 

The  bipolar  converters  employ  burst  and  gate  inputs  to  change  the  TTL  bursts 
to  "bipolar"  bursts.  Bipolar  pseudo-noise  burst  trains  were  used  because 
they  experience  miniman  transient  effects  when  applied  to  ac  coupled  ampli¬ 
fiers. 
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EXPERIMENTAL  RESULTS 


The  result  of  a  comparison  between  the  Impulse  response  of  two  transducers  in 
a  water  bath  as  measured  with  a  conventional  pulse-echo  ultrasonic  system  and 
the  UPSC  system  is  shown  in  Fig.  5.  The  similarity  of  the  responses  is  striking. 
It  is  noteworthy  that  the  Panametrics  unit  used  a  390-V  pulse,  while  the  UPSC 
system  used  a  2-V  p-p  burst  in  obtaining  the  response  shown. 

Figure  6  is  an  example  of  slgnal-to-sidelobe  ratio  improvement  obtained  by 
increasing  the  number  of  repetitions  per  correlation  point.  Figure  7  is  a 
plot  of  the  mflvimiim  sidelobe  magnitude  (relative  to  the  correlation  peak)  as 
a  function  of  repetitions  per  correlation  point  for  a  fixed  burst  length  of 
375  bits.  The  experimental  data  were  taken  using  the  UPSC  system  in  the 
through  transmission  mode  in  a  water  bath.  Two  raw  data  sources  were  used 
to  obtain  the  experimental  results  plotted  in  Fig.  7:  an  oscilloscope  and 
a  Northern  Model  NS-575  signal  averager.  Agreement  between  the  calculated 
and  experimental  results  is  excellent  up  to  100  repetitions.  Beyond  this 
point,  the  computer  model  predicts  a  more  rapid  decrease  in  sidelobe  level 
with  Increasing  repetition  number,  while  the  experimental  values  appear  to 
be  following  the  original  slope.  At  the  present  time  it  is  not  known  if  this 
discrepancy  is  due  to  experimental  errors,  since  the  sidelobes  are  near  sys¬ 
tem  noise  levels,  or  to  an  inadequate  computer  model  of  the  system  (e.g., 
the  model  does  not  take  into  account  transducer  effects).  For  the  375-blt 
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burst  length,  sldelobe  levels  at  least  47  dB  below  the  peak  height  have  been 
obtained  experimentally.  As  stated  In  Ref.  [11],  obtaining  the  -47  dB  ratio 
required  an  acquisition  time  of  almost  2  min.  It  should  be  noted  that  the 
sldelobe  level  can  also  be  reduced  by  Increasing  the  burst  length  for  a  given 
number  of  repetitions. 

Figure  8  Is  an  example  of  the  signal- to-nolse  ratio  improvement  obtainable 
with  the  UPSC  system. 

A  simple  experiment  was  performed  to  demonstrate  that  system  resolution  is 
degraded  by  the  presence  of  a  pair  of  transducers  in  the  reference-signal 
path,  as  is  the  case  for  the  random  signal  correlation  system  of  Ref.  [3]. 
Figure  9(a)  shows  the  Impulse  response  of  two  transducers  In  a  water  bath 
with  the  UPSC  system  operating  normally.  Figure  9(b)  shows  the  result  of 
the  same  experiment  but  with  two  additional  transducers  in  a  water  bath 
being  Inserted  into  the  reference  signal  path. 

CONCLUSIONS 

This  study  has  shown  that  a  UPSC  system  can  yield  results  comparable  to 
those  obtained  with  a  random  signal  correlation  system.  The  UPSC  system 
has  been  shown  to  have  better  range  resolution  than  a  random  signal  system 
which  makes  use  of  a  narrow  band  variable  delay  line.  Nevertheless,  the 
range  resolution  of  the  UPSC  system  can  be  no  better  than  that  of  an  ideal 
pulse-echo  system  since  both  are  limited  by  the  transducer  Impulse  response. 
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Like  a  pulse-echo  system,  the  range  resolution  of  the  UPSC  system  is  improved 
by  using  higher-frequency  transducers.  For  the  UPSC  system,  however,  considera¬ 
tion  must  be  given  to  the  power  spectrum  of  the  pseudo-noise  which  is  3-db  down 
at  about  one-half  the  pseudo-noise-generator  clock  frequency  and  has  nulls  at 
integer  multiples  of  that  frequency. 

UPSC  systems  could  be  applied  to  signal  detection  in  highly  attenuative  mate*' 
rials.  Such  systems  might  also  be  used  in  the  implementation  of  those  quanti¬ 
tative  nondestructive  evaluation  techniques  which  employ  spectral  analysis  and, 
therefore,  require  a  large  signal- to-noise  ratio.  If  long  acquisition  times 
are  a  problem,  faster  logic  and  parallel  processing  using  multiple  correlators 
might  be  used. 
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FIGURE  CAPTIONS 


Figure 

1  Block  diagram  of  experimental  ultraaonlc  pseudorandom  signal 
correlation  system.  PNG  #1  and  PNG  #2  are  the  transmitter  and 
reference  pseudorandom  binary  noise  generators,  respectively. 

2  (a)  Timing  diagram  of  the  ultrasonic  pseudorandom  signal 
correlation  system,  (b)  Correlator  output,  (c)  Actual 
receiving  transducer  signal  for  30-blt  burst. 

3  Simplified  block  diagram  of  the  logic  subsystem  of  the  ultrasonic 
pseudorandom  signal  correlation  system. 

4  Typical  experimental  setup  for  ultrasonic  pseudorandom  signal 
correlation  system. 

5  Impulse  response  of  two  transducers  in  a  water  bath,  (a)  Response 
obtained  using  a  Panametrics  pulser/ receiver,  (b)  Response  obtained 
using  the  ultrasonic  pseudorandom  signal  correlation  system. 

6  Correlator  output  for  burst  length  of  100-bits  (20  psec) .  (a)  Out¬ 

put  with  10  repetitions  per  correlation  point,  (b)  Output  with  100 
repetitions  per  point  showing  the  improvement  in  signal-to-sidelobe 
ratio. 
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Plot  illustrating  the  effect  of  varying  the  nuaber  of  repetitions 

per  correlation  point  upon  the  maximum  side lobe  magnitude  (relative 

\ 

to  the  correlation  peak)  for  a  fixed  burst  length  of  375-bits. 
Calculated  points  were  obtained  from  a  computer  simulation  of  the 
UPSC  system  for  an  autocorrelation  (no  transducer  or  water  effects 
are  included  in  the  simulation) .  The  solid  line  is  a  linear  fit 
to  the  calculated  data.  The  break  in  slope  occurs  when  the  product 
of  the  burst  length  (in  bits)  and  the  number  of  repetitions  per 
correlation  point  exceeds  the  m-sequence  length  (65  535-blts) . 

Example  of  signal-to-noise  ratio  improvement  obtainable  with  the 
UPSC  system.  The  received  signal  is  the  first  echo  through  a  highly 
attenuative  specimen,  (a)  Signal  from  Panametrics  Model  5052  PR 
pulser/ receiver  using  maximum  "energy"  and  maximum  receiver  gain, 

(b)  Signal  from  UPSC  system  using  the  receiver  section  of  the  same 
Pauametrics  unit  (with  the  same  gain)  as  a  preamplifier. 

Illustration  of  the  degradation  in  system  range  resolution  which 
occurs  when  the  reference-channel  bandwidth  is  reduced  by  the 
addition  of  an  acoustic  delay  path  in  the  reference  channel. 

(a)  Correlator  output  with  digitally  delayed  reference  channel 
(normal  UPSC  system  operation) ,  (b)  Correlator  output  with  fixed 
acoustic  delay  in  reference  channel.  This  delay  was  produced 
by  two  transducers  in  a  water  bath. 
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A  PSEUDORANDOM  BINARY  NOISE  NDE  ULTRASONIC  CORRELATION  SYSTEM 
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Abstract 


A'  pseudorandom  binary  noise  NDE  ultrasonic 
correlation  system  has  been  designed  which  retains 
the  signal  to  noise  improvement  capabilities  of  the 
random  signal  correlation  system  while  improving 
data  acquisition  speed  to  near  real  time.  A  pro¬ 
totype  instrument  has  been  constructed  using  a  65K 
bit  maximal  sequence  pseudorandom  noise  generator 
and  its  performance  has  been  evaluated.  One  impor¬ 
tant  consideration  for  the  application  of  this  tech¬ 
nique  is  the  presence  of  time  sidelobes.  To  date, 
using  5  MHz  immersion  NDE  transducers,  the  proto¬ 
type  system  has  shown  that  a  correlation  peak  to 
sidelobe  ratio  of  >  47  dB  can  be  achieved.  It  has 
also  been  shown  that  phase  modulating  a  cw  rf 
carrier  with  the  binary  code  allows  the  system  to 
be  used  to  measure  acoustic  transit  time  changes 
with  ppm  accuracy  or  to  drive  transducers  at  fre¬ 
quencies  above  the  10  MHz  limit  of  the  instrument 
without  the  moduiation.  Additional  details  of  sys¬ 
tem  performance  and  possible  applications  are  dis¬ 
cussed  . 

1.  Introduction 

Large  improvements  in  the  signal  to  noise  ra¬ 
tio  of  ultrasonic  signals  from  targets  typically 
encountered  in  Nondestructive  Evaluation  (NDE)  have 
been  achieved  in  the  laboratory  using  a  random  sig¬ 
nal  correlation  system1.  The  random  signal  system 
achieves  this  improvement  by  driving  the  trans¬ 
ducer  with  a  long  burst  of  white  noise  and  corre¬ 
lating  received  signals  with  a  replica  of  the  drive 
burst.  The  average  transmitted  energy  can  be  made 
much  greater  with  the  random  signal  system  than 
with  short  pulse  excited  conventional  systems  where 
electrical  breakdown  in  the  transducer  limits  the 
transmitted  energy.  The  bandwidth  limiting  which 
occurs  in  the  correlation  receiver  also  contributes 
significantly  to  the  enhanced  sensitivity. 

With  this  system,  the  only  way  a  replica  of 
the  transmitted  signal  can  be  derived  is  by  split¬ 
ting  the  signal  from  the  noise  source.  This  refer¬ 
ence  is  then  delayed  with  an  ultrasonic  delay  line. 
The  delay  line  consists  of  a  pair  of  ultrasonic 
transducers  in  a  water  bath  whose  separation  is 
mechanically  varied.  A  sample  In  the  signal  chan- 
nol  Is  scanned  by  varying  the  reference  delay  over 
1978  Ultrasonics  Symposium  I’roceedi ngs 
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the  range  in  which  delayed  signals  from  targets  in 
the  sample  are  expected  to  occur.  Present  stgnaL 
acquisition  times  are  greater  than  can  be  tolerated 
for  field  NDE  use,  due  to  a  large  degree  to  the 
slowness  of  the  mechanical  scan. 

This  paper  will  describe  a  system  whereby  the 
requirement  for  a  mechanically  scanned  delay  line 
can  be  eliminated  and  the^scanning  can  be  accom¬ 
plished  in  near  real  time  .  The  new  system  sub¬ 
stitutes  a  pair  of  pseudorandom  binary  noise  (PBN) 
sources  which  produce  identical  output  sequences 
for  the  white  noise  source  used  in  the  random  sig¬ 
nal  system.  Many  of  the  features  of  the  random 
signal  system  have  been  retained  and  the  range 
resolution  of  the  new  system  has  been  found  to  be 
superior  to  that  of  a  random  signal  system.  A  pro¬ 
totype  instrument  has  been  constructed  and  its 
capabilities  will  be  described. 

2.  System  Description 

A  block  diagram  of  the  system  is  shown  in 
Figure  1.  The  primary  difference  from  the  random 
signal  system  is  the  pair  of  PBN  generators.  Each 
generator  consists  of  a  sixteen  stage  shift  regis¬ 
ter  with  exclusive  NOR  feedback  chosen  to  produce 
a  65,535  bit  maximal  length  binary  sequence-*. 

Figure  2  shows  the  first  100  bits  of  the  sequence 
produced  by  the  generators.  The  control  logic 
determines  the  length  of  the  noise  burst,  the  re¬ 
lative  time  delay  between  the  reference  burst  and 
the  signal  channel  burst  and  the  rate  at  which  this 
delay  is  varied  for  scanning  of  the  sample.  The 
reference  delay  is  stepped  in  increments  equal  to 
single  periods  of  the  master  oscillator,  resulting 
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Figure  1.  System  block  diagram 
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Figure  2.  First  100  hits  of  the  pseudorandom 
binary  m-sequence  used  with  the 
present  instrument  . 


in  a  sampled  system.  To  insure  that  the  Nyquist 
criterion  for  adequate  sampling  is  fulfilled,  the 
master  oscillator  frequency  is  divided  by  two  for 
the  timing  of  the  PBN  generators,  making  a  single 
bit  time  of  the  sequence  twice  the  reference  delay 
step. 

The  correlation  receiver  uses  the  multiplier 
function  of  a  Phillips  Model  3265  oscilloscope  to 
multiply  the  signal  and  reference  channels  and  a 
low  pass  filter  on  the  scope  output  to  perform  the 
simulated  integration.  The  output  is  displayed  on 
a  storage  oscilloscope. 

Figure  1  illustrates  the  operation  of  the 
system.  In  the  signal  channel,  the  burst  length 
and  cycle  period  are  chosen  to  allow  returning  tar¬ 
get  signals  of  interest  to  be  fully  separated  from 
the  transmitted  signals.  The  reference  burst, 
which  is  identical  to  the  signal  burst,  is  initially 
delayed  slightly  less  than  the  first  return  signal. 
The  delay  is  then  incremented  in  steps  to  a  point 
past  the  delay  time  of  any  signal  of  interest  in 
the  signal  channel.  Each  time  the  reference  delay 
equals  a  signal  delay,  maximum  correlation  occurs. 

It  may  be  noted  from  the  figure  that  each  trans¬ 
mitted  burst  differs  from  the  previous  one.  In 
tact*  if  the  first  burst  represents  N  bits  of  the 
overall  sequence,  each  succeeding  burst  consists  of 
the  next  N  bits  of  the  sequence  beginning  at  the 
point  where  the  previous  burst  ends. 

There  is  one  feature  of  this  system  which  we 
believe  is  unique.  This  is  the  feature  which  al¬ 
lows  the  reference  delay  to  be  held  fixed  for  a 
selected  number  of  repetitions  at  each  step  to  im¬ 
prove  the  signal  to  noise  and  reduce  range  side- 
lobes.  Since  single  transducer  operation  is  a  de¬ 
sirable  option,  the  maximum  duty  cycle  is  limited 
to  less  than  50%  due  to  the  requirement  that  the 
transmitted  be  off  during  the  signal  return  time. 

I'he  maximum  burst  length  is  also  limited  to  less 
than  the  round  trip  transit  time  between  trans¬ 
ducer  and  target.  For  reasonable  path  lengths, 
this  restricts  the  hurst  length  to  a  maximum  on  the 
order  of  25-50  usee.  These  constraints  led  to  the 
design  decision  to  approximate  a  long  sequence 
which  has  good  correlation  properties  with  a  series 
of  short  bursts,  each  of  which  is  a  different  seg¬ 
ment  of  the  overall  sequence.  This  fixing  of  the 
reference  delay  for  many  repetitions  tends  to  aver¬ 
age  out  side  lobes  since  each  burst  produces  a 
different  side lobe  pattern. 


The  logic  speed  in  the  present  prototype  in¬ 
strument  limits  the  clock  rate  of  the  pseudonoise 
generators  to  approximately  13  MHz.  The  PBN  spec- 
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Figure  3.  Simplified  system  timing  diagram 
to  illustrate  operation. 

trum  exhibits  nulls  at  integer  multiples  of  the 
clock  frequencyV  This  allows  use  of  transducers 
with  center  frequencies  up  to  %10  MHz.  Using 
state  of  the  art  components,  it  should  he  possible 
to  run  the  logic  fast  enough  to  drive  transducers 
with  center  frequencies  up  to  50  MHz  directly. 
However,  there  is  a  well-known  radar  technique 
known  as  phase  shift  keying  w'ich  would  allow  high¬ 
er  frequency  transducers  to  be  driven  even  with  our 
preserc  slow  speed  prototype  Instrument.  The  tech¬ 
nique  consists  of  using  the  PBN  sequence  to  phase 
modulate  a  CW  carrier,  allowing  the  matching  of  the 
spectrum  center  frequency  to  the  frequency  of  the 
transducer . 

3.  System  Operation 

The  system  was  initially  tested  using  a  pair 
of  5  MHz  transducers  separated  by  approximately 
10  cm.  in  a  water  tank  to  generate  and  receive  the 
acoustic  signals.  The  purpose  of  these  tests  was 
to  determine  the  level  of  the  range  sidelobes  and 
the  resolution  potential  of  the  system. 

The  range  sidelobes  arise  from  the  fact  that 
the  system  operates  with  a  burst  mode.  It  is  known 
that  a  complete  maximal  length  sequence  transmitted 
as  a  burst  exhibits  range  sidelobes  which  have  a 
maximum  amplitude  on  the  order  of  \ /JW  times  the 
correlation  peak  amplitude  where  N  is  the  number  of 
bits  in  the  sequence^.  Since  our  system  uses 
segments  of  a  complete  sequence,  we  have  examined 
the  sldelobe  levels  experimentally  and  with  a  com¬ 
puter  model  to  determine  the  effect  of  the  number 
of  repetitions  per  reference  delay  point  on  the 
sldelobe  levels. 

Figures  U  and  5  show  examples  of  the  data 
which  must  be  analyzed  to  evaluate  system  perfor¬ 
mance.  Figure  4  shows  the  correlation  peak  and  the 
close-in  sidelobes  on  a  linear  and  log  plot  for  a 
short  hurst  (12.5  iisec)  and  a  small  number  (10)  of 
repetitions.  The  peak  sldelobe  levels  for  this 
case  are  on  the  order  of  20  dB  below  the  correla- 
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tii'ii  peak  height.  When  the  burst  length  Is  in¬ 
creased  to  SO  psec  and  the  number  of  repetitions 
to  99,  the  sidelobe  levels  are  reduced  to  appro¬ 
ximately  36  dB  below  the  peak  as  is  shown  in  Fig¬ 
ure  S.  To  determine  in  a  systematic  way  the  effect 
of  varying  the  repetition  number,  a  computer  model 


Figure  4.  Correlation  peak  for 

short  burst  length  and  small 
number  of  repetitions  per 
reference  delay. 
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Figure  3.  Correlation  peak  lor  a  longer 
burst  and  u  typical  number  of 
repetitions  per  reference 
delay . 


of  the  system  was  devised  and  the  autocorrelation 
function  of  the  signal  and  reference  bursts  was 
calculated.  Figure  6  shows  the  results  of  the 
calculation  for  a  burst  length  of  37S  bits  with 
repetitions  between  1  and  2000  per  delay  point. 

The  amplitude  of  the  maximum  sidelobe  is  plotted 
vs  repetitions.  In  addition,  experimental  results 
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Figure  6.  Experimental  and  calculated 
maximum  sidelobe  level  for  a 
375  bit  burst  length  as  a  function 
of  the  number  of  repetitions  at 
each  reference  delay. 


obtained  both  visually  from  the  oscilloscope  and 
with  a  Northern  Model  NS-575  signal  averager  are 
shown  in  the  same  figure.  Agreement  between  the 
calculated  and  experimental  results  is  excellent 
up  to  100  repetitions.  Above  this  point,  where 
the  entire  PBN  sequence  length  is  exceeded  at  each 
correlation  point,  the  computer  model  predicts  a 
more  rapid  decrease  in  sidelobe  level  with  in¬ 
creasing  repetition  number,  while  the  experimental 
values  appear  to  be  following  the  original  slope. 

At  the  present  time  it  is  not  certain  if  this  dis¬ 
crepancy  is  due  to  experimental  errors,  since  the 
sidelobes  are  near  the  system  noise  levels,  or  to 
an  inadequate  computer  model  of  the  system.  The 
present  model  does  not  take  into  account  the  ef¬ 
fect  of  the  transducer  on  the  signal  channel  hurst. 
It  should  be  noted  that  we  have  been  able  to  ac¬ 
hieve  sidelobe  levels  at  least  47  dB  below  the  pc *k 
height.  The  price  one  has  to  pay  for  such  sidelobe 
levels  is  long  signal  acquisition  times,  almost 
2  minutes  in  the  case  of  2000  repetitions,  making 
it  rather  impractical  to  use  such  large  numbers  ot 
repetitions  with  the  present  system.  One  possible 
method  of  improving  signal  acquisition  time  nl  1 
be  to  use  more  than  one  reference  PBN  generator 
and  parallel  process  the  outputs.  It  should  al 
be  noted  that  an  increase  in  pulse  width  with  the 
repetition  number  held  fixed  will  also  reduce  the 
sidelobe  level. 


The  width  of  the  correlation  peak  was  also  in¬ 
vestigated.  It  was  determined  that  this  width  was 
limited  only  by  the  transducer  response.  In  terms 
of  ability  to  resolve  target  signals  in  time,  there 
is  virtually  no  difference  between  this  system  and 
a  conventional  pulse  NDE  system.  A  comparison  was 
also  made  between  this  system  and  the  random  sig¬ 
nal  system.  An  ultrasonic  delay  produced  by  a  se¬ 
cond  pair  of  transducers  matched  to  the  first  pair 
was  added  to  the  reference  channel  to  approximate 
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ihe  random  signal  system.  Figure  7  illustrates  the 
broadening  of  the  correlation  peak  which  occurs 
when  the  reference  channel  is  also  bandwidth  lim¬ 
ited  as  is  the  case  with  Che  random  system. 

The  final  aspect  of  system  operation  which 
was  investigated  was  performance  in  the  phase  shift 
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Figure  7.  A.  Normal  PBN  system  correlation 
peak.  B.  Broadened  peak  when 
reference  channel  is  bandwidth 
l  iini  Led  to  simulate  the  random 
signal  system. 


to  the  signal  channel,  the  phase  of  the  PBN  se¬ 
quence  after  carrier  removal  on  the  receiver  side 
is  sensitive  to  the  phase  of  the  carrier.  With  no 
reflections,  the  phase  of  the  carrier  in  the  re¬ 
ceived  signal  is  determined  by  the  total  delay  time 
in  the  signal  channel  from  the  modulation  point  to 
the  demodulation  point.  Since  the  electronic  de¬ 
lays  are  negligible,  the  acoustic  delay  determines 

-o 


-20 


-40 


6  10  '  20  '  30 ' 40 

FREQUENCY  (MHz) 


keyed  mode  of  operation.  There  are  two  possLble 
ways  to  implement  such  a  system.  In  the  first, 
which  was  chosen,  a  cw  carrier  is  phase  modulated 
by  the  PBN  sequence  us  in*;  a  double  balanced  mixer 
and  the  resulting  signal  drives  the  transmitting 
transducer.  The  carrier  is  then  removed  from  the 
receiving  transducer  output  by  heterodyning  with  a 
signal  split  off  from  the  cw  carrier  signal  gener¬ 
ator.  The  second  way  is  to  produce  phase  modulated 
signals  in  both  signal  and  reference  channels**. 

Such  a  set-uj>  would  require  a  very  broad  bandwidth 
correlation  receiver  if  it  were  desired  to  operate 
.•t  frequencies  above  20  MHz. 

Figure  HA  shows  the  spectrum  of  an  ordinary 
PBN  sequence  operated  at  a  10  MHz  clock  frequency 
and  illustrates  the  nulls  at  multiples  of  the  clock 
frequency.  Figure  8B  illustrates  the  manner  in 
which  phase  shift  keying  centers  the  spectrum  at 
the  carrier  frequency,  allowing  high  frequency 
transducers  to  be  driven  by  a  low  frequency  PBN 
sequence . 

There  are  several  points  to  be  noted  with  re¬ 
gard  to  operation  of  the  system  in  the  phase  shift 
keyed  mode.  When  the  modulation  is  confined  solely 
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Figure  8.  Frequency  spectrum  of  PBN  se¬ 
quence  with  a  10  MHz  clock  fre¬ 
quency.  A.  Ordinary  operation., 

B.  When  PBN  phase  modulates  a 
20  MHz  cw  carrier. 

this  phase  and  the  phase  of  the  recovered  PBN  se¬ 
quence  since  heterodyning  signals  of  the  same  fre¬ 
quency  produces  an  output  proportional  to  the  phase 
difference.  It  is  possible,  depending  on  the  de¬ 
lay,  to  obtain  a  positive  or  negative  correlation 
peak,  or  if  the  acoustic  carrier  signal  and  the 
oscillator  signal  are  at  quadrature,  a  null  output. 
S'jch  a  response  is  analogous  to  that  obtained  when 
using  a  laboratory  technique  known  as  the  phase 
comparison  method  to  measure  small  (M  part  in 
10^)  changes  in  ultrasonit  velocity^.  With  the 
phase  comparison  method,  the  transmitting  trans¬ 
ducer  is  excited  with  a  tone  hurst  and  the  phase 
of  the  received  echoes  is  compared  to  that  of  a  co¬ 
herent  cw  reference.  In  the  laboratory,  tempera¬ 
ture  or  some  other  variable  is  changed  to  produce 
acoustic  transit  time  changes  (caused  bv  sound  ve- 


120 


AMPLITUDE  (dB)  AMPLITUDE  (dB) 


iocity  or  path  length  changes).  The  frequency  of 
ihe  carrier  is  then  varied  to  hold  the  output  phase 
at  null.  In  this  case 


Af_  _  AT  Av  _  A £ 
f  T  "  v  t 


where  f  is  carrier  frequency,  T  is  total  acoustic  delay 
time,  v  is  sound  velocity,  i  is  the  acoustic  path 
length  and  Af  is  the  frequency  change  needed  to  main¬ 
tain  the  null  condition.  Electrical  delays  are  assumed 
neglible  as  are  phase  variations  produced  by  the  trans¬ 
ducer  as  a  function  of  frequency. 

The  analogy  between  the  outputs  of  the  phase  shift 
keyed  PBN  system  and  the  phase  comparison  ultrasonic 
velocity  measuring  method  indicates  that  it  should  be 
possible  to  measure  small  changes  in  velocity  or  thick¬ 
ness  by  recording  the  change  in  carrier  frequency  nec¬ 
essary  to  maintain  the  null  during  a  scan  of  a  part. 
This  was  tested  experimentally  by  placing  a  pair  of 
9  MHz  transducers  10  cm.  apart  in  a  water  tank.  The 
output  was  initially  set  to  null  and  the  spacing  was 
then  changed  in  0.025  mm  steps  with  the  carrier  fre¬ 
quency  varied  to  return  the  output  to  null  at  each 
step.  After  50  steps  it  was  found  that  the  maximum 
error  in  the  measured  if  for  a  single  step  was  21%, 
with  the  average  value  for  all  points  being  within 
3.5%  of  the  expected  value.  The  largest  source  of 
error  was  vibration  in  the  water  tank  which  produced 
position  changes.  The  error  in  the  average  was  prob¬ 
ably  due  to  error  in  the  mechanical  step  distance 
and  slight  frequency  variation  of  the  transducer 
response.  The  experiment  itself  was  not  set  up  in 
an  optimal  fashion,  but  it  indicated  that  it  would 
not  be  unreasonable  to  expect  resolution  of  10  ppm 
changes  in  velocity  or  dimension  with  this  system. 


It  should  also  be  noted  that  the  other  way  to 
implement  the  phase  shift  keyed  system  which  is  de¬ 
scribed  in  Ref.  6  would  show  the  same  phase  sensi¬ 
tivity  unless  the  carrier  and  both  signal  and  refer¬ 
ence  tPBN  sequences  are  phase  locked  to  ensure  that  the 
burstjs  are  exact  replicas  of  one  another. 

4.  Conclusion 

r 

#The  primary  advantage  of  the  PBN  system  over 
conventional  NDE  ultrasonic  systems  is  the  increased 
sensitivity  of  the  system.  If  the  sample  environment 
is  sufficiently  clean,  i.e.  there  are  no  large 
scatterers,  it  is  relatively  easy  to  obtain  signal 
noistj  Improvements  of  20  to  40  dB  compared  to  the 
conventional  systems.  The  presence  of  spurious 
scatterers  such  as  grain  boundaries,  etc.  would  ob¬ 
viously  produce  interfering  signals  which  would  reduce 
or  eliminate  the  Improvement. 


||  NDE  application  areas  would  principally  be  those 
inspections  where  the  attenuation  is  large  and/or 
deep  penetration  is  required.  Additional  applica¬ 
tions  would  be  to  Improve  the  signal  to  noise  ratio 
for,|,apectral  analysis  of  the  signal  and  in  imaging 
systems  where  the  digital  nature  of  the  system  should 
provide  great  flexibility.  The  phase  shift  keyed  mode 
of  operation  provides  an  additional  sensitivity  to 
small  dimensional  and  velocity  variations  as  well 
as  permitting  the  extension  of  flaw  detection  and 


acoustic  attenuation  measurasients  to  higher  frequen¬ 
cies  using  the  same  basic  Instrument. 
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LOGIC  DESCRIPTION  OF  THE  CONTROL  LOGIC  OF  THE  ULTRASONIC 
PSEUDORANDOM-SIGNAL-CORRELATION  SYSTEM  USING  AN  ASM  CHART 


An  algorithmic  state  machine  (ASM)  Is  a  concise  and  very  useful  means  of 
modeling  a  logic  system  [1].  Such  a  machine  has  a  number  of  predefined 
states  and,  in  general,  the  next  state  of  the  machine  depends  on  its  pre¬ 
sent  state  and  certain  inputs  or  qualifiers.  In  a  given  state,  the  machine 
may  produce  certain  outputs.  If  these  outputs  are  a  function  of  that  state 
alone  they  are  called  state  outputs.  If  the  outputs  depend  on  the  state 
and  some  other  condition(s) ,  they  are  conditional  outputs.  Often  the  states 
of  a  machine  are  defined  logically  by  a  state  counter;  when  the  counter  is 
at  count  S,  the  machine  is  said  to  be  in  state  S.  The  logic  which  controls 
the  counting  sequence  is  thus  the  state  control  logic  for  the  machine. 

An  ASM  is  represented  graphically  by  an  ASM  chart.  The  chart  uses  the  same 
rectangular  box,  diamond  and  "rounded  box"  symbols  which  are  used  in  flow 
charts  and,  in  fact,  bears  a  superficial  resemblance  to  a  flow  chart.  How¬ 
ever,  in  an  ASM  chart  the  operations  described  within  a  rectangular  box  and 
all  branching  associated  with  that  box  are  executed  simultaneously,  not 
sequentially  as  in  a  flow  chart.  Each  rectangle  or  state  box  and  its  asso¬ 
ciated  branching  and  conditional  output  symbols  represent  a  single  machine 
state.  The  diamond  symbol  or  decision  box  is  used  to  represent  conditional 
branching  and  the  symbol  with  straight  parallel  top  and  bottom  and  semi¬ 
circular  ends  is  a  conditional  output  box.  State  outputs  are  written  in 
the  state  box  which  corresponds  to  the  state  in  which  they  occur.  The  state 
number  associated  with  a  particular  state  box  is  written  outside  the  box 


122 


•t  Che  upper  right  hand  corner.  A  Boolean  expression  representing  a  condi¬ 
tion  or  qualifier  la  written  in  the  condition  box.  The  two  exit  paths  of 
a  condition  box  are  labeled  H  and  L.  The  H  path  is  taken  if  the  qualifier 
Is  HIGH,  the  L  path  If  the  qualifier  is  LOW.  A  conditional  output  box 
placed  In  a  conditional  exit  path  Indicates  that  those  outputs  written  In 
the  box  will  be  asserted  only  when  that  path  is  taken. 

The  ASM  chart  for  the  control  logic  of  the  pseudorandom  signal  correlation 
system  in  shown  in  Fig.  1  .  The  "rest  state"  for  the  machine  is  state  zero 
(000).  In  this  state  the  operator  may  set  the  following  parameters:  burst 
length,  cycle  period  -  the  time  between  transmitted  noise  bursts,  initial 
correlation  delay  -  the  initial  starting  time  for  the  reference  noise  burst, 
maximum  correlation  delay  -  the  final  starting  time  for  the  reference  noise 
burst,  and  repetitions  per  point  -  the  nunfcer  of  noise  bursts  transmitted 
before  the  correlation  delay  is  increased.  The  values  for  these  parameters 
are  held  in  storage  registers. 

The  logic  term  (ER-L)  is  asserted  in  state  zero  to  enable  the  loading  of  th'- 
storage  registers.  The  terms  LNC-L,  LQC-L,  LPC-L,  LP*C-L,  and  LIC-L  are 
also  asserted  in  state  zero.  These  terms  cause  the  N,  Q,  P*,  and  I  counters 
respectively,  to  be  loaded  or  preset  to  the  nuafcer  stored  in  the  registers. 
The  P  counter  is  loaded  from  the  P*  counter  which  is  itself  loaded  from  a 
storage  register.  0A-L  is  used  to  set  the  two  pseudorandom  binary  sequence 
generators  to  zero.  Switching  the  RUN/HALT  switch  to  RUN  sets  the  GO  flip- 
flop.  GO-H  causes  the  ASM  to  advance  to  state  1  (001)  on  the  following 
clock  transition.  Note  that  the  state  counter  for  this  ASM  (as  well  as  the 
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Figure  1. 
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I,  N,  P,  P*,  and  Q  counters)  is  clocked  by  the  system  20  l«z  clock.  This 
■sans  that  all  state  transitions  will  be  synchronised  with  this  clock.  Here¬ 
after,  unless  otherwise  stated,  the  term  "clock"  refers  to  the  20  MHz  system 
clock  and  "one  clock  tine"  is  one  period  of  this  clock. 

In  state  1,  (RUN  1-H)  enables  pseudorandom  binary  sequence  generator  number  1 
(PRBS  #1).  PRBS  #1  serial  output  is  "anded"  with  (RUN  1-H)  and,  after  level 
shifting  and  amplification,  becomes  the  transmitted  noise  burst.  In  state  1, 
none  of  the  counter  "load  terms"  (LNC-L ,  LQ-L,  etc.)  are  asserted.  The  con¬ 
trol  logic  for  the  N,  P,  and  Q  counters  is  such  that  these  counters  count 
down  when  their  load  terms  are  not  asserted.  Thus,  in  state  1  PRBS  #1  is 
nmn^ng  and  the  N,  P,  and  Q  counters  are  counting  down  from  their  preset 
counts.  When  the  N  counter,  which  controls  the  noise  burst  length,  counts 
down  to  the  count  of  2,  NE2-H  is  asserted  and  one  clock  time  later  the  machine 
enters  state  2  (010) . 

State  2  is  entered  with  the  N-counter  count  equal  to  1.  In  state  2  LNC-L  is 
asserted  so  that  the  zero  count  of  N  counter  never  occurs,  but  the  counter  is 
agalq. preset  to  the  number  contained  in  the  "N  register"  NR  (burst  length). 

In  state  2  (RUN  1-H)  is  no  longer  asserted,  therefore,  PRBS  #1  stops  running, 
ending  the  transmitted  noise  burst.  The  gating  term  for  PRBS  #1  is  also 
forced  LOW  in  state  2.  This  ensures  that  the  gated  noise  output  is  LOW  even 
if  PRBS  #1  was  stopped  with  a  ONE  or  HIGH  output  condition  (PRBS  #2  is  gated 
in  a  similar  manner) .  In  state  2  the  N  counter  is  held  at  its  preset  value, 
but  the  P  and  Q  counters  are  still  counting.  When  the  P  counter,  which 
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controls  the  starting  tine  of  the  reference  noise  burst (  reaches  the  count 
of  2,  PE2-H  Is  asserted  and  state  3  (011)  is  entered  one  clock  tine  later. 

In  state  3,  LPC-L  causes  the  P  counter  to  preset  to  the  nunber  held  in  the 
p*  counter  and  to  hold  this  count.  RUN  2-H  enables  the  reference  pseudo- 
randon  binary  sequence  generator  PRBS  *2  and,  since  LNC-L  is  no  longer  asserted, 
the  N-counter  starts  counting  down.  The  ASM  remains  in  state  3  for  one  clock 

time  and  then  enters  state  4  (100) .  (Due  to  a  logic  design  change,  state  3 

Is  now  actually  not  needed  since  the  logic  which  was  previously  unique  to 
state  3  is  now  implemented  elsewhere.) 

i- 

In  state  4  (100) ,  PRBS  #1  is  still  running  and  the  N  and  Q  counters  still 
counting  down.  Just  as  for  the  transmitted  burst,  the  N  counter  controls 
the  reference  burst  length.  When  NE2-H  signifies  that  the  N  counter  has 
counted  down  to  2,  the  state  changes  to  5  (101) .  RUN  2-H  is  no  longer 
asserted  so  the  reference  burst  ends,  the  N  counter  is  again  loaded  with 
the  contents  of  register  NR,  and  neither  the  N  nor  the  P  counter  is  counting. 

Now  that  the  reference  burst  has  ended,  both  PRBS  #1  and  PRBS  #2  should  hold 

the  same  number.  The  parallel  outputs  of  the  shift  registers  used  in  PRBS  #1 
and  PRBS  #2  are  fed  to  a  comparator.  If  the  contents  of  the  registers  are 
the  same,  AEQB-H  is  asserted.  After  a  one  clock  time  stay  in  state  5  to 
allow  for  comparator  propagation  delays,  the  ASM  advances  to  state  6  (110). 

In  state  6,  if  AEQB-H  is  not  asserted  the  ASM  enters  state  7  (111)  which  is 
the  "error  state".  In  state  7,  an  error  indicator  on  the  instrument  front 
panel  is  lighted  and  KGO-H  resets  the  00  flip-flop.  When  the  RESET  button 
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la  depressed,  RST-L  causes  the  ASM  to  return  to  state  zero.  If  AEQB  Is 
asserted  in  state  6,  state  7  Is  not  entered.  Instead,  the  ASM  waits  for  '• 
QE2-H  to  indicate  that  the  Q  counter,  which  controls  the  time  between 
transmitted  noise  bursts,  has  counted  down  to  2.  When  QE2-H  Is  asserted, 
IE0-H,  which  is  HIGH  when  the  I  counter  has  counted  down  to  zero,  determines 
which  of  two  sets  of  conditional  outputs  will  be  asserted.  The  I  counter 
was  preset. (In  state  zero)  to  the  nunfcer  of  "repetitions  per  correlation 
point".  (Note  that  one  "repetition"  has  already  taken  place  before  the 
I-counter  is  checked  so  the  counter  should  be  preset  to  a  number  which  is 
one  less  than  the  desired  number  of  repetitions.)  If  IE0-H  is  not  asserted, 
the  I-counter  does  not  contain  zero  and  the  L  branch  of  the  decision  box 
Is  taken.  In  this  case,  EIC-L  is  asserted  enabling  the  I-counter  to  count 
down.  If  the  I-counter  contains  zero,  IEO-H  is  asserted,  the  H  branch  is 
taken  and  LIC-L  reloads  the  I  counter  with  the  contents  of  the  I  register 
IK.  Also,  EP*C-L  enables  the  P*  counter  to  count  up.  Since  QE2-His  HIGH 
for  only  one  clock  time,  EIC-L  or  EP*C-L  can  be  asserted  for  only  one  clock 
time  which  assures  that  the  I  or  P*  counter  will  change  by  only  one  count. 
Since  LPC-L  Is  asserted  for  the  duration  of  state  6,  the  P  counter  will 
have. its  required  one  clock  time  to  preset  to  the  updated  contents  of  the 
P*  counter  before  QE1-H  causes  the  ASM  to  exit  state  6. 

The  term  PMAX-H  is  asserted  when  the  contents  of  the  P*  counter  equals  the 
contents  of  the  register  PMR.  PMR  Holds  the  final  or  maximum  value  of  the 
correlation  delay.  When  QE1-H  is  asserted,  LQC-L  enables  reloading  of  the 
Q  counter.  If  PMAX-H  is  not  asserted  and  the  GO  flip-flop  is  still  set 
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(Indicating  that  the  RUN/HALT  switch  la  still  In  RUN) ,  state  1  (#01)  Is 
re-entered  and  RUN  1-H  starts  a  new  transmitted  burst.  If  the  P*  counter 
was  incremented  In  state  6,  PE2-H  will  be  asserted  one  clock  time  later 
than  It  was  for  the  previous  state.  This  means  that  the  exit  from  state 
2(  hence  the  start  of  the  reference  burst,  will  occur  one  clock  time  later. 
In  summary,  when  the  I  counter  has  counted  the  required  number  of  repeti¬ 
tions  for  a  given  preset  value  of  the  P  counter  (correlation  delay) ,  the 
P*  counter — from  which  the  P  counter  is  loaded —  Is  incremented  by  one 
count  thus  increasing  the  delay  by  one  clock  time. 

If  PMAX-H  were  asserted  in  state  6,  this  would  indicate  that  the  maximum 
correlation  delay  had  been  reached.  In  this  case,  KGO-H  resets  the  GO 
flip-flop  and  the  ASM  returns  to  state  zero.  In  state  zero  all  the  counters 
are  reloaded,  PRBS  #1  and  PRBS  #2  set  to  zero  and,  if  the  RUN/HALT  switch 
is  in  RUN,  the  GO  flip-flop  is  set  and  the  ASM  goes  to  state  1. 
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MNEMONICS  INFORMATION 


NOTE:  Each  mnemonic  consists  of  a  character  string  followed  by  a 

hyphen  and  then  an  H  or  an  L.  The  H  Implies  that  the  logic 
signal  is  asserted  when  HIGH:  the  L  means  the  signal  is 
asserted  when  LOW. 


Mnemonics  Used  on  ASM  Diagram 


Mnemonic _ Description 


AEQB-H 

The  contents  of  the  two  pseudorandom-noise- 
generator  shift  registers  are  numerically 
equal. 

EIG-L 

Enable  the  I  counter. 

EP*C-L 

Enable  the  P*  counter. 

(ER-L) 

Enable  loading  of  data  registers  (actually 
Implemented  as  a  buffered  [inverted]  decoded 
state-count  zero) . 

GO-H 

GO  flip-flop  is  set. 

IE0-H 

I  counter  count  equals  zero. 

KGO-H 

K  input  to  GO  flip-flop. 

LIC-L 

Load  I  counter. 

LNC-L 

Load  N  counter. 

LPC-L 

Load  P  counter. 

LP*C-L 

Load  P*  counter. 

LQC-L 

Load  Q  counter. 

NE2-H 

N  counter  count  equals  2. 

PE2-H 

P  counter  count  equals  2. 

PMAX-H 

P*  counter  count  equals  contents  of  "PM" 
register. 
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Mnemonics  Used  on  ASM  Diagram  (Cont'd) 


Description 


QE1-H 


Q  counter  count  equals  1. 


RST-L 


Reset.  Clears  GO  flip-flop.  Also  a  qualifier 
for  exit  from  ERROR  state  (111)  to  return  to 
state  000. 


RUN1-H 


RUN  command  for  pseudorandom-noise  generator 
number  1. 


RUN2-H 


RUN  command  for  pseudorandom-noise  generator 
nunber  2. 
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MNEMONICS  INFORMATION 


Registers 


Mnemonics  Description 

IN  I  register.  Holds  the  count  to  which  the  I 

counter  Is  preset.  This  Is  the  number  of  noise 
bursts  to  be  transmitted  for  each  correlation 
point. 

NR  N  register.  Holds  the  count  to  which  the  N 

counter  is  preset.  This  represents  the  noise 
burst  length. 

PR  P  register.  Holds  the  initial  value  of  the  count 

to  which  the  counters  PC  and  P*C  are  preset. 

This  is  the  initial  starting  time  of  the  reference 
burst  or  the  initial  correlation  delay. 

PMR  P-Max  register.  Holds  the  final  value  of  the 

correlation  delay. 

QR  Q  register.  Holds  the  count  to  which  QC  is 

preset  which  is  the  time  between  transmitted 
noise  bursts. 
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MNEMONICS  INFORMATION 


Counters 


Mnemonics _ Description _ 

IC  I  counter.  Presettable  down  counter.  Determines 

number  of  noise  bursts  generated  per  correlation 
point.  "P*C"  counter  is  incremented  by  one  count 
each  time  I  has  reached  zero  and  the  P*C  count 
determines  the  starting  value  for  the  "P"  counter. 

NC  N  counter.  Presettable  down  counter.  Determines 

noise  burst  length.  Counts  upon  negation  of  load 
signal. 

PC  P  counter.  Presettable  down  counter.  This 

counter  is  preset  with  the  current  value  of  "P*" 
counter.  The  P  counter  determines  the  start 
time  of  the  reference  noise  burst  or  the  "corre¬ 
lation  delay". 

p*C  P*  counter.  Presettable  up  counter.  Holds  the 

current  value  of  the  correlation  delay  (the  start 
time  of  the  reference  noise  burst).  This  counter 
is  incremented  by  one  to  increase  the  correlation 
delay . 

QC  Q  counter.  Presettable  down  counter.  Determines 

time  between  transmitted  noise  bursts.  Counts 
upon  negation  of  load  signal. 
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MNEMONICS  INFORMATION 


Generators 

Mnemonics  Description 

PRBS  #1  Pseudorandom  binary  sequence  generator  #1 

PRBS  #2  Pseudorandom  binary  sequence  generator  #2 

(reference) . 


133 


-4  4  4  >  '■» 


C.  R.  Clare,  Designing  Logic  Systems  t 
IT:  McGraw-Hill  Book  Company,  1973. 


13*» 


